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Abstract 

In this article we use the gauge/gravity duality to investigate various properties of strongly cou- 
pled gauge theories, which we interpret as models for the quark-gluon plasma (QGP). In particular, 
we use variants of the D3/D7 setup as an implementation of the top-down approach of connecting 
string theory with phenomenologically relevant gauge theories. 

We focus on the effects of finite temperature and finite density on fundamental matter in the 
holographic quark-gluon plasma, which we model as the M = 2 hypermultiplet in addition to the 
A/" = 4 gauge multiplet of supersymmetric Yang-Mills theory. 

We use a setup in which we can describe the holographic plasma at finite temperature and either 
baryon or isospin density and investigate the properties of the system from three different viewpoints. 

(i) We study meson spectra. Our observations at finite temperature and particle density are in 
qualitative agreement with phenomenological models and experimental observations. They agree 
with previous publications in the according limits. 

(ii) We study the temperature and density dependence of transport properties of fundamental 
matter in the QGP. In particular, we obtain difi'usion coefficients. Furthermore, in a kinetic model we 
estimate the effects of the coupling strength on meson diffusion and therewith equilibration processes 
in the QGP. 

(iii) We observe the effects of finite temperature and density on the phase structure of fundamental 
matter in the holographic QGP. We trace out the phase transition lines of different phases in the 
phase diagram. 

1 Introduction 

The entire content of matter and radiation in the universe is a manifestation of the energy unleashed 
in some unknov^rn process which we commonly refer to as the big bang. This event is thought of as 
the moment of the creation of matter, space and time — the universe. From that moment on energy 
existed in various manifestations. At an order of magnitude of 10"^'^ seconds after the big bang, quarks 
formed. Today we experimentally detect these particles together with leptons and the force mediating 
gauge bosons as the fundamental constituents of all visible matter. The interaction of these particles 
is described incredibly accurately within two different theories. Processes taking place at energy levels 
below the TeV scale involving electromagnetic, weak and strong interaction are described accurately 
by the standard model of particle physics, although the strong force is hard to exploit theoretically at 
low energies for mathematical reasons. The fourth of the known forces, gravity, is described within the 
separate framework of general relativity. 

While many aspects of the particles that make up our world are well understood, others remain a 
mystery. Among the latter is the behavior of matter under conditions that must have existed shortly 
after the big bang. The earliest period of the universe that can either be described by theoretical models 
and numerical simulations or probed by experiments is ranging from about 10"'^^ seconds after the 
big bang when quarks and gluons emerged until approximately 10~^ seconds after the big bang when 
hadronization of quarks set in. During this early phase matter existed in conditions of extremely high 
density and temperature. Under these conditions quarks are not confined and do not form hadrons. 
Instead they are moving independently and interact with each other predominantly via the exchange of 
gluons, which mediate the force of strong interactions. As typical for plasmas, the freely moving quarks 



allow for (color) charge screening. Matter in this phase is therefore referred to as the quark-gluon plasma 
(QGP). 

From the beginning on, the universe expanded and cooled down. After 10 ^ seconds at a critical 
vahic of the temperature of approximatcily IGO 190 MeV the energy dependent coiipling constant of 
the strong force rose to a value that let confinement set in. Eventually quarks combined to bound 
states and formed the hadronic matter that is now composing the galaxies visible in the universe. Only 
very particular regions of the present universe come into consideration for providing conditions extreme 
enough to contain matter in the phase of the quark-gluon plasma. Such barren places are the cores 
of neutron stars. These are remnants of supernova explosions of stars of about 20-30 solar masses. 
Due to the extremely high gravitational pressure, the hadronic matter existing on the planets surface 
in deeper layers is squeezed together to such an extent that electrons and protons combine to neutrons 
(thereby emitting neutrinos). From the surface towards the core of these objects the pressure increases. 
In the inner layers even neutrons are not stable anymore. Instead the quarks and gluons may interact 
individually to appear as the QGP. Other temporary habitats of the quark-gluon plasma seem to exist on 
earth: The experiments conducted at heavy ion colliders are dedicated to monitor the processes occurring 
at collisions of heavy nuclei at energies high enough to produce a fireball of extremely hot and dense 
matter. Such experiments are hosted at the Super Proton Synchrotron (SPS), the Relativistic Heavy 
Ion Collider (RHIC) accelerating gold nuclei, and in future also the Large Hadron Collider (LHC) which 
can be used to accelerate lead nuclei, as well as future SIS experiments at the Facility for Antiproton 
and Ion Research (FAIR) . The state of matter observed at RHIC is a strongly coupled system composed 
of deconfined quarks and gluons, the strongly coupled quark-gluon plasma (sQGP). 

To get an impression of the processes occurring during the first moments after the creation of the 
universe — including the interactions of quarks and gluons which lead to the genesis of the hadronic 
matter that composes our world — it is necessary to understand the properties of the quark-gluon plasma. 
This will eventually allow for deeper insight into the process of hadronization and the phase transition 
from the quark-gluon plasma to the hadronic phase. Further knowledge about the nature of matter may 
also allow for progress in finding the answer to questions about the nature of dark matter and dark 
energy, the majority of the energy content of our universe — by far greater than the contributions visible 
matter can account for. 

Still a manageable theoretical description of the interaction of quarks and gluons in the strongly 
coupled systems observed at experiments is not straightforward, although the standard model contains a 
theory of quarks and gluons, known as quantum chromodynamics (QCD). It is the strong coupling that 
impedes the solution of the equations of motion of QCD at low energies. Analytical answers from QCD 
are obtained from perturbation expansions in the coupling constant, which do not converge at strong 
coupling. Therefore today there is no analytic description of the formation of bound states of quarks 
or the interaction of quarks and gluons in the sQGP from first principles. One successful alternative to 
obtain results at strong coupling is lattice gauge theory, which tries to simulate the dynamics of QCD 
numerically on a number of discrete points in spacetime. While this approach gave answers to numerous 
questions, by nature it cannot produce analytical results which would lead to conceptual insights. It 
moreover approximates spacetime as a coarse grid and so far has to incorporate some simplifications of 
QCD. 

A completely different approach to the quark-gluon plasma can be pursued from a point of view that 
also motivates the work presented in this work. A possible alternative to the description of strongly 
coupled quarks and gluons may be given in terms of string theory. This theory assumes strings to 
be the fundamental degrees of freedom, from which all matter is composed. The different elementary 
particles we know are thought to arise as the different oscillation modes of the strings. Initially, around 
1970, it aimed to explain the relation between spin J and mass m of the resonances found in then 
performed collision experiments, J = ao + a'w?, with a' known as the "Regge slope". The idea was 
to describe the force between quarks as if a string of tension 1/a' holds the particles together. Despite 
modeling the Regge behavior, the theory failed to describe the observed cross sections correctly and was 
successfully displaced by QCD wherever applicable. Nevertheless, the understanding and interpretation 
of string theory evolved to a great extent, especially the possibility to describe quantized gravity attracted 
interest. Today it is the most promising candidate for a unified description of all known forces of nature 
within one single theory. In this sense it can be thought of as a generalization of the successful standard 
model by including a description of gravity. The world of strings appears as a stunningly complex system 
that may give answers to such fundamental questions as the origin of the number of spacetime dimensions 
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we live in, and allow for a formulation of quantum gravity. Still much of the theory has to be understood 
and almost no predictions lie within the reach of experimental verification. 

However, during the past dozen years evidence mounted that indeed there are connections between 
string theory and gauge theories, like QCD. In the mid 1990s, during the so called second string theory 
revolution, it was discovered that string theory not only features strings as degrees of freedom. In 
addition, there are higher dimensional objects, called branes as an allusion to membranes. Branes and 
strings interact with each other. In this way branes influence the degrees of freedom introduced by the 
string oscillations. As the understanding of string theory grew, it was discovered that certain limits of 
string theory contain the degrees of freedom of particular non-Abelian gauge theories. These insights 
heralded a new era of applications of string theory to problems in gauge theory. It began in 1997 with 
Juan Maldacena's discovery of analogies between the classical limit of so called type IIB string theory, 
including branes, and the A/" = 4 supersymmetric Yang-Mills quantum gauge field theory. It is possible 
to establish a one to one mapping between the degrees of freedom of both theories. Maldacena therefore 
speculated that both of them are different descriptions of the same physical reality. The formulation of 
this conjecture is known as the AdS / CFT correspondence [6] . As we will discuss later, this correspondence 
between type IIB string theory in Anti-de Sitter space (AdS) and the A/" = 4 supersymmetric non-Abelian 
conformal field theory (CFT) is especially suitable to describe the strongly coupled regime of gauge 
theories. An astonishing feature of the AdS/CFT correspondence is that it conjectures the equivalence 
of a classical theory of gravity and a quantum field theory. In some aspects this field theory resembles 
the properties of QCD. Moreover, it relates the strongly coupled regime of the quantum field theory to 
the weakly coupled regime of the related gravity theory. One therefore can obtain strong coupling results 
of field theory processes by means of well established perturbative methods on the gravity side. Finally, 
the AdS/CFT correspondence allows to interpret the quantum field theory to be the four dimensional 
representation of processes in string theory, which is defined in ten spacetime dimensions. Because 
of these properties the correspondence is more generally also referred to as gauge/gravity duality and 
is said to realize the holographic principle. So far there is no mathematically rigorous proof for the 
correspondence to hold. Nevertheless, in all cases that allowed for a direct comparison of results from 
both theories, perfect matching was found. 

The AdS/CFT correspondence is considered as one of the most important achievements in theoretical 
physics of the last decades. However, by now the string theory limit which exactly corresponds to QCD 
is not known. Albeit the direct way ahead towards a comprehensive analytical description of strongly 
coupled QCD is not foreseeable, numerous cornerstones where already passed and some junctions and 
connections to the related physical disciplines where found. Examples are deeper insights into the 
connection of black hole physics to thermodynamics, the relation to finite temperature physics, and the 
discovery of quantities like the famous ratio rj/ s oi shear viscosity to entropy density that are universal 
for large classes of theories. The motivation to use the AdS/CFT correspondence to explore the strongly 
coupled quark-gluon plasma therefore is twofold. On the one hand side there is the attempt to provide 
a description of strongly coupled quarks and gluons as a supplement to QCD. In this way string theory 
might contribute to further understanding of gauge theories. On the other hand a phenomenologically 
relevant application of string theory can be used as a benchmark to evaluate the capabilities of string 
theory in describing nature. In this way string theory might benefit from the exploration of new regimes 
of QCD, so far described predominantly by quantum field and lattice gauge theories. The ability to 
produce the sQGP in collision experiments for the first time may allow to check predictions from string 
theory. There is well-founded hope that the quark-gluon plasma can provide a link between string theory 
and experiment. 

We will make use of the AdS/CFT correspondence to investigate strongly coupled systems. The 
models we use for this purpose will be various modifications of the gauge/gravity duality that allow for 
the description of quantum field theories that feature certain aspects known from QCD. Which aspects 
and which parameter regions we can cover with this approach will be pointed out in the introductory 
section on AdS/CFT and in those sections where we introduce the models. It is interesting in its 
own to see how far the correspondence may be extended and which facets of quarks and gluons can 
be modeled at all. However, it is even more fascinating to see that already today some properties of 
phenomenological relevance can be captured by a so called holographic description of the sQGP via the 
AdS/CFT correspondence. Such results allow for comparison with lattice gauge theory and effective 
field theories. A vast number of attempts to apply the correspondence to the dynamics of quarks and 
gluons has been under investigation during the past years. The questions pursued in this work are the 
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following. 

• Can quarks and gluons combine to form hadrons inside the quark-gluon plasma? How does the 
spectrum of bound states of quarks, esp. of mesons, look like inside the sQGP? 

• How do these spectra and the lifetime of mesons depend on temperature and quark density? 

• How do quarks and their bound states move through the plasma? 

• What effects has the strong coupling? 

The answers we obtain are by part of qualitative nature, or can be expected to receive corrections, which 
can be calculated as soon as progress in the field allows to relax some limiting assumptions. Nevertheless, 
it is amazing to see that the gauge/gravity duality can give answers to these questions in terms of a 
minimal number of input parameters. We do not follow the so called "bottom-up" approach, also known 
as AdS/QCD. There the goal would be to find gravity duals to phenomenological gauge theories which 
incorporate certain desired aspects of QCD. Instead we pursue the "top-down" approach. This means 
that we are aware of the fact that the AdS/CFT correspondence is a phenomenon discovered in string 
theory. We try to construct models which are consistent solutions of string theory and observe the 
consequences on the gauge theory side. From this point of view, the results obtained by means of the 
AdS / CFT correspondence can be interpreted as a sign of the predictive power of string theory. 

This article is organized as follows. Section [5] gives a brief introduction to the gauge/gravity duality, 
the extensions which are of relevance for the derivation of our results, and a short discussion of the 
application to QCD and the quark-gluon plasma. However, we will not try to give an introduction 
neither to string theory nor to quantum field theory, supersymmetry or general relativity. Nevertheless, 
these theories are the basis of this work. Especially string theory is on the one hand the basis of this 
work, on the other hand too rich to provide a broad background in detail here. Therefore, we will provide 
the necessary theoretical arguments and details wherever needed in a hopefully adequate manner. The 
remaining sections deal with the answers of the questions mentioned above. Each of them contains a 
brief introductory section and one or more technical sections which lead to results that will be discussed 
at the end of each section. In particular, section [3] deals with meson spectra at finite temperature and 
particle density, and discusses the infiuence of these parameters on the spectra. In section 2] transport 
coefficients of quarks and mesons in the plasma are calculated. Comparison with weak coupling results 
enables us to estimate the effects of strong coupling. In section[S]we examine the lifetime and stability of 
mesons at different temperatures and particle densities and in this way get new insights into the structure 
of the phase diagram of the dual field theory. A summary and discussion of the results is finally given in 
section ini Appendices clarify conventions and notational issues and present some calculations in a more 
detailed form than the main text allows for. 

2 The AdS/CFT correspondence and extensions 

Conventional holograms are able to encode truly three-dimensional information on a two-dimensional 
surface. Analogously, in particle physics and quantum gravity, the equivalence of information contained in 
a theory defined in some lower dimensional space and a different theory on a higher dimensional domain, 
is referred to as the holographic principle. One of the first observations of such kind of holography was 
the discovery that the information captured inside the horizon radius of a black hole, i. e. the entropy 
given by the number of possible microstates, can be described in terms of the horizon surface area alone 
[7] . This observation suggests the existence of a holographic realization of quantum gravity. 

Another observation of a holographically realized connection between gravitational physics and quan- 
tum mechanics was made by Juan Maldacena at the end of the so-called second string theory revolution. 
He then conjectured the equivalence of a supergravity theory in Anti-de Sitter spacetime (AdS) and a 
certain type of conformal field theory (CFT) This discovery triggered an enormous amount of efforts 
to establish the long sought connection between quantum gauge field theories and gravity, which did not 
abate so far. The fact that the theory on the AdSside of the correspondence can be expressed as a low 
energy limit of string theory, which naturally incorporates gravity, is widely interpreted as a support 
of the claim of string theory to offer a formalism which allows for a unified description of all known 
fundamental forces of nature. 
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In this section we briefly review the Maldacena conjecture and some of the extensions invented during 
the last decade. Instead of giving an exhaustive review, we merely draw a hopefully concise and consistent 
sketch of the whole picture. In doing so we emphasize those features of the correspondence that are most 
important for the developments in the subsequent sections. Classical reviews which deal with the subject 
in depth are refs. [8l|9]. For the sake of clarity, we restrict explicit calculations to a minimum here. 



2.1 The original AdS/CFT correspondence 

The AdS/CFT correspondence as it was formulated by Maldacena in 1997 relates two particular theo- 
ries which we introduce subsequently. Afterwards we rephrase the conjecture of Maldacena, before we 
comment on the relaxation of the underlying assumptions, other generalizations and the applicability to 
QCD and the quark-gluon plasma in the following sections. 



2.1.1 J\f — 4 Super- Yang-Mills theory 

One of the two theories related by the AdS / CFT correspondence is Af — A super- Yang-Mills theory in 
four spacetime dimensions of Minkowski topology. It is a supersymmetric quantum field theory with 
SU(A/' = 4)r R-symmetry, which rotates the four supercharges into each other. All fields are arranged in 
one supersymmetry multiplet. The on-shell field content is given by six real spacetime scalars X' with 
i = 1,2, ... ,6, one spacetime vector field A and four two component spin 1/2 left Weyl fermions A" with 
a = 1,2,3,4. Under R-symmetry transformations the six scalars transform as an antisymmetric 6 of 
rank two. The Weyl fermions represent a 4 and the vector field is a singlet. 

With respect to gauge symmetries, all the fields constitute one single multiplet, called the A/" = 4 
gauge multiplet. They transform under the adjoint representation of the gauge symmetry group SU(iV), 
where the integer TV is left as a parameter for now. Later we will interpret it as the number of color 
degrees of freedom. 

The according gauge indices labeling the elements of the gauge symmetry generators T'^ with k — 
1,2,..., N"^ — 1 are suppressed in our notation, e. g. the notation X* for a spacetime scalar is the short 
form of X^^^'' where T'^ as a matrix has elements {T'')™'^ labeled by m, n = 1, 2, . . . , TV. One would 
write out the elements of X"^ as Ar*™„, where i labels the index which transforms under the 6 of the 
R-symmetry and m, n are the indices transforming under SU(A'^) gauge symmetries. We label spacetime 
directions by and v. With this convention and the field strength tensor F — dA + A A A the unique 
Lagrangian ^ reads as 

L = Tr - J^F^.i^f- + ^F^, *F^- - I'X^a^D^K - D^X' + 5,,, CfXa [X' , A,] 

(2.1) 

+ 5vM aab'x" [x\ 'X'] + ^ [X\ X^ ' 

where the trace is performed over the suppressed gauge indices and D is the gauge covariant derivative. 
The symbol Oj denotes the real valued instanton angle, the Cf and Ciab are related to the structure con- 
stants of the R-symmetry group, and there is one dimensionless coupling constant g^i in this Lagrangian. 
The energy dimensions of the operators and constants are given by 

[A] = [X^] = 1, [^^] = l, [5vm] = N-0, (2.2) 

so the action S = /d^x L is scale invariant. In fact J\f — 4 theory is invariant under transformations 
generated by the conformal symmetry group S0(4, 2) = SU(2, 2) composed by Poincare transformations, 
scaling and so-called superconformal transformations as well as under the above mentioned R-symmetry 
group SU(4)p{ . These transformations compose the global symmetry group denoted by PSU(2, 2 | 4). Note 
that these symmetries are realized also in the quantized theory and not broken by anomalies. Moreover, 
the Lagrangian of A/" — 4 super- Yang-Mills theory is unique. In contrast to other supersymmetric theories, 
which allow for different choices of the potential for the superfields, the form of the action is completely 
determined by the demand for renormalizability. 
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There even is a further symmetry. With respect to the complex combination of the couphng g^M and 
the instanton angle 9i given by 

(2.3) 

27r 



:fYM 



the action is invariant under r i~->- t + 1. Generalizing this symmetry, the Montonen-Olive conjecture 
states that the theory is invariant under SL(2,K) transformations acting on the complex coupling t, 
this symmetry is denoted as S-duality. It includes a transformation r i-> — 1/r which indicates that the 
theory describes a duality between strongly and weakly coupled regimes. We will not make use of this 
duality, though. 

The large N limit and the connection to string theory 

The Af = 4: super Yang-Mills theory introduced here is strongly related to string theory, which becomes 
visible in the limit of asymptotically many colors, — >■ oo, while the effective coupling A = g^^N is 
kept fix. This is the so-called 't Hooft limit [TU]. The motivation to consider the large iV limit is to find 
a parameter which allows for perturbative calculations in strongly coupled gauge theories, namely 1/N. 
To illustrate this in a simplified way for the theory given by (j2.ip . we note that one can schematically 
write the interaction terms of this Lagrangian as 

L ^ Tr + 5™c'^'=$.$j$fc + , (2.4) 

where $i are any of the bosonic fields Xi or A (and the fermions are related to them by supersymmetry). 
Note that the three-point vertices are proportional to (7ym while the four-point vertices are proportional 
to ^YM- After the introduction of = gYM^i the Lagrangian acquires the form 



^Tr 

5?M 



(2.5) 



Remember that all fields $i of A/" = 4 SYM theory transform in the adjoint representation of the gauge 
group SU(iV). So can be written in a matrix notation, where the elements of the matrix are denoted 
by ($1)^5, with a,h = 1, 2, . . . , transforming in the fundamental and antifundamental representation, 
respectively. In a Feynman diagram a propagator for some particle then corresponds to a double 
line, with one line corresponding to the upper and one to the lower index a, b of the gauge group, see 
figure O 

In this double line notation we can now order diagrams in an expansion parametrized by 1 /N to see 
that the contributions to gauge invariant processes may be ordered according to the Euler characteristic of 
the Feynman diagram, i. e. they are ordered according to topology of the diagram. As an example consider 
the diagrams in fig. 12.11 For the amplitude corresponding to some Feynman diagram, a propagator 
introduces a factor of g^^ = X/N while the above Lagrangian (|2.5p shows that vertices pick up a factor 
of l/g^M = N/X. From the double line notation it is clear that each closed line therein represents a 
loop which introduces a factor N. Moreover, think of the diagrams as describing polyhedrons which are 
characterized by vertices, edges (propagators), and faces which are the regions separated by the edges. 
We observe that the factors of N for a diagram with V vertices, E edges and F faces (i. e. loops of lines 
in the double line diagrams) appear in powers of 

^V-E+F^E-V ^ ^XXE-V ^ ^2-2gxE-V_ (3.6) 

The number F — i? + F = x = 2 — 25is the Euler character of the polyhedron described by a Feynman 
diagram. The genus of the corresponding Riemann surface is given by g. From this dependence on N we 
see that diagrams with smallest i.e. planar diagrams with 5 = 0, contribute with highest order, while 
diagrams with topologies of higher genus g are suppressed by factors of N'^^ relative to the planar ones. 
In this way any process in the field theory can be decomposed into diagrams ordered by their genus g 
in the double line notation. The amplitude of a given process may then be obtained by a sum of the 
contributions from all relevant Feynman diagrams, 

oo 

X=^7V2-29;^(A). (2.7) 
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Figure 2.1: Feynman diagrams (left) can be translated to double line diagrams (middle), which in turn 
can be interpreted as Riemann surfaces of well defined topology (shaded) . These surfaces (deformed to 
the shape on the right) can be interpreted as stringy Feynman diagrams. 

This type of expansion is exactly the same as the one obtained by performing an expansion of diagrams 
describing the interaction of closed oriented strings, the type II string theories, upon recognizing the 
parameter 1/iV = 9^^/^ ^ being proportional to ^s, the string coupling constant [8 . The hope is, that 
the topologies of Feynman diagrams reflect the contribution of the string theory diagrams with the same 
topology. From the standard examples shown in figure ETTI we see that the genus g represents the number 
of loops in the associated string theory diagram. 

The large N limit corresponds to weakly coupled string theory, as gs oc X/N. In this limit we only 
have to consider the leading diagrams with genus 3 = 0. These are the gauge theory processes described 
by planar diagrams, corresponding to tree level diagrams in string theory. 

These arguments are of heuristic nature. For instance, there are effects like instantons in a gauge 
theory, which can not be treated in a 1 /N expansion. Such effects therefore should match the according 
non-perturbative effects in string theory. 

2.1.2 Type IIB supergravity 

The preceding section introduced the quantum field theory, which represents one of the two theories 
connected by the AdS/CFT correspondence. The second theory is type IIB supergravity. Supergravity 
theories are supersymmetric gauge field theories containing a spin 2 field identified with the graviton, 
the quantum field of gravitation. Supergravity thereby is an attempt to combine supersymmetric field 
theory with general relativity. 

Even though Supergravity is an interesting field to study on its own right, it can be embedded in 
a larger and more general framework. In fact supergravity is a certain limit of string theory. A brief 
comment on this perception will follow below. As string theory revealed that a consistent description of 
the forces and matter of nature requires ten spacetime dimensions, we will be interested in a formulation of 
supergravity in ten-dimensional backgrounds. There are different supergravity theories in ten spacetime 
dimensions, which can be constructed from compactifications of a unique causal unitary 11-dimensional 
supergravity theory [5] . 

The one formulation we will make use of throughout this monograph, and which is the one most 
intimately connected with the AdS/CFT correspondence is the so-called type IIB supergravity in ten 
spacetime dimensions. This theory is a A/" = 2 supersymmetric theory with a field content given by 
the bosonic fields G (symmetric rank 2, the metric), C (scalar, axion), $ (scalar, dilaton), B (rank 2 
antisymmetric, Kalb-Ramond field), A2 (rank 2 antisymmetric), (self dual rank 4 antisymmetric). 
The fermions of the theory satisfy Majorana-Weyl conditions and are given by two t/i^ (/ = 1,2, spin 
3/2 gravitinos of same chirality) and two fields A''^ (J = 1, 2, spin /1/2 dilatinos of same chirality, which 
is opposite to the chirality of the gravitinos) . This theory is chiral in the sense that it is parity violating 
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The action of type IIB supergravity may be written down in terms of the field strengths 



Fi = dC, 

^3 = dA2, 



H3 = dS, 

F5 ^ F5 ~ A H3 + A F3 



(2.8) 
(2.9) 



F5 = dAi 



(2.10) 



and then reads 



'S'lIB 



+ 




(2.11) 




A4 A H3 A F3 + fermions, 



where k is the Newton constant and TZ is the Ricci scalar. Additionally, at the level of the equations of 
motion one has to impose the self-duality constraint 



Type IIB supergravity as a string theory limit 

Starting from the Polyakov action to describe string world sheets, tachyonic string modes were discovered 
in the derived spectrum. These tachyons indicate an instability of the theory. To arrive at a stable and 
causal theory, one should remove these tachyonic excitations from the spectrum. To do so, one may 
modify the Polyakov action by introducing supersymmetry, and truncate the spectrum of physical states 
in a consistent way, a procedure called GSO-Projection (after the inventors Gliozzi, Scherk and Olive). 
This projection exactly leaves a spacetime supersymmetric spectrum. 

This procedure not only removes the tachyonic ground state from the closed string spectrum but 
additionally demands a number of 10 spacetime dimensions to preserve causality. The lowest remaining 
modes after the GSO projection represent the ground state of the remaining theory. For instance, in the 
Neveu-Schwarz sector of the string excitations this ground state happens to be described by massless 
excitations of strings. The so-called level matching condition demands this state to be generated from a 
vacuum |0) by the action of two creation operators, a left- and a right-moving one, a^lj^alj^ |0). These 
excitations may be described by a tensor valued field M with components M^'''. This field in turn 
decomposes into a symmetric part with components G^'' (describing the degrees of freedom of the 
graviton), antisymmetric components B^'^ (the B-field) and the scalar $ (the dilaton) that determines 
the trace of AI. 

Computing the masses of string excitations generated by more than two creation operators acting on 
the vacuum, unveils that these excitations describe fields which represent particles of finite positive mass 
proportional to 1/a'. In a low energy theory compared to the energy scale of inverse string length, or 
equivalently on length scales that do not resolve the stringy nature of the fundamental theory one may 
approximate the strings by pointlike particles, effectively described as a' — ^ 0. In this limit, however, 
all massive modes gain infinite masses and will not effect the low energy dynamics. The low energy 
theory may therefore only contain particles described by the massless supersymmetry multiplet to which 
the fields B, G and $ belong. The action remaining for these fields exactly describes the supergravity 
action. We are interested in the sector of closed string excitations with same chirality for the left and 
right moving excitations, which is called type IIB string theory and leads to type IIB supergravity in the 
low energy limit. 

In the action (|2.1ip above, the fields G, B and $ can be found in the first line, they originate from 
the Neveu-Schwarz sector (NS-NS) of the string theory fields, while the second and third lines contain 
the Ramond sector (R-R) contributions. 

Extremal p-brane solutions and D-branes 

Solutions to the supergravity equations of motion with non-trivial charges of (p-f l)-forms Ap^i are called 
p-branes. These solutions exhibit Poincare-invariance in {p+ 1) dimensions, their name thus stems from 
the number p of spatial dimensions included in this symmetry group. In this sense these solutions are 
higher dimensional generalizations of membranes, which one would denote as 2-branes in this context. 



*F5 = F5. 



(2.12) 
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Note that the flux / of the field strength — dAp+i througli some surface E is conserved since 

d/ = dFp+2 — 0, as is an exact and thus closed form. Moreover, the electric coupling of the 

p-form to the p-brane with worldvolume Sp+i of spacetime dimension p + 1 can be described by the 
diffeomorphism invariant action 

5p=Tp /^p+i. (2.13) 



The proportionality constant Tp denotes the tension of the p-brane. It has the interpretation of the 
energy or mass per unit area of the worldvolume, 

^''^5s(27r4)P+i- ^^-^^^ 

In type IIB supergravity, there are 0-forms, 2-forms and 4-forms, allowing for the following p-brane 
solutions. The 0-forms allow for (— l)-brane solutions, so-called D(— l)-instantons. Then there are 1- 
branes, charged under and thus coupling to the according solutions of the B-field. The two dimensional 
1-branes are identified with the worldsheet of the fundamental strings of the underlying string theory. 
They are called Fl-strings. The 1-branes which couple to the Ai field are called Dl-strings, and the 
3-brane solutions according to the field are called D3-branes. 

The magnetic analogon to the electric couplings are given by the Hodge dual field strengths. The 
magnetic dual field strength to in a ten-dimensional background, is the (10 — n)-form which has 
a (9 — n)-form field as its potential. This in turn couples to a (8 — n)-brane. In this way the type IIB 
field strengths F\ and F-^ allow for magnetic couplings to D7-branes and D5-branes. 

The naming of branes as Dp-branes we just saw, arises from string theory. As we can see, the Dp- 
branes are coupling to fields in the Ramond sector. The letter D is short for the Dirichlet boundary 
conditions such a brane imposes on the dynamics of the endpoints of open strings. Namely, in string 
theory D-branes are identified with the surfaces on which open strings end IJl [T2]. The endpoints 
of these strings then have a well defined position in the direction perpendicular to the brane, namely 
the position of the brane. Such a specification of a certain value for a actually dynamical quantity is 
known as the imposition of Dirichlet boundary conditions. It is believed that the p-brane solutions in 
the supergravity limit of string theory may be identified with Dp-branes in full string theory. 

The fact that p-branes are {p + l)-dimensional Poincare invariant imposes restrictions on the metric. 
For instance, some d-dimensional spacetime which supports a p-brane will include a Poincare invariant 
subspace with symmetry group IR^'^^ x SO(l,p). Additionally, one can always find solutions which are 
maximally rotationally invariant in the {d — p — l)-dimensional space transverse to the brane. Thus, in 
particular a ten-dimensional spacetime supporting D3-branes has an isometry group of x S0(1, 3) x 
S0(6). 

Analog to Reissner-Nordstrom black holes in general relativity, the possible solutions of supergravity 
backgrounds can be parametrized by the mass M of the p-brane solution and its RR charge iV, which are 
functions of two parameters r± , which can be interpreted as horizons of the solution. The case r+ < r_ 
exhibits a naked singularity and is therefore regarded as unphysical. In the limit of r+ ~ r_ the brane 
is said to be an extremal p-brane, while it is a non-extremal black brane for r+ > r_, with an event 
horizon. For details refer to ref. [8 . We restrict our attention to the case of extremal p-branes. 

The most general form of an extremal p-brane metric can be written in terms of a function iJ as [3] 

ds2 ^ H{yy^ v^,dx^'dx'' + H{y)^ df. (2.15) 

Here, the coordinates of the vector y parametrize the space transverse to the brane, and 77 is the {p -I- 1)- 
dimensional Minkowski metric. 

Supported by the insight that D-branes are dynamical objects of the theory ITT], one can adopt the 
point of view that the above geometry is generated by a stack of € N branes placed in an initially fiat 
d-dimensional Minkowski spacetime at a positions iji, with i — \,2^ . . . , N . Asymptotically far away from 
the stack one can therefore expect the whole spacetime to become flat again. String theory calculations 
then restrict the function H{y) to 

a ('47r^(5-p)/2r (I:iH^ rv'(''-''-^^/2 

g(j/)=i+E '^^ ^ "^^i:^ — , (2.16) 



1=1 



\y-yi\ 
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where gs is the string couphng constant and a' parametrizes the string tension. 

Of special interest for this article are D3-branes and D7-branes. For the introduction of the AdS/CFT 
correspondence, it is useful to look at D3-branes first. D7-branes will become an important ingredient 
for generalizations of the correspondence. 



D3-branes and Anti-de Sitter space 

There are several aspects which make D3-branes especially interesting. First of all, D3-branes by defini- 
tion introduce four-dimensional Poincare symmetry, the resulting ten-dimensional geometry for p = 3 is 
regular. Moreover, the solution for the axion and dilaton fields (C with Fi = dC and $ in (|2.1ip ') can be 
shown to be constants. In addition, the field strength is self-dual. For our considerations the metric 
will be a central quantity. Especially the case of N coincident D3-branes located at a position y^s in a 
spacetime of dimension d ~ 10 will be important. From (j2.16l) we see that the function H{y) in this case 
is given by 

Hiy)^l + p^. (2.17) 

We introduce the quantity R simply as an abbreviation, 

i?4 ^ 4TrgsNa'^. (2.18) 

However, a few lines below we will see that this parameter has a crucial geometric interpretation. With 
R we write H{y) as 

H{y) = l+ (2.19) 

\y~VnA 

By a coordinate shift we may always denote the position i/o3 as the origin of the y coordinates and set 
it to zero. The distance from the brane will be denoted by r = |^|. The metric (|2.15p generated by a 
stack of D3-branes therefore may be written as 

_\ \ 

/ d4\ 2 / 2 



ds" + dx^^dx^' + (dr' + r'dni) . (2.20) 

Far away from the stack of branes, at large r R, where the influence of the branes on spacetime will 
not be sensible, the metric is asymptotically flat ten-dimensional Minkowski spacetime. However, in the 
limit of r — >■ the metric appears to be singular. This limit is therefore known as the near horizon limit. 
In fact spacetime is not singular in this limit but develops constant (negative) curvature. Because space 
is flat at large r, but has constant curvature at r — this limit is also referred to as the throat region. 
In the near horizon limit at small r ^ R the metric asymptotically becomes 

^2 n2 

ds^ = — da;.da;'' + — dr^ + R^d^l (2.21) 
R'' 

This is the product space AdSs x , where the first two terms describe what is known as five-dimensional 
Anti-de Sitter space, or AdSs for short. The parameter R is called the radius of AdSspace. The last 
term represents the familiar five-dimensional sphere, of radius R as well. The geometry of Anti-de Sitter 
space is crucial for the gauge/gravity duality. To discuss some properties we introduce the coordinate 
z — R^ /r and write the metric as 



n2 

ds^ = _ (dx^da;'' + dz^) + R^d^l. (2.22) 

The metric (|2.22p can be derived as the induced metric of a five-dimensional hypersurface which is 
embedded into a six-dimensional spacetime with metric 

4 

dsl = -dXl + J2 dXf ± dXl (2.23) 
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where the Xi parametrize the six-dimensional space and the choice of the ambiguous sign depends on 
whether we aim for a metric on AdSs with Euchdean or Minkowski signature. Originally, the AdS/CFT 
correspondence was conjectured for Euclidean signature. The hypersurface which defines AdSs obeys 

Xf.X^' = (2.24) 

We now parametrize this surface with so-called Poincare coordinates z > and e R with fi — 0,1, 2, 3, 
such that Xij,x^ = ±(a;'')^ -I- (x^)^ -I- (x^)^ -I- (x"^)^ with the sign corresponding to the one in (|2.23l) and 



^0 


+ 


zU 






2z 


; 


Xi 




z^- 


dj l_Ldj 




2z 






= R-, 

z 




1 = 1,2,3, 


X5 


^0 

= R-. 

z 







(2.25) 



The hypersurface parametrized by z and the fulfills (|2.24l) and therefore represents AdSs. It has the 
induced metric 

7-)2 

ds2 = (dxndx^' + dz^) (2.26) 
z-' ' 

which appears as the first factor of the product spacetime (12.221) . Note that for the Minkowski signature 
background the restriction z > leaves only one of the two separate hyperboloids described by p. 241) . 
The other half is parametrized by 2: < and is a clone of the part we use. The spacetime coordinates 
parametrized by x suggest to be related to four dimensional Euclidean or Minkowski spacetime, depending 
on the choice of sign in (j2.23p . The coordinate z on the other hand is called the radial coordinate of AdS 
space. When we establish the AdS/CFT dictionary we will pay special attention to the behavior of fields 
near the so-called conformal boundary of AdS space. It is defined as the projective boundary which lies 
at z — > in the coordinates at hand. In the embedding space introduced above the boundary would be 
infinitely far away from the origin of the coordinate system. The metric (|2.26p . however, is diverging at 
the boundary, except we rescale it [T3]. A scale factor f{z) with a first order root of / at z = will 
exactly cancel the divergence after rescalings 

ds^ f{z)ds^. (2.27) 

As we are free to choose the function f{z) as long as we do not introduce new roots in f{z) or change 
the order of the root at z = 0, we can choose between a family of rescaling functions /, which are related 
by some arbitrary function w{z) as 

/(z)^/(z)e-(^). (2.28) 

This freedom therefore expresses the fact that the boundary of Anti-de Sitter space is only well defined up 
to conformal rescalings. Then the boundary at z = represents four-dimensional Euclidean or Minkowski 
spacetime, defined up to conformal rescalings. 

For later reference we point out the isometry group of AdSs ^ Sphere. The Lorentzian version with 
negative sign in (|2.23l) clearly displays an SO(2,4) rotational invariance in the AdSs subspace, while the 
isometry group of the five sphere is SO (6). 

We interpreted supergravity as a limit of string theory in the last paragraphs. Moreover, we will 
work in the near horizon limit from now on. Consequently, we will work in a background spacetime with 
the topology of AdSs x S^. The non- vanishing curvature of AdSs spacetime can be characterized by the 
Ricci scalar 

7^=-. (2.29) 

String theory, however, is not solved in curved backgrounds so far. To allow for a good approximation 
of type IIB string theory by working in the supergravity limit, one should therefore arrange spacetime 
curvature to be small. A large AdS radius leads to small curvature. Note that by (|2.18|) the relation 
of R to the string scale is = ^fd depends on two parameters of the theory. These are N and the 
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string coupling constant gs = e , which can be tuned by specifying a value of the arbitrary constant 
dilaton field $. We thus see that the supergravity approximation seems to be valid only for ^s-^ ^ 1- 
This guaranties £s, such that the radius of the string theory background is large compared to the 
string length ig. In this way gsN 1 ensures that the strings do not resolve the curved nature of 
the background, and type IIB string theory can be trusted as a good approximation to string theory on 
AdS5 X 

2.1.3 The Maldacena conjecture 

In a famous publication from the year 1997, Juan Maldacena pointed out that there exists a connection 
between certain quantum field theories and classical supergravity theories ^6j. In particular, the degrees 
of freedom found in type IIB supergravity on AdSs x S''^contain the large coupling limit of the A/" = 4 
SYM theory in four dimensions. 

As a generalization, consider full string theory instead of the supergravity limit, and relax the limit 
of large coupling on the quantum field theory side. Maldacena then conjectured the equivalence of two 
theories, formulated as the AdS/CFT correspondence. We summarize it as follows: 

Computations of observables, states, correlation functions and their dynamics yield the 
same result in the following two theories, which may therefore be regarded as physically 
equivalent. 

On the one side (AdS) there is 10-dimensional type IIB string theory on the spacetime 
AdSs X S^. The 5-form flux through the given by the integer N, and the equal radii R of 
AdSs and are related to the string coupling constant gs by i?^ = ^■Kg^Na''^ . 

On the other side (conformal field theory, CFT) of the correspondence there is a confor- 
mally symmetric four-dimensional N = A super- Yang-Mills theory with gauge group SU(A^) 
and Yang-Mills coupling ^ym, related to the string coupling by g^M = 27r(7g . 

This equivalence is conjectured to hold for any value of N and (Jym- 

It is a remarkable feature of this correspondence that it relates a theory containing gravity to a quantum 
field theory, which otherwise lacks any description of gravity. In the supergravity limit, a 10-dimensional 
classical theory of gravity matches a four-dimensional quantum theory. In fact a dictionary between 
operators of the quantum field theory and the supergravity fields can be established. We will comment 
on this below. However, the correspondence in its strong form, given above, is very general and thus allows 
hardly any applications. For instance, so far there is no formulation of string theory on curved spaces, 
such as AdSs x S^. Nevertheless, there are interesting non-trivial limits in which explicit computations 
can be performed. 

The 't Hooft limit is defined as considering a fixed value of the 't Hooft coupling A — g^^N while 
N oo. This yields a simplification of Feynman diagrams of the field theory. As we saw in section l^.l.ll 
in this limit only planar diagrams contribute to physical processes. Note that a fixed value of A in the 
large N limit implies weak coupling on the string theory side as the coupling constants are related by 
27r(7s = g^jj. So on the string theory side this results in the limit of a classical string theory (no string 
loops) on AdSs x S^. 

The Maldacena limit implements a further restriction. Starting from the 't Hooft limit, we let A —> oo. 
This of course prohibits perturbative computations on the field theory side, since here A is the effective 
coupling parameter. On the string theory side, though, this limit results in a' / R^ — > 0. So the curvature 
of the string theory background becomes small compared to the string length, which allows for consistent 
applications of the classical supergravity limit of string theory, which does not resolve the stringy nature 
of the fundamental building blocks of matter. 

Thus, working in the Maldacena limit not only allows to describe a quantum field theory in terms of a 
classical theory of gravity. It also allows to investigate the strongly coupled regime of the quantum field 
theory by performing calculations in the weakly coupled regime of the dual theory, where perturbative 
methods are applicable. 

The conjecture is not an ad hoc statement, but rather results from string theory arguments. Consider 
a stack of N coincident D3-branes which interact with open strings. In the low energy limit a' — >■ we 
have to consider infinitely short strings since £^ = a'. These strings may end on any of the N branes 
on the stack. As all the branes are coincident we can not distinguish between them, which implies an 
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U(A^) ^ U(l) X SU(A^) symmetry of the theory, where the U(l) factor gives the position of the brane and 
does not play a role here. It can be shown that the D3-branes' solutions exhibit N = A supcrsymmetry. 
Therefore in the low energy limit this theory describes precisely the conformal N — A SU(A^) gauge 
theory. Our special interest is the behavior of the strongly coupled regime of this theory, which is not 
accessible by perturbation theory. Instead of first taking the low energy limit and then the large coupling 
limit, we look at what happens if we proceed in reverse order. Starting from the stack of branes we are 
now interested in the strong coupling limit. From section r2.1.2l we know that the near horizon geometry 
in this case will have the topology of AdSs x S^with radius i?^ = AngsNa''^ , so we are forced to consider 
string theory on curved backgrounds. We also mentioned that the low energy limit of string theory is 
captured by supergravity. If we adopt the attitude that the physics of our system should be the same 
regardless of the order in which we impose the limits, then in the Maldacena limit we have to consider 
strongly coupled gauge theory and supergravity as two descriptions of the same physical setup. 

The quantum field theory may be interpreted as a description of the dynamics of open strings end- 
ing on the D3-branes. In the low energy limit a' — >■ the degrees of freedom (strings) are confined 
to the domain of the D3-branes. In the AdSs x S^geometry of the string theory background (I2.22p 
this domain is parametrized by the coordinates along the boundary of AdSs. So we can say that the 
AdS/CFT correspondence describes how a four-dimensional field theory defined on the boundary of 
five-dimensional AdS space encodes the information of a higher dimensional theory. In analogy to con- 
ventional holograms which encode three-dimensional information on a lower dimensional hyperspace 
(namely a two-dimensional surface), the AdS/CFT correspondence is said to realize the holographic 
principle. 



2.1.4 An AdS/CFT dictionary 

So far we recognized that the gauge/gravity duality allows for the reformulation of some problem defined 
in a gauge theory in terms of a gravity theory. In order to obtain quantitative answers, it is necessary to 
identify the corresponding quantities in both theories. The supergravity theory is formulated in terms of 
classical fields on a ten dimensional background, while the Af — A SYM theory describes the dynamics of 
operators acting on quantum states in four spacetime dimensions. The relations between the parameters 
of the theories were introduced with the correspondence in section [2. 1.31 For the coupling constants ^ym, 
gs and A, as well as the AdSradius R, the string tension a' and the number of colors N, they are 

^AirgsNa'^, 2^3^ = 5™, )^ = 9LN. (2.30) 

Observables, however, are expressed in terms of correlation functions of gauge invariant operators of 
the quantum field theory. It is possible to translate correlation functions of the field theory to expressions 
in terms of supergravity fields. A precise prescription of how to accomplish this was given in two seminal 
papers from 1998 by Edward Witten fT3^, and Gubser, Klebanov, Polyakov [TT. As a result it is possible 
to establish a complete dictionary, which translates quantities from on side of the correspondence to the 
other. 

Since the domain on which the field theory is defined can be identified with the boundary of AdSs 
space, one can imagine supergravity fields (j) in AdSs to interact with some conformally invariant operator 
O on the boundary. We denote the boundary value of the supergravity field by — hnigAdSs 4>- A 
coupling would look like 

5i„t= jd^xMx)0{x). (2.31) 

9AdS5 

In this sense the boundary value <j)^ of the supergravity field acts as the source of the operator O in 
the field theory. Such an interaction term appears in the generating functional for correlation functions, 
which we write schematically as 

/expf(j>oo\ . (2.32) 

Witten's proposal was to identify the generating functional for correlation functions of operators O with 
the partition function Zgugra of the supergravity theory, which is given by 

^sugra[0o] = exp (-S'sugraM) , (2.33) 

4>=<l>o 
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where S'sugra is the supergravity action. So the ansatz for the generating functional of correlation functions 
of operators of the field theory can be written as 



exp 



= exp (-S'sugra [(/>]) 



OAdSs 



'CFT 



(2.34) 



Correlation functions for O can then be obtained in the usual way by evaluating the functional derivative 
of the generating functional with respect to the source ipo of the operator. Explicit calculations will be 
performed in later sections. As a general example, some two point function would be obtained by solving 
the supergravity equations of motion, plugging these solutions into the action S'sugra, then expressing the 
result in terms of solution 0o on the boundary, and eventually evaluating 



(O(x)O(y)) 



S(t)oix) S<j)oiy) 



exp(-Ssugra[0]) 



(2.35) 



't>o=0 



The remaining question is which operators are dual to which fields. The fact that there exists such a 
dictionary relies heavily on the symmetries of the two related theories. The symmetry of a theory reflects 
the transformation behavior of the field content and by the Noether theorem accounts for the conserved 
quantities (charges). We can expect that two equivalent theories share the same amount of degrees of 
freedom, which must be reflected in their symmetries. 

To ensure a gauge invariant field theory action, including the source term (j2.3ip . we have to restrict 
our attention to operators O which are gauge invariant. The local SU(A^) gauge symmetry of the quantum 
field theory in fact has no counterpart on the supergravity side in the Maldacena limit. The parameter 

is translated into the number of D3-branes on the string theory side of the correspondence. The stack 
of D3-branes merely accounts for the emergence of the AdS^ x S^spacetime, see section [2. 1.2 1 Moreover, 
the following arguments strictly only apply to BPS states. 

Comparing the remaining symmetry groups of A/" = 4 SYM theory and type IIB supergravity we 
indeed observe a matching of symmetries. In section 12.1.11 we noted the symmetry group of the gauge 
theory to be PSU(2, 2 | 4). The bosonic subgroup of this is SU(2, 2) x SU(4)i^ ^ S0(2, 4) x S0(6). These 
are precisely the isometry groups of AdSs x S^, where SO(2,4) is the isometry group of the AdSs part, 
while the five-sphere is invariant under S0(6) transformations. The fermionic symmetries can be shown 
to coincide as well, leading to the overall symmetry group PSU(2, 2 | 4). 

In fact the isometrics of AdSs x S^act as the conformal group on the boundary [K]. Any gauge 
invariant field theory operator O does transform under some representation of the conformal group. 
Since the boundary theory is invariant under conformal transformations, the supergravity field in the 
source term (I2.31|) has to transform in the dual (conjugate) representation. Conformal invariance of 
the theory restricts the field further. For instance, in the coordinates where the boundary is located 
at M = the supergravity equations of motion for a scalar 4> have two linear independent solutions at 
asymptotically small u = e near the boundary. 



(j)o{x) e 



d-A 



(2.36) 



Here d denotes the number of dimensions of the boundary, which in our case is d = 4. Generically, the 
value of A for a scalar supergravity field depends on the mass of the field |I3j as 



A(A-d), 



(2.37) 



with A > 0. The second term of (j2.36p vanishes at the boundary e — > while the first term may diverge. 
The existence of a well defined boundary value 4'o{x) tells us that this function has scaling dimension 
d — A, i.e. <f)Q{x) 0o(x)/e'^~^ on rescalings. From the interaction term of the conformally invariant 
action (|2.3ip we thus see that the boundary value of 4){x,u) acts as the source to an operator 0{x) of 
scaling dimension A. 

In summary, to identify the supergravity field <j) dual to an operator O we have to spot all supergravity 
fields transforming in the dual representation to that of the operator under consideration. The conformal 
weight A of the operator determines the mass of the supergravity field by (|2.37p . The mass spectrum 
of A/" = 2 supergravity compactified on AdSs x S^has been computed [TB], and therefore the field can 
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be identified uniquely. Examples of computations of dual field-operator pairs can be found e. g. in refs. 

We also identified the boundary of AdS space with four-dimensional Minkowski spacetime. This 
spacetime was only defined up to conformal transformations, and we will identify it from now on with 
the domain of the conformally invariant A/" = 4 SYM theory. Notice that near the boundary all processes 
occurring in the field theory directions can be thought of as being scaled in such a way that all lengths 
of the AdS theory, even long distance or IR phenomena, are mapped to short scales, i. e. the UV limit on 
the conformal field theory side. To see this, consider the metric (j2.2ip in the near horizon limit. Then 
distances ds^prp in the field theory, which are measured along x appear with a warp factor relative to 
the distance ds^^jg in AdS space, 

'^^FT = ^d4ds- (2-38) 

At large r ^ 1 (IRacIs) close to the boundary, short scale phenomena of the CFT (UVcft) match the 
events in AdS space, while at small r <C 1 (UVAds)5 far from the boundary, long scale phenomena in the 
CFT (IRcft) match the AdS distances. The radial coordinate in this way sets the renormalization scale 
of the field theory which is holographically described by the supergravity theory. This phenomenon is 
called the UV/IR duality. 

In fact the behavior of correlation functions under renormalization group flows can be computed 
holographically. The UV divergences known from field theory translate into IR divergences on the gravity 
side. The procedure to incorporate scale dependence and renormalize n-point correlation functions is 
known as holographic renormalization. We will not review the procedure in detail here, but rather give 
an idea of the procedure and state some results. A nice overview which also addresses some subtleties 
can be found in ref. [T7] . 

The correlation functions (|2.35p in general suffer from IR divergences, i. e. divergent terms at large 
values of the radial coordinate. Analogous to quantum field theory renormalization, they can be cured 
by analyzing the behavior of the field solutions near the boundary and adding appropriate counterterms 
Set to the action S which do not alter the equations of motion but render the resulting correlators finite. 

To analyze the field behavior near the boundary it is convenient to work in coordinates u as in (I2.26P 
where the boundary is located at m — >■ 0. The solution for the second order equation of motion of any 
field can then be expanded in a series around u — 0. In general there are two independent solutions 
scaling as u™ and near the boundary. The general solution can be written as 

T{x, u) = (x) + f^^'Hx) + . . . + + /(2") In + . . .) (2.39) 

in a well defined manner where the coefficients /'•"■'(a;) carry the dependence on the other coordinates. 
The values of m and n are determined by the mass of the supergravity field and related to the conformal 
dimension of the dual operator, as in the example above. The coefficient /^"^ determines the boundary 
behavior of the two independent solutions for the equation of motion of Solving these equations order 
by order in u determines the relevant coefficients /'■'^^ for fc < 2n as functions of /^°\ which thereby can 
be used as the initial value of the first of the two linearly independent solutions. The second parameter 
needed to define the full solution of the second order equation of motion to T is the coefficient /^^"\ 
which in turn determines the remaining higher order coefficients. It then is possible to extract the 
divergent terms in the regularized action ^rcg, which is given by the on shell action with respect to the 
u dependence of the solution, evaluated at the cutoff e ^ 1, 

(2.40) 

The coefficients a^"'' now are functions of the coefficient /(°), and the v > solely depend on the scale 
dimension of the operator in the conformal field theory. Defining the counterterm action as 

Set — —divergent terms from 5rog (2-41) 

The renormalized action is given by 

5ren = lini (S-rcg + Set) ■ (2.42) 



Sic 



d^x a(") u- 
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Supergravity on AdSs x 



4 dim. TV = 4 CFT 



boundary of AdSs 

isometry of AdSs 

isonietry of S'^ 

weak coupling in gs 

variations in radial coordinate 

field boundary value 4>o 

field mass m 

IR normalizable mode 

IR non normalizable mode 

quantum corrections of 0{gs) 

stringy corrections of 0{a') 



renormalization group flow 



source of dual operator 
corrections in 



source for operator O 
conformal weight A 
(C) of dual operator 



field theory domain 
conformal symmetry 



strong coupling in A 



corrections in 1/ A 



R-symmetry 



Table 2.1: A few examples for entries of the AdS/CFT dictionary. Precise operator-field pairings can 
be found e. g. in refs. il [H [13] • 

Finding the renormalized action therefore involves a careful analysis of the equations of motion. An 
extremely useful result of holographic renormalization is the fact that the solutions to the equations of 
motion of a supergravity field J- can directly be related to the source and the vacuum expectation value 
of the dual operator in the field theory |17l I18j . In particular holographic renormalization unveils that the 
mode /(°^ is proportional to the source of the dual operator, while the mode J^^") is proportional to the 
vacuum expectation value of the same operator. In general, the mode that is proportional to the source 
scales in a non- normalizable way with u, while the mode proportional to the vacuum expectation value 
is normalizable. An example is given in (|2.36p for the case of a scalar field. The integers m and n are 
determined in terms of the supergravity field's mass, which in turn translates to the conformal dimension 
of the dual operator. In (j2.36l) the value of 4>a is proportional to the source and 0i is proportional to 
the vacuum expectation value of the dual operator. We will encounter an explicit example for the source 
and vacuum expectation value when we introduce the prominent pair of a D-brane embedding function 
and its field theory dual operator in section [^.2.21 

2.1.5 Tests and evidence 

Although a rigorous mathematical proof of the correspondence is not derived so far, there is increasing 
evidence for the Maldacena conjecture to hold. Soon after the discovery of the correspondence correlation 
functions of field theory operators where computed from gravity. A direct comparison to results obtained 
from field theory calculations is not straightforward since the results from gravity calculations are valid in 
the strongly coupled regime of the gauge theory while the gauge field theory computations are performed 
in the perturbatively accessible regime of weak coupling. 

Nevertheless, it was early realized that the answers obtained by gravity calculations gave the correct 
scaling behavior of n-point correlation functions, which is dictated by conformal invariance [13j . More- 
over, certain correlation functions satisfy non-renormalization theorems, which state that the results are 
independent of the coupling constant. Examples are the two- and three-point functions of 1/2-BPS oper- 
ators, which show a perfect matching of gauge and gravity results [191 [20]. Also the conformal anomaly 
of A/" = 4 SYM theory which is present in curved background spacetimes could be reproduced exactly 
from AdS/CFT, which even provided methods to obtain the anomaly in six-dimensional field theories 
for the first time pTj . 

Generalizations of the correspondence, which will be partly discussed below, relate different back- 
grounds to different gauge theories. Such modifications allowed for comparison of finite temperature 
A/" = 4 SYM with calculations obtained from the gauge gravity theory. These calculations show agree- 
ment of correlation functions in the hydrodynamic limit of low frequency/long distance [1U[13]. The 
recently most enthusiastically discussed hydrodynamic result from gauge/gravity calculations is the ob- 
servation of one universal value for the lower bound on the ratio rj/s of shear viscosity 77 over entropy 
density s for all known gauge theories with gravity duals [24] [25l [26] . In the large N limit the result in 
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SI units is 

! >- A, 

If it should turn out that QCD has a gravity dual and is in the same universality class as the known gen- 
eralizations of the correspondence, this result (including large A'' corrections) could be the first prediction 
from string theory within reach of experiment. Possibly the experiments at the RHIC and future collider 
experiments will give answers on the value of rj/s for QCD. So far the data seems to be in agreement 
with the above bound, which means that the QGP appears as the most perfect fluid ever observed. 



2.2 Generalizations and extensions 

Before we discuss the features of QCD and the sQGP which can be described by holographic duals, 
wc introduce those generalizations of the correspondence which are most relevant for this work. These 
generalizations are necessary to incorporate features which are missing in the jV = 4 SYM field theory 
described by supergravity on AdSs x S^. Most important for a description of the strongly coupled quark- 
gluon plasma are the inclusion of fundamental degrees of freedom and a way to describe systems at finite 
temperature. 

The introduction of finite temperature allows to model interesting qualitative features like the confine- 
mcnt/deconfinement phase transition at some critical temperature. In holographic models the transition 
occurs at different temperatures for the gauge fields and the fundamental degrees of freedom. The re- 
alization of finite temperature is achieved by a modification of the background geometry of the gravity 
theory. Moreover, wc will make some comments about the subtleties that are related to the computa- 
tion of Green functions at finite temperature. While correlation functions at zero temperature can be 
obtained in Euclidean spacetimes and a subsequent Wick rotation, this procedure generally cannot be 
applied at finite temperature. 

The inclusion of fundamental degrees of freedom is a generalization in the sense of additional fields 
we add to the gauge theory. All fields in J\f = 4 SYM transform in the adjoint representation of the 
gauge group and therefore rather account for the gauge degrees of freedom (gluons) than for the quarks. 
The added quark degrees of freedom will be represented by additional D-branes on the supergravity side. 
The incorporation of quarks into the theory furthermore allows to investigate the spectra of their bound 
states, which the subsequent section is devoted to. 



2.2.1 Finite temperature and AdS black holes 

At finite temperature T any quantum mechanical system can be found in one of the possible states of 
energy E with the probability distribution in equilibrium described by the density matrix 

where (3 = 1 /T is the inverse of the temperature and H is the Hamiltonian of the system. The statistical 
partition function of the ensemble of systems at temperature T can then be defined as 

Z,tat = Tr e-f^" = Y.{M e-^" \<t>n) , (2.45) 

n 

where the form a basis of the state space of the system. The statistical partition function defines 
the weight of each state that contributes to ensemble averages. Expectation values of some observable 
A in a thermal ensemble are calculated with respect to Zgtat by 

Tr (Ae-^") 

(A) = Tr(M) = — V -■ (2-46) 

■^stat 

In the quantum mechanical formalism of path integrals, transition amplitudes are given by 

<i>=<t>f 

{MtfMiU)) = {hitf)\ e-^(*^-*')« im)) = j^4> e*^f*l. (2.47) 
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The idea is to sum up all possible paths that evolve from the initial configuration (j>i to the final (pf. 
The complex phases give the weight for each possible configuration that contributes to the evolution. If 
we would not only consider one initial and one final state but sum over an ensemble of many possible 
states, we should therefore recover the sum of all weights, the partition sum. The action S[4i\ above is 
defined as 

*/ 

S[(t)]^ [dt [d'^-'^xL{t,x,(t)). (2.48) 



In quantum mechanics the standard method to obtain expectation values is the evaluation of functional 
derivatives of generating functionals, which are commonly also referred to as partition functions. They 
are defined for some functional >S'£;[(/)] by 

Zge„ = J Dc/^e-^-I"^!, (2.49) 

where we need to specify which functions (j) we have to integrate over and what the functional Se[(P] is. 
The imaginary time formalism gives a prescription which exactly reproduces the thermal equilibrium 
probability weights with the Boltzmann factor given in (j2.44l) . The prescription is to analytically con- 
tinue the time coordinate into the complex plane, such that t in (j2.48l) integrates over complex times. 
Additionally, we introduce a new time coordinate r = it as a Wick rotation of t. If we now restrict 
the system to such field configurations (f> that are periodic (in fact fermionic fields would have to satisfy 
anti-periodicity) along the imaginary axis in complex time t with tf — U — —il3 and /? G R between the 
initial ti and the final tf, we can reproduce the Boltzmann weights by setting 

= /dr /d'*-ia;L^(r,a;,0). (2.50) 







The index E refers to the fact that we use the Euclidean, i. e. the Wick rotated version, of the action. 
Then the integration from to /3 translates into integration over complex times ti to tf and therefore 
introduces the factor tf — ti = —if3 in (j2.47p . As we restrict to periodic states on the integration intervall, 
the final states {(j)f\ match the initial states \(j)i). Thus the path integral resembles a trace [27] . 

^gc„= / D(/. e-^-I"^! = Yl {Me~^"\<l>p)^Z,te.t- (2.51) 



Adding source terms to the action, which are set to zero after functional derivation, yields the Boltzmann 
factors as weights from the Euclidean generating functional. The imaginary time formalism in this way 
trades time for temperature and imposes boundary conditions. The result no longer depends on a real 
valued time interval but only on the purely imaginary time interval tf — ti ~ which we interpret as 
the temperature T by identifying T = . Abandoning time dependence is in accordance with the fact 
that we investigate a system in equilibrium, where expectation values do not change with time. Another 
consequence of the periodic boundary conditions is that any solution admits a discrete spectrum in its 
Fourier transformation. A propagator G(t) can be decomposed according to 

G{T)^-Y.e-'^--^GM (2.52) 

n 

The frequencies are called Matsubara frequencies. 

The link between field theory at finite temperature and gravity is known to be related to black 
hole physics, which has many parallels to thermodynamics. From the moment of the discovery of the 
AdS / CFT correspondence it was expected that the finite temperature description of the field theory must 
be given by gravity in the AdS black hole background [B]. We will now introduce the AdSs black hole 
background and give arguments for the relation between the horizon radius and the temperature of the 
dual field theory along the lines of arguments in refs. [251 [53] ■ The generalization of the AdSs x S^metric 
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(|2.2ip to the black hole solution with a horizon at r = Tq is given by, 



(2.53) 



The signature of the metric again depends on our choice of working in either Lorentzian or Euclidean 
AdS space. To draw the connection to finite temperature imaginary time formalism we again work in 
the Wick rotated coordinate r = it, where the metric has Euclidean signature. Euclidean signature 
however is only given outside the horizon, inside we would introduce negative signs from f(r) (in the 
Lorentzian case the signs of t and r would change). On the other hand it is known that the spacetime 
can be continued beyond the horizon and therefore the spacetime can be regularized at Tq. 

The idea is to show that periodicity in the Euclidean time r that leads to thermal probability distri- 
butions in field theory corresponds to a regularization of the Euclidean spacetime at the horizon. The 
period /3 — 1 /T is then identified with the inverse temperature as in field theory. We concentrate on the 
T and r coordinates in the above metric and observe how they behave near the horizon at r « ro, 

ds2 = ^{r- To) dr^ + ^{r~ ro)-'dr^. (2.54) 

We know that the Euclidean spacetime is well defined only on and outside the horizon and therefore 
introduce a new coordinate = r — Tq , which casts the metric into the form 

ds' = ^(^dg' + g'^dr'y (2.55) 

The factor in front is merely a constant. The metric in parentheses is the metric of a plane in polar 
coordinates, ds^ = dg^ + g^dO^. The angular variable in our case is 6* = TlrojB? . The space in polar 
coordinates, however, is regular only for an angular variable that is periodic with period 27r, otherwise a 
conical singularity is located at = 0, which is the horizon of the AdS black hole. Periodicity of Q with 
period 27r then translates into periodicity of r with period /3 by 

||r + 2.^||(r + /3). (2.56) 

Where we used the symbol ~ to denote the equivalence relation of identified points. Regarding identified 
points as equal and using /3 — \/T we obtain the relation between the temperature of the field theory 
and the horizon radius on the gravity side, 

To = TttR^. (2.57) 

This result exactly reproduces the expression for the Hawking temperature of the AdSSchwarzschild 
black hole described by the metric (|2.53p . 

The AdS black hole background allows to investigate strongly coupled gauge theories at finite tem- 
perature by performing calculations in the gravity theory. Moreover, finite temperature defines an energy 
scale (or length scale ro on the gravity side) in the theory and therefore breaks global scale invariance, 
which will have consequences for the property of (de)confinement of the field theory, which we discuss in 
section [^751 Despite the detour to a background of Euclidean signature with time t, the field theory we 
want to use as a model for real world QCD is defined on a background of Minkowski signature with time 
t, which we will use predominantly from now on. Some subtleties and consequences regarding Minkowski 
signature AdS / CFT are discussed in the following. 



Thermal real time Green functions 

The AdS/CFT correspondence was originally formulated and successfully applied in backgrounds of 
Euclidean signature, i. e. in the imaginary time formalism. Results in real time can be derived by 
subsequent Wick rotation. However, many common situations require the formulation of a problem and 
its solution in real time. One mathematical reason for the need of a real time formulation arises from 
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Figure 2.2: Integration contour in the complex time plane for the imaginary time (red) and the Schwinger- 
Keldysh (green) formalism for finite temperature field theory. 

simplifications which are often introduced by deriving solutions only for certain limits of the parameter 
space. In many cases for instance, solutions are obtained in the hydrodynamic limit of low frequency /long 
distance physics. In this case only the low Matsubara frequencies are known and therefore analytic 
continuation of results to real time is somewhere between difficult and impossible. Physical arguments 
against the imaginary time formalism arise whenever deviations or even far from equilibrium scenarios 
are considered. We discussed that the imaginary time formalism mimics thermal equilibrium probability 
distributions. For systems out of equilibrium the restriction of the path integral to periodic paths is not 
justified. Moreover the solutions considered in the imaginary time formalism are periodic, it is doubtful 
whether such solutions can model long time evolutions. Summing up, it is desirable to have a real time 
prescription for computations of correlation functions in Minkowski space. 

In field theory such a prescription is given by the Schwinger-Keldysh formalism |27| . The difference 
to the imaginary time formalism is a modification of the time integration in the complex time plane. 
The time coordinate is still defined on the complex plane. However, instead integrating from ti to ti — ij3 
along the negative imaginary axis, this time a detour along the real axis is taken. The path first proceeds 
along the real axis from ti to /3, which marks the end of the physical real valued time interval of interest. 
Then, the integration contour enters the negative half plane arbitrarily far and leads back below the real 
axis to Ret — Ret^ to continue parallel to the imaginary axis and finally end in — ti — i(3. Figure [52] 
shows the integration paths of the imaginary time and Schwinger-Keldysh formalisms. We are interested 
in the analog of the latter prescription in the gauge/gravity duality. 

A recipe for the derivation of holographic Minkowski space Green functions in real time was derived 
by Son and Starinets together with Herzog and Policastro [321 1221 ISO]- The resulting prescription is 
concise and amounts in only small changes from the prescription given by (I2.35|) . The difference is given 
by the way we compute the on shell action S'sugra- The action is commonly obtained by writing the 
solution of a field (j) such that the "bulk contributions" in radial direction r factorize from the "boundary 
contributions" along the field theory directions x on the boundary rb. Usually we work in momentum 
space with momentum k instead of position space with coordinate x and write 



In the coordinates introduced so far the boundary was located aX ri, — oo. The on shell action can then 
be written as _ 



Here, we carried out the integration over radial coordinates from the black hole horizon to the boundary. 
The two 4>Q arise from the kinetic term and the function T{r, k) collects the remaining factors of /(r, k) 
and drf{r, k) and possibly other factors appearing in the action under consideration. A detailed and 
explicit calculation can be found in section [3] where we apply the recipe, and in refs. [231121]. The two 
point Green function of the operator dual to (p would then be given by the second functional derivative 
of the action. 



,/,(r, k) = /(r, k) (fi^'^^ik), with lim /(r, k) = 1. 



(2.58) 




(2.59) 



G = - J-(r, k) 



F{r,-k) 



(2.60) 
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So far we followed the method we already introduced for zero temperature. The only difference is that now 
To ^ 0. The evaluation of this expression is mathematically possible, but gives physically wrong answers. 
For example, the resulting Green functions would be real functions, opposed to physical solutions, which 
are in general complex valued. This behaviour is due to the boundary conditions that have to be 
imposed on the fields. In the Schwinger-Keldysh formalism they arise from the periodicity of the fields 
and from the orientation of the integration contour. This introduces a contour ordering prescription, 
which translates to time ordering in physical processes. Loosely speaking, a causal propagator in the 
AdSblack hole background describes propagation of a field configuration that has to obey the infalling 
wave boundary condition at the black hole horizon. This boundary condition imposes the physically 
given fact that at the horizon, positive energy modes can only travel inwards, while negative energy 
modes only travel outwards. It can be shown that upon imposing these boundary conditions, the time 
ordered retarded part of the propagator in momentum space is determined by the boundary behaviour of 
the fields alone |30j . The prescription of Son and Starincts for fields obeying the infalling wave boundary 
condition then reads 

G''^-2F{r,k) . (2.61) 

The contributions from the horizon are neglected. This method was used to great extent in subsequent 
publications. At zero temperature it agrees with the analytic continuation of Euclidean results |23) . 



2.2.2 Fundamental matter — adding flavor 

We so far recognized that closed string excitations in the vicinity of the stack of N D3-branes gives rise 
to fields transforming in the adjoint representation of SU(iV), which was identified with the color group. 
We could criticize a lack of fundamental degrees of freedom, such as quarks in QCD. Karch and Katz 
introduced a way to add fundamental fields to the theory . The cure can be obtained from modes of 
open string excitations with one end of the string on the stack of N D3-branes and the other end on a 
different stack of Nf Dp-branes. In this work we will restrict to stacks of coinciding branes. We cannot 
distinguish between the branes of a stack on which a string ends. We thus encounter a U(iV) and a 
\]{Nf) symmetry which reflects the invariance of the theory under the exchanges of the branes [5^ I12| . 
The modes of these strings transform in the fundamental representation of SU(A^) C U(iV) and SU(iV/), 
respectively. We interpret these as the color and flavor groups. So the fundamental fields (quarks) of 
our gauge theory correspond to strings that have one end on the stack of N D3-branes, and the other 
on an additional stack of Nf Dp-branes which may be separated from the color branes. 

Throughout this work we will consider the so-called probe limit in which Nf <^ N . This ensures 
that the backreaction of the additional branes on the near horizon geometry of the D3-branes can 
consistently be neglected. In this way we do not have to worry about how the new Dp-branes might 
alter the background geometry but stick to AdSs x S^. The Maldacena limit of infinitely many colors 

— > oo is then also called the probe limit, since we add some neglectable amount of Nf probe branes. 

The string modes stretching from the D3-branes to the probe Dp-branes also transform under the 
fundamental representation of the probe branes' gauge group SU(A^/). However in the Maldacena limit 
with N oo the 't Hooft coupling Xf — 2TrgsNf of the stack of probe branes can be neglected with 
respect to the color gauge group coupling A = 2TTgsN. The probe brane gauge group in this way decouples 
from the color gauge group. We will identify the probe gauge group as the fiavor group and interpret 
strings stretching from the stack of Nf Dp-branes to the stack of N D3-branes as fundamental matter 
which comes in A^ varieties of color and Nf flavors. The additional Dp-branes are therefore also called 
flavor branes. For finite Nf, the large A" limit then is the equivalent to the quenched limit of lattice 
QCD, which allows to neglect fermion loops in all amplitudes relative to effects of the glue. 

The global U(A^/) flavor symmetry of the field theory translates into a gauge symmetry on the super- 
gravity side. The conserved currents of the field theory are dual to the gauge fields on the supergravity 
side. We will elaborate on this issue further when we introduce finite particle density. For now we 
only stress that the introduction of Nf flavor brane accounts for a gauge field on these branes which 
acquires values in a V{Nf) Lie algebra. We denote the field strength tensor of this gauge field by F. 
The components of this tensor are labeled by i^^^, where the fi and v denote spacetime indices while 
a = 1,2,..., Nj is an index in the vectorspace of the U(A^/) generators. 

The remaining issues then are, what dimensions the flavor branes should have and how they have 
to be positioned with respect to the D3-branes. Generically, D-branes couple to the field strengths of 
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type IIB supergravity, cf. (|2.64l) . Karch and Randall showed that there are stable probe brane solutions 
which span topologically trivial cycles and are determined by the DBI action alone [33] . There are several 
such solutions which then give rise to fundamental degrees of freedom in the dual field theory [3T] . 



Dirac-Born-Infeld action 

The dynamics of Dp-branes is crucial for the calculations performed in the following. Here, we introduce 
an action which allows to derive the equations of motion for D-branes. Later we will deal with stacks of 
Dp-branes, for now we consider the simpler case of a single brane. 

The interpretation of a single D-brane as the surface on which the endpoints of strings lie implies 
Dirichlet boundary conditions for the positions of these points. It is the Polyakov action that describes 
the dynamics of the strings. In the presence of background fields, a generalization of this action is 
given by a non- linear sigma model jl2| . The extremization of such an action respecting the Dirichlet 
boundary conditions is equivalent to the extremization of the Dirac-Born-Infeld action ^34j . This action 
captures the low energy dynamics of the string mode corresponding to the open string excitations of the 
Dp-brane. For a single Dp-brane with a worldvolume A4 parametrized by worldsheet coordinates ^* with 
i = 0, 1, . . . ,p, the DBI action is given by 

SuBi = -Tp J df e-*v/|det {P[g + B] + 2na' F)\ . (2.62) 

M 

Here g{£,), B{^) and F{^) — dA(^) are the background metric, the Kalb-Ramond i?-field and the gauge 
field strength tensor on the brane. The operator P[-] denotes the pullback on the brane worldvolume. 
The field $ is the dilaton. The brane tension Tp was given in (|2.14p . We will make extensive use of the 
DBI action. 

The DBI action is a low energy effective action that includes stringy corrections in a' up to arbitrary 
order. An expansion of the DBI action in powers of a' reproduces the Maxwell action in order and 
introduces higher powers of F as corrections. However, this action does not include any powers of the 
derivative of the field strength and therefore is strictly valid only for constant field strengths. For a 
DO-brane the DBI action resembles the worldline action of a pointlike particle. 

In cases where the Ramond-Ramond sector contributes non- vanishing n-forms C„ the full action for 
a Dp-brane is given by 

S = SoBi + '5'wzj (2.63) 

where 5'wz is the Wess-Zumino action 



Swz — Fp j P 



e^^^"'^. (2.64) 



However, all problems discussed in this work restrict to cases where there are no contributions from the 
Wess-Zumino action, 5wz = 0. For the case of the DO-brane, the WZ action resembles the coupling of a 
pointlike particle to an electromagnetic field. 



The D3/D7 setup 

Obviously, in the quark-gluon plasma the fundamental degrees of freedom move freely throughout the 
directions which are interpreted as Minkowski spacetime. Thus we will restrict our attention to probe 
branes which span at least all the Minkowski directions, i. e. we consider spacetime filling Dp-branes 
with p > 3. There are two heavily used models on the market, the D3/D7 setup and the Sakai-Sugimoto 
model. Throughout this work we will use the D3/D7 setup, which we introduce here. The Sakai-Sugimoto 
model will be addressed in a novercally short section afterwards. 

In the D3/D7 configuration, the background is generated by a stack of D3-branes in the way intro- 
duced above, which is then probed by TV/ flavor D7-branes. The D3-branes account for the background 
geometry and in the near horizon geometry give rise to the closed string excitations of type IIB super- 
gravity, accounting for the dual field theory A/" = 4 SYM. The additional degrees of freedom introduced 
by open string oscillations of string stretching between the D3-branes and the Nf probe D7-branes give 
rise to an A/" = 2 hypermultiplet in the fundamental representation of SU(Af) ^3lj|. The fermionic fields 
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in this multiplet, which we will denote by are interpreted as the analogon to quarks in QCD. The 
dynamics of quarks and mesons will therefore depend on the dynamics of D-branes in the holographic 
dual. An explicit realization of D7-branes embedded into AdSs x S^is given below. Applications of 
such flavor branes to describe quarks and mesons, their spectra and stability will be the subject of the 
following sections. 

Other D3/Dg setups 

Type IIB string theory exhibits D3, D5, D7 and D9 branes. We will investigate the D3/D7 model 
subsequently. Constructions with other types of branes can be interesting. However, with the quark- 
gluon plasma in mind, we do not investigate other probes than D7-brane for the following reasons. 

D9-brane can not be separated from the D3-branes since the former ones span the entire background 
spacetime and therefore cannot be separated from the D3-branes in order to generate massive quarks. 
D3-brane and D5-brane branes do not have this caveat. On the other hand they still have to span a 
certain cycles in the dimensions transverse to the D3-branes. Consider for example the embedding scheme 
of the D3/D5 setup. Here, we split the directions of the = ]R x transverse to the D3-branes, given 
by the radial coordinate of AdSspace and the directions along the S^, into a radial coordinate g and a 
cycle on the D5-brane, and (p, (p and L transverse to all branes. 



In this setup the ends of the string cannot move freely in the X3 direction of the Minkowski space- 
time. Configurations like the D3/D5 and D3/D3 setup therefore describe defect theories, in which the 
fundamental degrees of freedom are confined to lower dimensional hyperplanes. 

The Sakai-Sugimoto model 

The Sakai-Sugimoto model, introduced in refs. [23 El], describes the gravity dual to a Yang-Mills field 
theory in 3 -I- 1 dimensions where the gauge fields transform in the adjoint representation of the color 
group SU(Af), supplemented by Nf additional chiral fermions and Nf antichiral fermions which transform 
in the fundamental representation of the U(iV) and in the fundamental representation of a U(A^/) fiavor 
group. Supersymmetry is completely broken in this theory. 

The geometric realization of this setup is given by a D4/D8/D8 construction. A stack of N D4-branes 
in the near horizon limit gives rise to the background geometry of the (type HA) supergravity theory, 
analogous to the D3 setup. This time however one of the directions along the D4-branes has to be 
compactified in order to avoid a conical singularity in the resulting background. The matter fields are 
introduced by a number of Nf <^ N probe D8-branes and anti-D8-branes. These branes introduce the 
chiral symmetry groups \J{Nf)fi and V^Nf)^ which account for fermions of opposite chirality. A caveat 
of this model is that the bare quark masses of these fields are vanishing. 

Embedding D-branes 

As a concrete realization of the D3/D7 setup we now consider the embedding of a D7-brane into AdSs x 
S^and its thermal generalization, the AdSs x S^black hole background. We will perform the calculation 
in the black hole background and can obtain pure AdSs x S^solutions as the limit of vanishing horizon 
radius, Tq — > 0. Intuitive expectations would lead to embeddings which are infiuenced by the attractive 
gravitational force of the black hole as drawn in the cartoon of figure 12. 3[ which we will quantify now. 

The action of a probe D7-brane is given by the DBI action (|2.62p . For now we consider the case of 
vanishing field strengths 27Ta'F — B — 0. Therefore we are left with 



The constant prefactor Tp (cf. (|2.14p ') is not important for the following discussion. The action is 
determined by the induced metric G(^) = P[g{x)] on the worldsheet, where ^ are coordinates of the 
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(2.65) 
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Figure 2.3: Sketch of brane embeddings in the directions transverse to the D3-brane for different values 
of the temperature (relative to quark mass). Left: zero temperature, center: small temperature, right: 
high temperature, here the brane crosses the horizon. 



worldvolume of the D7-brane. The elements of the induced metric are given by 

GM-Q^Q^9a,. (2.66) 

Here g is the metric of the AdSs x S^black hole background with coordinates a;°(^) into which we embed 
the D7-brane. Two of these coordinates can be interpreted as functions which determine the position 
of the eight-dimensional worldvolume of the probe brane in the two directions transverse to the brane. 
These functions have to be determined in order to minimize the action (j2.65p . 

It is convenient for this purpose not to work in the coordinates of (|2.53p but to change to a new 
radial coordinate g given by 

(?2 = + y^r^ - r4. (2.67) 



The metric now is 



(2.68) 



In this way we can identify the transverse part to the D3-branes as nothing else than R and we write 
it as 

6 

dg^ + g^dnl = V dg"^ = dw^ + w^dnj + di^ + L^dcj)^ . (2.69) 
7^1 " ' ^ ' ' 

R*(ei 4) R=(e5,6) 

In these coordinates we parametrized the domain of the D3-branes by t and the three spatial coordinates 
X. The part of the spacetime transverse to it is parametrized by the six coordinates gi, with radial 
coordinate g — Cl2gfy^^- Equivalently, we wrote the transverse space as a product space of a four- 
dimensional in polar coordinates with radial coordinate w and a two-dimensional with radial 
coordinate L, such that g^ — w'^ + L'^ . 

An embedding of the eight-dimensional worldvolume of the D7-brane into AdSs x S^is then given by 
two functions which describe the positions in the two dimensions transverse to the brane. Stability of 
the D7-brane solution demands that the brane spans a trivial three-cycle in the transverse direction to 
the D3-branes [Slj. We thus embed the brane along the following directions. 

t Xi X2 Xj, w L (j) 

D3-brane 

D7-brane ^^^^^^^^^^^^^^ 



The worldvolume of the D7-brane is then parametrized by coordinates ^, 

^0 ^ ^1,2,3 ^ ^1,2,3^ ^4 ^ ^5,6,7 ^^Qj^g (2.70) 
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they determine the position of the D7-brane by the embedding functions and However, to 

ensure Poincare invariance the embedding functions cannot depend on ^o,...,3_ Moreover, the rotational 
S0(4) symmetry along the directions of the internal R"* of the worldvolume results in embedding functions 
which only depend on = w. The induced metric (j2.66p on the D7-brane then reads 



d4, = [-C dt^ + /dx^) + (d.^ + vj^dnl) , 

4 4 



(2.71) 



2 ■ 



This metric is AdSs x at asymptotically large w. Note that the embedding function (t){w) does not 
appear. This reflects the rotational symmetry of the setup in the space perpendicular to the brane. 
Further inserting the result into (|2.65p allows to derive the equation of motion for the embedding L{w), 



yi + (9^i)2 



v/1 + {d^Lf 



(w2 + £2)5 ' 

(2.72) 



^ (u;2 + L2)4^ 

From the asymptotic form of the equation of motion we see that the solution near the boundary at large 
w behaves as 

i = m^ + -^ + ... . (2.73) 

The embedding profile L asymptotically tends to a constant value TOi, = limu,_j.oo L, which we use as a 
free parameter of the setup. Together with the demand for smooth embeddings the boundary conditions 
for the solutions L{w) can be written as 

lim L{w) ^ TTir^, d^L{0) = (2.74) 
for embeddings that reach w = Q and 

lim L{w) = mj^, L{w) _L horizon (2-75) 

w^oo horizon 

for embeddings which enter the horizon. Any other boundary condition than orthogonality to the black 
hole horizon would lead to a transverse component of the gravitational force on the brane, which would 
deform the embedding until orthogonality is reached in the final equilibrium state. 

The differential equation (j2.72p generally has to be solved numerically 37|. At zero temperature, 
however, where Tq = 0, as well as in the limit of large g the equation of motion is solved analytically by 
a constant embedding function. Some brane profiles are shown in figure 12.41 The embeddings which do 
not touch the horizon have a regular worldvolume metric. They are called Minkowski embeddings. Note 
that these embeddings do not span the whole range of the coordinate g in AdSs x S^, since g^ — L'^ +w'^ 
and ming = min L{g) > Vq. From the induced metric (j2.7ip we see that these branes "end" at finite g 
before reaching the black hole horizon, since the S'^ wrapped by the D7-brane probe shrinks to zero size 
as in ref. [21] ■ Those embeddings that end on the horizon exhibit a black hole on their worldvolume and 
are therefore called black hole embeddings. 

Note that the supergravity scalar L is part of the radial coordinate r of AdSspace. At asymptotically 
large values of the radial coordinate, where L — m^,, the relation is 2r^ = + w'^. So at fixed w near 
the boundary L (x r. In the inverse radial coordinate u it scales like = u"^"^. According to the 
AdS / CFT dictionary, the dual operator of the d = 4 dimensional field theory therefore is of dimension 
A = d — 1 = 3. This operator is the bilinear ■0V'. From holographic renormalization we learned that the 
mode of the solution L scaling like u'^~^ is proportional to the source term of this operator. From the 
field theory Lagrangian L = mqipip + ... we see that this source is the mass of the "quark field" tp. 
An exact calculation relates the parameter m^, to the quark mass bjQ 

m . 

(2.76) 
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^Equation 1 12. 761 1 does not look like the formula for rriq given in the original paper, ref. |31l . where the concept was intro- 
duced. The relative factor of \/2 arises from the different coordinate systems used here and in ref. 1311 . The transformation 
between them introduces that factor in the embeddings and therefore also in the quark masses, cf. appendix 1X1 
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Figure 2.4: Solutions to (|2.72p yield Black hole (blue lines) and Minkowski embeddings (brown lines) 
of D7-branes in the AdSs x S^black hole background. The jump between black hole and Minkowski 
embedding at m^, = 1.3 induces a change of the worldvolume topology, reflected in a first order phase 
transition of the dual field theory, see refs. [571 138) . 

The embedding function in this way determines the mass of the quarks in the dual field theory [31] . 

Accordingly, the coefficient c in (|2.73p that scales like — is proportional to the vacuum expec- 
tation value Cc = {^ip), known as the chiral condensate. We will not discuss this quantity in detail in 
this work. 

2.3 Holographic quantum chromo dynamics 

In this section we want to point out some of the most important features of QCD and whether they can 
be described by gravity duals or not. We will see that the conditions at which the quark- giuon plasma 
exists, in particular finite temperature, allow for an at least qualitative description of many aspects via 
gravity duals. 

The gauge/gravity correspondence is a remarkable tool for the investigation of the strongly coupled 
regime of gauge theories. Depending on the choice of parameters, the features of the gauge field theory 
will be more or less close to what we expect from quantum chromodynamics. 

There is hope that the AdS/CFT correspondence may provide some pre- and postdictions even though 
the exact gravity dual to QCD is not known. During the last years it turned out that there exist some 
quantities, like the celebrated ratio rj/s = l/(47r) of shear viscosity to entropy density, that are universal 
in the sense that they do not depend on a particular supergravity background. Instead they are valid 
for all theories that have a gravity dual. If QCD is within this universality class, the results from other 
gauge theories than QCD may be applied to quantum chromodynamics as well. 

Field content and supersymmetry 

On the field theory side of the setup we stated that there is the gauge multiplet of A/" = 4 super Yang- 
Mills theory supplemented hy Nf Af ~ 2 hypermultiplets. The gauge degrees of freedom show up in a 
multiplet together with scalar and fermionic superpartner particle fields. In the Maldacena limit these 
multiplet exists in infinitely many colors, i. e. the rank of the gauge group is infinite. The fundamental 
fermionic flavor degrees of freedom on the other hand exist only in a flnite number Nf of flavors, mostly 
we will restrict to Nf — 2. The models studied in the following do not explicitly break supersymmetry 
further, while other setups may allow to break supersymmetry completely |39j . 

Compared to QCD we deal with infinitely times more degrees of freedom. However, as we saw the 
large N limit simplifies the theory drastically in the way that it allows us to neglect string loops. In fact 
corrections of results obtained from large N expansions with expansion parameter 1 /N tend to be small, 
as e. g. lattice calculations have shown. 

Moreover, at all finite values of temperature T ^ supersymmetry is broken spontaneously. Unlike 
most other spontaneously broken symmetries it is not restored at high temperatures |40) . For a system 
at equilibrium this can be seen from the fact that the fermionic degrees of freedom contained in a 
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supermultiplet obey antiperiodic boundary conditions along the imaginary axis in the complex time 
plane, while bosons are periodic. Therefore the Fourier decompositions and therewith the masses of these 
fields will differ. Fields with different masses however are not related to each other by supersymmetry 
transformations. This shows that supersymmetry is broken at finite temperature. 

To great extent we will be interested in the behavior of fundamental matter and the bound states of 
the fundamental fields. As we have done above when we identified the mass of the quarks with the mass 
of the fundamental fermions in the hypermultiplet, we will think of these fields as the cousins of the 
quarks in QCD — keeping in mind that gauge invariant operators receive contributions from the scalar 
superpartners. 

Conformal symmetry 

By definition the AdS/CFT correspondence relates gravity to a conformal, i.e. scale invariant, quantum 
theory. We introduced it as A/" = 4 SYM theory. This is in vast contrast to QCD where we know several 
scales which break conformal symmetry. 

The masses of the quarks which have to be regarded as fundamental parameters of QCD break 
scale invariance explicitly. The analogon to quark mass is realized geometrically in terms of the brane 
embedding profile. 

The dynamically determined momentum scale Aqcd at which the coupling constant diverges arises 
from quantization effects and therefore is a manifestation of scale anomalies. This scale has its dual in the 
background geometry of the supergravity theory. Temperature is geometrically realized by introducing 
a black hole into the spacetime of the gravity theory. The radius Tq of the black hole introduces the 
scale dual to the finite temperature in field theory. In the way a finite temperature can be interpreted 
as a lower bound on the momentum of particles, the black hole horizon radius introduces a cutoff. We 
interpreted the radial coordinate in AdSas the scale of a renormalization group flow. The horizon radius 
then works as a momentum cutoff. Geometrically, it introduces a scale and thereby explicitly breaks 
global conformal invariance. 

Bound states of quarks 

One of the great successes of the gauge/gravity duality is the possibility to derive spectra of bound states 
of fundamental matter from first principles. The string tension a' is a fundamental parameter of the 
theory and determines the quark mass, which naturally appears as a parameter in the spectra of bound 
states. Another parameter is the 't Hooft coupling A. Since we work in the low energy limit where we 
expect not to resolve the string scale we cannot explicitly assign any value to a', and in the Maldacena 
limit we cannot assign a finite value to A. We are thus unable to derive numerical values for meson 
masses. Nevertheless, the ratio of meson masses we obtain from holographic models can be compared to 
observations from experiment. Experimentally observed ratios of meson masses are reproduced with an 
accuracy of about 10% [4T . With respect to the various limits in which the according calculations are 
performed, this is an astonishing accuracy. 

Baryons The possibilities to model baryons is very limited to this day. This is partly due to the 
fact that we are restricted to the limit of infinitely many colors. Since baryons are colorless composite 
particles made out of Nc quarks, we would have to describe an object made of an infinite number of 
particles. 

There are however Skyrmion like solutions in the Sakai-Sugimoto model [36j. Recently, baryon like 
operators where considered in Chern-Simons- matter field theory derived in an AdS4/CFT3 model [32]. 

Mesons As mesons are composed operators containing one quark and one antiquark field, they trans- 
form in the adjoint representation of the flavor gauge group S\J{Nf), which can also be expressed as 
a bifundamental representation with one index in the fundamental representation Nf and the other 
in the antifundamenal Nf, as we did in section [2.1.11 for the color gauge group. This transformation 
property on the string theory side is given by a string that has both ends on the probe D7-branc. So the 
mesons of the fleld theory are dual to the excitations of D7-D7 strings. The endpoints of these strings 
determine position of the probe D7-brane. Consequently the excitations of the D7-D7 strings describe 
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fluctuations of the probe branes. The meson masses can then be obtained from the solutions to the 
Unearized equations of motion of these fluctuations of the probe branes around the embedding. 

As a very short sketch of the procedure and for later reference we outline the calculation of the 
spectrum of scalar mesons at zero temperature, first published in ref. [43]. Above, the embedding of the 
flavor D7-branes at zero temperature was shown to be described by constant functions L^w) — niL and 
constant (j}{w). We now allow for small deviations from this embedding by adding small fluctuations 
'fL,<i)(C) to the embedding functions, 

L^L + ^l{0. ^^4> + M0- (2-77) 

We want to consider small deviations from the brane profile and therefore may restrict our attention to 
the linearized equations of motion for the fluctuations ^l,4>- In the same way in which the embedding 
is determined by the equations of motion obtained from the DBI action, we can derive the linearized 
equations of motion for the fluctuations from the same action (|2.65p . Analogous to the derivation of 
the equation of motion (j2.72p for the embedding functions the equations for the fluctuations around L 
in the zero temperature case of Tq =0 are obtained by plugging in the ansatz (12.771) into the action. 
The resulting linearized equation of motion for the fluctuation (p was calculated in ref. |43j for zero 
temperature. To make the connection to the computation of the derivation of the embedding functions, 
we stick to the coordinates (I2.68P with (|2.69p in which the fluctuation equation for fluctuations (p^ are 
derived to be ^ 

= , d'd,^ +^d^ {w^d^^) + \ V^V,^. (2.78) 

Here i is summed over the Minkowski directions, Va are the covariant derivatives along the directions of 
the spanned by the probe D7-brane, and the radial coordinate of AdSis given by = w'^ + L^. Use 
the ansatz 

<fi^ip{w)e-'''^Y^{S^) (2.79) 

with y'(S'^) as the spherical harmonics along the three sphere, such that V°Va<^ = + 2)<f. Here 
Z = 0, 1, 2, ... is the angular momentum number on the S^. However, in this work we only consider the 
solutions with / = 0. The plane wave factor is responsible for d^di^p = —k'^cp, with momentum vour 
vector k which determines the meson mass M by = —k"^ . Therefore the above ansatz transforms 
(|2.78p into an ordinary differential equation for the radial part (p{w), 

= , 2^rr2A2 ^-^'^(w^) + ^d^ {w'd^,ip{w)) . (2.80) 
It can be solved in terms of hypergeometric functions. However, normalizable solutions only exist for 
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M„ = ^J{n + l)(n + 2) , n = 0,l,2,..., (2.81) 

where the quark mass enters through the embedding L by (|2.74p and (j2.76p , with constant L = . 
This is the mass spectrum of mesons at zero temperature and vanishing particle density. In fact this is 
the form of the spectrum for scalar, pseudo scalar and vector mesons [43j . We will compare later results 
at finite temperature and finite density to this formula. 

Various aspects of meson spectroscopy have been under investigation, among these are the discrete 
meson spectra of stable quark-anti quark mesons at zero temperature |43j , the decreasing stability and 
melting of these states at finite temperature and finite particle density [H [3l [44] , and the investigation 
of the spectra of heavy-light mesons [45l |46] . 

We will derive meson spectra for various purposes. On the one hand side we are interested in the 
dependence of the spectra under variation of temperature and particle density in order to understand 
the behavior of bound states of quarks. On the other hand we will observe the influence of external flelds 
on the mass spectra to derive the polarizability of the mesons, which in turn influence their diffusion 
behavior inside the quark-gluon plasma. 

For mesons in the Sakai-Sugimoto model we again refer to refs. [351 136) . 
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Confinement / deconfinement 

The probably most prominent feature of QCD is the running couphng constant, meaning the change in 
the value of the coupling constant of the strong interaction under variations of the energy scale of the 
interactions. Processes involving high momentum transfer are influenced less by the strong interactions 
that those which occur at low momentum. Mathematically, the value of the coupling constant even 
diverges at a momentum scale known as Aqcd • Such a running of the coupling constant with respect to 
the energy scale is obviously only possible in the absence of conformal invariance, which would forbid the 
existence of a characteristic scale. As a result, quarks at low energies, e. g. low temperature, are confined 
to bound states which appear as colorless entities to a far away located observer. At high energies/high 
temperatures, the quarks may escape from these states and travel through spacetime independently. 

Experiments show us that at energies above approximately 175 MeV quarks and gluons start to enter 
the deconfined regime. The exact value depends on various parameters. So far there is no analytic 
proof for these properties of quarks and gluons. To great extent this is due to the fact that the confine- 
ment/deconfinement transition occurs in the strongly coupled regime of the gauge theory. Traditional 
perturbative methods may not be applied here. 

In the framework of the gauge/gravity duality, however, one can hope to see effects of the confine- 
ment/deconfinement transition, since we can work in the strongly coupled regime of the gauge theory. 
The original correspondence contained adjoint matter fields, given by the gauge multiplet of A/" = 4 
SYM theory. At finite temperature, the gauge fields undergo a first order phase transition at a temper- 
ature Tgaugc- It coincides with the Hawking-Page temperature and can be interpreted as the confine- 
ment/deconfinement transition ITT] . 

The fundamental matter existing in probe brane setups also exhibits a phase transition at finite 
temperature, which though occurs at a different temperature Tfund than the transition of the gauge 
fields. Various different models exist where the fundamental degrees of freedom indeed undergo a phase 
transition from stable bound states to dissociating ones [371 UHl HI Hll EH] • We will come back to this 
transition when we discuss the QCD phase diagram in section [5l 

It is interesting to note a difference between holographic and lattice models. The deconfinement 
temperature for fundamental matter associated with the destabilization of mesons derived from holo- 
graphic models is proportional to the mass of the constituent quarks of the meson, Tfund oc m^. Lattice 
results in the quenched approximation, in contrast, suggest a scaling of the transition temperature of 
meson destabilization with the transition temperature for the gauge fields Tfund oc Tgauge [491 150] . This 
in principle allows for interesting comparison of lattice and holographic models with experimental data. 

Chiral symmetry 

Soon after the introduction of fundamental matter to the gauge/gravity duality it was shown that various 
probe brane setups are capable of realizing chiral symmetry breaking at finite temperature holographically 
[571 [551 [Sn [SH . The order parameter for the transition between the chiral symmetric phase and the 
phase of spontaneously broken chiral symmetry is given by the chiral condensate Cc for massless quarks, 
i.e. the vacuum expectation value of the bifundamental Cc = In some models this transition 

coincides with the confinement/deconfinement transition 1551. 

In the Sakai-Sugimoto model, chiral symmetry breaking is realized by the merging of the embeddings 
of the D8 and D8 in the low temperature phase. In this way the flavor groups U(iV/)i and \J{Nf)j^ 
originating from string excitations on the respective brane combine to a single vector subgroup U(A^/)v'. 

3 Thermal vector meson spectra at finite particle density 

In this section we address the first two questions raised in the introduction. We wondered whether bound 
states of quarks can be observed in holographic models of the thermal quark-gluon plasma, and how they 
are influenced by the medium. From experiment we know that temperature and particle density influence 
the interaction between particles, a prominent example is the transition of quarks and gluons from the 
confined to the deconfined phase at increasing interaction energy, or equivalently at high temperature. 

Temperature and quark density are the most important parameters of the model we make use of 
in this work. Together with the mass of the fundamental fields they define the axes of the phase 
diagram of fundamental matter in the holographic QGP. In the context of gauge/gravity duality, there 
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has been an intensive study of the phase diagram of A/" = 4 supersymmetric SU(A^) Yang-Mihs theory 
in the large N hmit, with fundamental degrees of freedom added by considering the AdS-Schwarzschild 
black hole background with D7-brane probes [SU |3H1 EZl ESI 113 HI]- Another approach was pursued 
by studying string worldsheet instantons [57] ■ Subsequently, particular interest has arisen in the more 
involved structure of the phase diagram when the baryon chemical potential is present, taking finite 
density effects into account ^58^. We contemplate the phase diagram and its parameters more detailed 
in section [5j In this section we concentrate on the mesonic bound states of quarks and their dependence 
on temperature and density. 

The aim of this section is the combination of both, finite temperature and finite density effects in 
the description of a thermal holographic plasma. In J\f — A SYM theory with finite baryon density, we 
relate our work to the phase diagram shown in figure 2 of ref . [SI] , reproduced below in figure 13.11 We 
restrict to setups with non-vanishing particle density. Here, fundamental matter is described solely by 
probe branes with the geometry of black hole embeddings [55] [SD] . The holographic realization of finite 
particle density is discussed below, where we introduce the setup. 

The concept of mesonic bound states in Af — 4 SYM theory in the Maldacena limit at vanishing 
temperature and particle density, and an outline of how to obtain their spectra was sketched in section [^751 
Here, we extend these calculations to incorporate in-medium effects of finite temperature and particle 
density, giving rise to non-vanishing baryon or isospin chemical potential. The motivation to do so stems 
from the possibility to conduct experiments at non- vanishing isospin density [61] . as well as the better 
accessibility by lattice methods of the finite isospin region in the phase diagram of QCD compared to 
finite baryon chemical potential. The work presented here restricts to the calculation of vector meson 
bound states with vanishing spatial momentum. At finite momentum the vector mesons couple to scalar 
mesons. Extensions of our work to finite momentum can be found in refs. [621 163] . 

The meson spectra will be represented in terms of spectral functions. These functions of an energy 
variable will exhibit resonance peaks of finite width, corresponding to decay rates, at energies corre- 
sponding to the meson masses. The necessary extensions of the setup to finite particle density and 
the concept of spectral functions are introduced in the subsequent sections. In order to determine the 
spectral function at finite temperature and finite baryon density, we make use of the methods developed 
in the context of AdS/CFT applied to hydrodynamics, c.f. for instance refs. [531[MlinS]. For vanishing 
chemical potential, a similar analysis of mesons has been performed in ref. [35] ■ There it was found 
that the mass spectrum is discrete for quarks with masses significantly above the energy scale set by the 
temperature. At lower quark mass, a quasiparticle structure is seen which displays the broadening decay 
width of the mesons. As the mass decreases or temperature rises, the mesons are rendered unstable, 
reflected in broad resonance peaks. These excitations dissipate their binding energy into the plasma. 
Note that for this case, there are also lattice gauge theory results [5S] . 

The achievements of the work presented in this section are the successful incorporation of either 
baryonic or isospin chemical potential at finite temperature. Before the results of this section where 
published as refs. [JJ |3] these aspects where investigated separately in the literature. As we will see, 
the simultaneous incorporation of both temperature and particle density leads to spectra which can be 
compared to previous publications consistently in the appropriate limits. In particular, we find that at low 
temperature to quark mass ratio, i. e. close to the Minkowski phase, where the characteristic energy scale 
of the system is given by the quark mass, the spectrum is asymptotically discrete and coincides with the 
zero-temperature supersymmetric meson mass formula found in ref. |43| and rephrased in equation (j2.8ip . 
However, away from this regime the dominant energy scale is either the finite temperature or the chemical 
potential. Here the observed spectra differ qualitatively from the above in some respects and resemble 
aspects of mesonic excitations in QCD. They also show interesting similarities to phenomenological 
models. We elaborate on the physical characteristics of our results in the summary of this section. 

In the case of an isospin chemical potential, previous work in the holographic context has appeared in 
refs. [Ml in?] . In this case, two coincident D7-brane probes are considered, which account for fundamental 
matter of opposite isospin charge. We find that spectral functions quantitatively deviate from the 
baryonic background case. A triplet splitting of quasi-particle resonances in the spectral function is 
observed, which depends on the magnitude of the chemical potential. 
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3.1 Spectral functions 

The spectral function 9l{uj, q) of an operator J describes the probability density in (w, q)-space to detect 
the quantity encoded in the eigenvalues of the operator J at given energy w and spatial momentum q. 
In our case, we want to describe quarkonium states and are interested in the mass/energy spectrum of 
the stable bound states and their lifetimes. In other words, we want to compute the spectral functions 
$H(w, q) of a quark-antiquark operator corresponding to vector mesons. This operator appears in the 
field theory as the flavor current J{x) = 'il}{x)jip{x) of fundamental fields ip{x) (and their superpartners) . 
For simplicity, let us restrict to the case of vanishing momentum, q = 0, where the remaining parameter 
is the energy w alone. Peaks in the spectral function at an energy tu indicate that there is a large 
probability to find a quark-antiquark state, which is denoted as a quasiparticle if the width of the peak 
is small compared to the height. The position tu of the peak gives the energy or mass of the quasiparticle 
while the width of the peak translates into the lifetime of this particle in position space. According to 
Fourier transformation, a broad peak, which is a large object in momentum space, corresponds to an 
event of short lifetime in position space, and vice versa a narrow peak in the spectral function is a signal 
for a particle with a long lifetime. 

We very briefly comment on how to derive the spectral function from two point functions, and how 
to extract the relevant information from them. See textbooks like ref. [68 for details. The formulation 
of spectral densities in terms of two point Green functions is convenient because we can compute the 
latter holographically. 

We think of J being the operator that describes the free mesonic quasiparticle as an excitation of 
one of the possible QGP many-particle states \n). There are infinitely many different of such states in 
the thermal ensemble that represents the QGP. The probability to occupy one of them is given by the 
density matrix p, described in section 12.2.11 These states form a basis of the Hamiltonian H of the 
ensemble, such that J2n l'^) ('^l — ^■ 

The probability of propagation from an initial spacetime point Xi, which we define as Xi = (0, 0), to 
some final point Xf — (t,x) is given by the time ordered Green function 

Gnit,x) = ~i{n\0{t) J{xf)f{0)\n) , (3.1) 

where the step function 0{t) accounts for time ordering. The index n shall remind us that this probability 
is not an ensemble average. It just gives the probability for the event if the QGP is in the state \n). 
By switching to the Schrodinger representation of the meson operators and denoting the momentum 
operator by k, 

J{t, x) = e-'^^''-^*) je*(kx-fft) ^ (3 2) 

and insertion of a full set of eigenstatcs \n') we arrive at 

G„(t,x) = -«^0(i)e*(^"-^"')*-^'^"' {n\j\ n') {n'\j^n) . (3.3) 

n' 

We denote the energy difference of the excited system to the QGP ground state as the energy w„„' of the 
mesonic excitation. This energy certainly depends on the state of the plasma with the excited mesonic 
state \n') , and the state \n) it was created from, w„„' = £"„ — We write (?^| J| jt-') (ti'] n) = 

|<^n'| J^^l and perform a Fourier transformation with respect to energy and momentum (a;,q), 

G„(c.,q)=5: Jii^i^- iH^i-')r^ (3-4) 

n' 

where the small e S R accounts for proper convergence. The delta function reflects that the momentum 
is conserved in the multiparticle system. To get the probability for the detection of a meson with energy 
uj and momentum q in the QGP, we have to perform the ensemble average. This eventually leads to the 
relation 

q) = ^ E J ^ ''rrl - ^) I I' ■ (^-^^ 

Zj ^nn' — Lu ^ Ic 

It is convenient to write this Green function as 

G(.,q)= /d.'^i^^. (3.6) 
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Here, we defined a weight function , k) for the propagation of the meson state, which assigns different 
probabilities to the propagation according to the Green function G'(w,q) = — w). The probability 

density is called the spectral density or spectral function. From p.Sp and p.6p we see that the spectral 
function is given by 

miLj',k) = ^J2H^'- ^nn')S (k - q) |(n'| J^l (l + e-^"') . (3.7) 

This notation reflects the physical interpretation we gave earlier. The probability to observe the multi- 
particle system in a state \n') with a mesonic excitation, created by acting with on the initial QGP 
state |n), obeys energy and momentum conservation, and depends on the temperature, given by /3. 

Most important for our purpose is the retarded Green function G^. We where not explicit about the 
retarded and advanced contributions to the Green function in the above discussion. Nevertheless, the 
Sokhatsky-Weierstrass theorem in complex analysis allows to derive the following relation between the 
spectral function and the retarded Green function that we will make use of, 

$n(w,k) -2ImG"(w,k). (3.8) 

The large probability density for the propagation of a quark-antiquark pair with the right energy content 
to form a bound state directly translates into an excess of the spectral function at that particular value 
of UJ. In the rest frame of the particle, which we are restricting our attention to, the energy can directly 
be translated into the mass M of the meson by w = M. The calculation of the Green functions of flavor 
currents J in this way yields information about the quasiparticle spectrum of a given theory — the meson 
spectrum. 

From the relation (j3.8p . we see that a convenient way to obtain the spectral function is to compute the 
retarded Green function of the mesonic operator. A way to achieve this was sketched in section [XTI We 
see that all information about the spectrum is contained in the correlation function C. The correlation 
function in turn is determined by the residues of its poles in the complex plane. From field theory we 
know that the poles of the correlation function in the complex w-plane can directly be translated to the 
energies of the states. We will consider spectral functions at vanishing spatial momentum q, determined 
by the energy w alone. 

The standard example in field theory is Klein-Gordon theory which amounts to the equation of motion 
for a scalar field 4> given by 

(□-m2)0 = O. (3.9) 

In terms of the formalism of Green functions the evolution of a delta-shaped initial perturbation of (j) is 
given by the inverse of the differential operator (□ — m^) . The modes of this solutions are then obtained 
from the Fourier transform (with □ in our example of a particle at rest), 

G{uj) (X . ^ „ . (3.10) 

The Green function exhibits poles at w = ±m, corresponding to modes with the energy uj of the stable 
particle at rest. These real valued poles are less frequently referred to as normal modes. The solution to 
more complicated systems than Klein-Gordon theory, where we have unstable excitations which dissipate 
energy, we encounter quasinormal modes £7 which are complex valued. Quasinormal modes (QNM) where 
introduced in the context of metric fluctuations in black hole background. The black hole geometry 
accounts for the attenuation since any amount of energy that crosses the event horizon irreversibly 
disappears from the system outside. This mechanism also works in the case we will investigate. The 
difference to the simple example of Klein- Gordon theory is that we consider correlators of gauge field 
fiuctuations in an AdSs x S^black hole background rather than a scalar field in fiat space. The fluctuations 
of the gauge field on the brane can transport energy into the black hole, but no energy can escape from the 
horizon. This introduces dissipation, which is then described by the imaginary part of the quasinormal 
modes. 

For demonstration, presume the solution (l){t) of the equation of motion for a mesonic excitation 
may be decomposed into quasinormal modes K{VL) for complex frequencies VL. Suppose we only excite 
one mode for a single complex £7', such that K{Vl) = k5{Vl — Q,'). This mode describes an attenuated 
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oscillation as long as the imaginary value of 17 is negative, 

=y"df7if(f7)e-*"* = A:e-*^'=("')*e^'"("')*. (3.11) 

For positive imaginary parts of the quasinormal modes, we encounter the unphysical case of infinite 
amplification of any fiuctuation of the field. Therefore, in a physical setup one may find singularities of 
the retarded Green functions G^(w,q) only in the lower half of the complex w-plane. Those with the 
lowest absolute value of the imaginary part are referred to as the hydrodynamic poles of the retarded 
real-time Green function since they determine the long time behavior. Consider the made up example 
in which a Green function exhibits a pole at fl = ujq — iT , 

G«(a.) (X (3.12) 

The pole emerges as a peak in the spectral density of real valued energies to, 

mioj) = -2ImG«(c^) cx , (3.13) 

(a; - cjo)^ + r-^ 

located at loq with a width given by F. These peaks are interpreted as quasi-particles if their lifetime 
1/F is considerably long, i.e. if F ^ and thus the peaks in the spectral function are narrow. 

We compute the spectral function for real and complex values of the energy uj in this section. However, 
we focus on the meson spectra, i. e. the 9\{uj, 0) for w G H, and postpone the discussion of the physical 
consequences and a detailed analysis of the behavior and location of the quasi normal modes to section [5j 
In-depth analytical and numerical investigations of the behavior of quasinormal modes in gauge/gravity 
duality can be found in refs. [5^ 144) . 



3.2 Holographic setup 

For an explicit calculation of the flavor current Green functions, we have to find the solutions to the 
equations of motion for the supergravity field, which is holographically dual to the flavor currents in the 
field theory. This field, which we called <j> in the general example (|2.58p . is the gauge field A on the 
probe brane [33] . This means that we have to include the contributions of the gauge field strength tensor 
F = dA in the DBI action. 

Because of the non-linearity of the DBI action, the embedding functions and gauge fields couple. We 
no longer can expect to determine the embedding functions in terms of the quark mass alone. Instead, 
we first introduce a suitable coordinate system to describe the background geometry, then derive the 
equations of motion and solutions of the background fields, which are the embedding and the gauge field 
on the D7-branes. In the subsequent section we investigate the fluctuations of the gauge flelds on this 
background to eventually derive the meson spectra from them. To account for flnite temperature, we 
consider finite values rb in the AdSblack hole metric. 

We compute the functional dependence of the gauge field A{p, k) and the embedding functions nu- 
merically in the limit of vanishing spatial momentum q — > for the fluctuations. In this limit the 
momentum four vector simplifies, k — (oj, q) = (a;, 0) and there are no couplings between the vector and 
scalar mesons. 



3.2.1 Background geometry and supergravity action 

We work in the D3/D7 setup introduced above, i.e. we consider asymptotically AdS^ x space-time 
which arises as the near horizon limit of a stack of N coincident D3-branes. More precisely, our back- 
ground is the AdS black hole geometry discussed in section [2. 2. 11 which is the geometry dual to a field 
theory at finite temperature [22 . In this background, we encountered D7-brane embeddings of Minkowski 
type as well as black hole embeddings. The phase transition between both classes of embeddings is of 
first order [371 138) . The analysis of the meson spectrum shows that it corresponds to a transition between 
a phase of stable bound states of the fundamental degrees of freedom and a phase in which these mesons 
have finite lifetime. In physical parameters from the field theory point of view, we are in the deconfincd 
phase at high temperatures at which mesons are unstable and said to be melting. For a well defined 
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notion of high versus low temperatures, we need to compare the temperature to some energy scale. In 
our setup, the only available scale for comparison is the quark mass. Whether we are in the stable or 
in the melting phase, i. e. whether one or the other type of embedding is realized on the gravity side, 
depends on the ratio of quark mass to temperature T, c. f. figure 12.41 This can be seen from the 
equation of motion for the embeddings (j2.72l) . It is invariant under scale transformations by a factor 
of a, resulting in L i— > aL, w i— > aw and Tq h- aro. Scaling L by a amounts to scaling the quark mass 
ruq by the same factor a, while scaling of Tq is equivalent to scaling the temperature T by this factor. 
Because of the scale invariance of the equation of motion the functional behavior of the embedding, and 
therewith the physics of the D3/D7 setup, is identical for all setups with the same ratio of quark mass 
and temperature. From (j2.76p and (j2.57p we infer 

rn^^rn^VX^ (3.14) 
T To 2 

The free parameters of our setup appear on the right hand side as the asymptotic value TOj, of the D7- 
brane embedding and the black hole horizon Tq . The ratio of quark mass and temperature is defined by 
the ratio of these parameters, which we will henceforth denote by the dimensionless quantity 

2mq 

It was found that the transition to the melting meson phase occurs at a value of approximately m = 1.3 
[3 7) . At this value there is a change in the topology of the probe brane, which changes between the black 
hole type with a singularity and the regular Minkowski embedding. We demonstrated this in figure 

We use m as the parameter which defines whether we are in the regime of high or low temperature. 
However, this will not affect the topology of the embedding in our setup. Below, we follow an argumen- 
tation which reveals that we may restrict to black hole embeddings, since this is the thermodynamically 
favored configuration in setups with finite particle density. Black hole embeddings are conveniently 
described in the coordinate system (|A.2bp . derived in appendix \X\ and also used in ref. |58) . 

f Pie) ^,2 , ^/„^J„2^ , o2/^df?2 ^ ,„2 , (, „2N-2 2 



d.^ = ^ (^-^ dt^ + /(,) dx^ j +R^^^ + {l- + (1 - xT' dx^ + X^d02 j . (3.16) 

with 

hQ) = ^ + -,- (3-17) 

Q Q 

In the following, some equations may be written more conveniently in terms of the dimensionless radial 
coordinate p = q/to, which covers a range from = 1 at the event horizon to p — >■ oo, representing the 
boundary of AdSs space. 

As in section [2.2.2l we embed Nf D7-branes in this spacetime, such that they extend in all directions of 
AdSs space and along the directions of the three-sphere S'^, which is part of the S^. Due to the symmetries 
of this background, the embeddings depend only on the radial coordinate g and are parametrized by the 
functions xio)- Due to our choice of the gauge field fluctuations in the next subsection, the remaining 
three-sphere in this metric will not play a prominent role. The induced metric on the D7-brane is given 

by 

d^T = ^ (-J + /dx^) + f ^^f-^ V + R\l X') dnl (3.18) 

Here and in what follows we use a prime to denote a derivative with respect to g (resp. to p in dimen- 
sionless equations). We write yJ—G to denote the square root of the determinant of the induced metric 
on the D7-brane, which is given by 



,/,/ 



i± {l-X^) ^l-X^ + g^x'^. (3.19) 

Note that in general the branes are not necessarily coincident, and thus there will be one embedding 
function per brane, i. e. a = 1, . . . , Nf. We will make use of the DBI action to derive the embedding 
profiles x°'{q)- However, we postpone this task to the following subsection because the action also depends 
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on the gauge field on the probe brane. We wih see that the embedding function couples to this field on 
the brane. 

Each of the branes features a U(l) "flavor gauge field" A"-, with a — 1, . . . , Nf. This gauge field is 
arises from the fluctuation modes of an open string with both ends attached to the probe brane. For 
branes at arbitrary positions we therefore have an overall U(l)^-'' symmetry which is promoted to an 
U(A'^) in the case of coinciding branes. This symmetry enhancement comes from the fact that we then 
can no longer distinguish the branes and therefore cannot tell on which brane a string ends. Each of the 
two ends of a string can be assigned a label, also refered to as a Chan-Paton factor, which identifles the 
brane on which the string ends. There are N'j possible conflgurations, matching the degrees of freedom 
of the non-Abelian symmetry group U(A^/). The correct action for such a conflguration of coinciding 
branes is given by the non-Abelian DBI action [70] . 

S = -Tp jd^i sTr [det^Q det (P [E + E.,{Q-^ - 1)'^ E^.] + 27ra'F)] ^ (3.20) 

with • as a placeholder for a spacetime index, and 

Efiy = gfj,v + Bf^^, (3-21) 
Q'j = 5] + i2Tia' , ^'^j Ekj . (3.22) 

The Greek indices label the background spacetime coordinates, while Latin labels fc denote the 
directions perpendicular to the brane. Note that Q is a matrix with labels denoting solely these directions, 
and the determinant operator detj^ acts with respect to them. The determinant det acts with respect 
to the labels of the directions along the brane. Additionally to this Lorentz structure, all the operators 
in the action above are elements of the U(7V/) algebra on which the symmetrized trace sTr acts. The 
operators are operator valued analogon to the scalar embedding functions. 

The non-Abelian nature of the embedding functions $ introduces non-commutativity of the spacetime 
coordinates. The physical consequences of the non-Abelian DBI action are not entirely understood by 
now. However, non-commutative spacetimes are candidates for the cure of UV divergences of quantum 
field theories and are applied in M-theory to describe spherical D- brane configurations [7T]. The non- 
commutative contributions are manifest in p.22p . and hidden in the gauge covariant puUback of g, which 
introduces gauge covariant derivatives of the embedding, 

P[g]ab = Gab = Va X^" Vfc x"^ 

= gab + 5a»Vfc$' + 5b^Va«'' + gijVa^' Vfe*^' 

with 

Va$* = da^' + I 2Tia' [Aa , , (3.24) 

where the indices i,j are transverse and a, 6 are along the worldvolume of the probe brane. The non- 
Abelian DBI action features commutator terms [A , $] and [$ , <&] of gauge fields and embedding functions. 
These commutator terms can be thought of as corrections to the Abelian DBI action, which is reproduced 
if the commutators in (j3.22p and (I3.24p are vanishing. 

All setups we consider, feature a symmetry in the directions transverse to the D7-branes which allows 
to set one of the two to zero, i. e. the embedding function in this direction is constantly zero. Thus, 
the commutators in (j3.22p vanish. Moreover, we restrict to background configurations, arising from fields 
which are part of the Cartan subalgebra of lJ{Nf). As a justification for this restriction, we claim our 
freedom to define a basis in the vector space of the IJ (Nf) algebra and choose the non trivial embedding 
to define the direction of the generator 1, which is the generator of U(l) C U(iV/). The embedding 
matrices $ thereby are diagonal. A usual interpretation of the Eigenvalues on the diagonal is that they 
give the embedding functions for each of the Nf branes. The generators of the U(l) symmetry affect all 
flavor branes in an identical manner. The U(l) is therefore interpreted as the symmetry associated to 
baryon charge. 

Recall that the embeddings determine the quark masses of the dual field theory. The construction 
of coinciding branes in our setup therefore implies that the flavor eigenstates coincide with the mass 
eigenstates of the particles in the dual fleld theory. Generalizations to distinct bases in the flavor and 
mass vector spaces should be possible. 
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In the following, wc arc especially interested in two different setups, both of which feature a vanishing 
Kalb-Ramond field B = but different Nf. The case of A^/ = 1 allows investigation of the effects of 
baryon charge and the corresponding baryonic chemical potential. The symmetrized trace in the action 
is trivial in this case. The choice of Nf = 2 features the diagonal generator proportional to diag(l, —1) 
which can be interpreted to charge the two flavors oppositely and therefore models isospin symmetry. 
The non-Abelian DBI action simplifies to 

S = -Tpj sTr V|det(G + 27ra'i^)| . (3.25) 

According to the arguments above, the U(A^/) matrix structure of G and F for Nf = 2 is given by 

G = Ga°, (3.26) 
F = (7° + F'a^, (3.27) 

where we use the Cartan subalgebra of the U(2), given by two Pauli matrices 



^' = [o 1 j' -' = io -1 )■ 

All operators in the action therefore are diagonal. The two probe brancs are coincident and the diagonal 
entries of the field strength tensor F determine the net charges of the branes. We have 

p_(F^+F' \_[fW \ 

F^-F' )-[ ) ■ (^-29) 

In any case, the restriction to the diagonal Cartan subalgebra of V{Nf) simplifies the non-Abelian 
DBI action, e. g. symmetrization of the trace is trivial in the sense that all commutators vanish and 
sTr = Tr. Expansion of the square root and evaluation of the trace with subsequent restoration of 
the square root eventually leads to the following action for two D7-branes in the AdS^ x S^black hole 
background with vanishing B field, 

Nf 

S = -TrY, I ^det(G + 27raW). (3.30) 

k=l ■' 

We will concentrate on the cases of Nf = 1 and Nf = 2, and separately switch on either the baryonic 
C/(l), parametrized by field strengths F^, or the isospin subgroup SU(2) along the direction of a^, 
parametrized by the field strengths F^ . Thus, the action for each brane is the same as long as we do not 
switch on both fields simultaneously. 

To sum up, the background geometry described so far is dual to thermal M = A supersymmetric 
SU(A'') Yang-Mills theory with Nf additional Af = 2 hypermultiplets. These hypermultiplets arise from 
the lowest excitations of the strings stretching between the D7-branes and the background-generating 
D3-branes. The particles represented by the fundamental fields of the Af = 2 hypermultiplets model the 
quarks in our system. Their mass rTi,^ is givcm by the asymptotic value of the separation of the D3- and 
D7-branes. Since the physics of the thermal D3/D7 setup is determined by the ratio of quark mass to 
temperature, we use the parameter m, which is proportional to this ratio. 



3.2.2 Background gauge fields — finite particle density 

In addition to the parameter m, we aim for a description of the system at finite baryon density n^,, 
which in turn accounts for a finite chemical potential /x. In the thermal SU(A^) gauge theory, a baryon 
is composed of N quarks, such that the baryon density Ub can be directly translated into a quark 
density Uq = Ug N . The thermodynamic dual quantity of the quark dc^nsity is the quark chemical 
potential fi, which is realized by a non-dynamical time component of the gauge field. The chemical 
potential is the source of the charge density operator Jg, i. e. the time component of the current J, of the 
particles charged with respect to the potential under consideration. The time component of the current 
of fundamental spin fields cj) and the spin 1/2 fields tp in the Af = 2 hypermultiplet is given by 

Jo = ^lotp + 4>do4>- (3.31) 
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In the dual holographic formulation the source fi of this charge density then corresponds to the non- 
rcnormalizable mode of the according holographically dual field in the supergravity theory. This field is 
the time component of the supergravity gauge field on the probe brane, which we denote by ^o- The 
normalizable mode will yield the expectation value of particle density. We consider a constant chemical 
potential in space and time, i. e. there is no spacetime dependence. Instead we work with a gauge field 
background Ao{p) which only depends on the radial AdScoordinate. 

We do not rederive the dictionary entries here, but rather rephrase what is important for the following 
developments. The holographic interpretation of the embedding x of the probe D7-branes was discussed 
(in a different coordinate system) in section 12.2.21 The probe branes account for holographic duals 
of fundamental quarks with mass niq, determined by the non-renormalizable mode of the embedding 
function. The asymptotic form of the fields xip) ^-^id Aq{p) can be found from the equations of motion 
in the boundary limit p — ?> cx). The expansion coefficients in an expansion in powers of p are given by 

Ao = fi,-^- \ ' +■■-, (3.32) 

X=- + ^ + ---. (3.33) 
P P^ 

Here fj,q is the quark chemical potential, riq is the quark density, m is the dimensionless quark mass 
parameter given in p.lSp and Cc is the quark condensate mentioned earlier (but irrelevant in this work) . 
We made use of the dimensionless p-coordinate that runs from the horizon value p = 1 to the boundary 
at p (X). The chemical potential and density of baryons are simply 

..^f, n. = ^. (3.34) 

Once we have found the solutions Aq to the equations of motion for the gauge field, the value fiq of 
the chemical potential in the dual field theory can be extracted as 

Pq = lim Ao{p) = TT^Pq, (3-35) 
p-i-oo lira 

where we introduced the dimensionless quantity /t for convenience. We apply the same normalization to 
the gauge field and distinguish the dimensionful quantity A from the dimensionless 

-- 2T:a' - 

Ao = Ao 3.36 

(we save the symbols without diacritics for later use). Analogously, the solutions of the embedding 
functions carry information about the quark mass parameter m, 

m — lim px{p)- (3.37) 

p— f oo 

We mentioned that for non-vanishing baryon density, there are no embeddings of Minkowski type, 
and all embeddings reach the black hole horizon. This is due to the fact that a finite baryon density in 
an infinite volume of Minkowski spacetime requires an infinite number of strings in the dual supergravity 
picture. These strings have one end on the stack of D3-branes and the other on the stack of Nf probe 
D7-branes. These strings pull the brane towards the black hole 58j. Such spike configurations are 
common for configurations in which branes of different dimensionality connect |33j . 

Very recently, however, it was found that for a vanishing baryon number density, there may indeed be 
Minkowski embeddings if a constant vacuum expectation value of is present, which does not depend 
on the holographic coordinate [SDl IMl IZH IZSl IZl]- The phase diagram found there is reproduced in 
figure [33] In the shaded region, the baryon density vanishes (ne — 0) but temperature, quark mass and 
chemical potential can be nonzero. This low temperature region only supports Minkowski embeddings 
with the brane ending before reaching the horizon. In contrast, the unshaded region supports black hole 
embeddings with the branes ending on the black hole horizon. In this regime the baryon density does not 
vanish (jib 7^ 0). At the I0W-/1 end of the line separating Ug — Q from 7^ in figure [231 there exists 
also a small region of multivalued embeddings, which are thermodynamically unstable [59) . In the black 
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Figure 3.1: The phase diagram of fundamental matter in the D3/D7 setup. Horizontal axis: Chemical 
potential normalized to the quark mass. Vertical axis oc T/rUq. In this work we analyze the white 
region of finite particle density d, for which we show some lines of constant values for d. The relation 
between d and is given by p.43p 

hole phase there is a phase transition between different black hole embeddings , resembling the meson 
melting phase transition for fundamental matter at vanishing density. This first order transition occurs 
in a region of the phase diagram close to the separation line between the two regions with vanishing 
(shaded) and non- vanishing (unshaded) baryon density. This transition disappears above a critical value 
for the baryon density rig given by 



In this work we exclusively explore the region in which > 0, i. e. we examine thermal systems in 
the canonical ensemble. For a detailed discussion of this aspect see refs. [501 [5U] . 

To determine the solutions of the supergravity fields on the probe branes we have to extremize the 
DBI action (j3.30L we write shortly as 




The induced metric G(^) on the stack of Nf coincident branes is given by p.lSp . F is the dimensionless 
field strength tensor of the gauge fields on the brane. 

For now we consider the simpler case of a baryonic chemical potential modeled by the U(l) subgroup 
of U{Nf). In this case, the sum amounts to an overall factor of Nf. In ref. [5^ the dynamics of such a 
system of branes and gauge fields was analyzed in view of describing phase transitions at finite baryon 
density. Here, we use these results as a starting point which gives the background configuration of the 
probe branes' embedding function and the gauge field values at finite baryon density. To examine vector 
meson spectra, we will then investigate the dynamics of fluctuations in this gauge field background. 

In the coordinates introduced above, the action Sdbi for the embedding xis) ^-nd the field strength 
F is obtained by inserting the induced metric and the field strength tensor into (j3.39l) . From now on 
we make use of the dimensionless coordinates and reproduce the action found in ref. [SS]- To do so, 
we remember that the only non-vanishing component of the background field is the p-dependent time 

~ (k) " (k) 

component. Therefore, the only non- vanishing components of the field strength tensor are F^q = —Fq^ . 
We evaluate the determinant and arrive at 



25/2 



with critical d* = 0.00315 . 



(3.38) 





(3.39) 




(3.40) 



with components G"*" of the inverse metric G ^. After inserting these components we get 



(3.41) 
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where F40 = dpAo is the field strength on the brane. The background fields x and Aq depend solely on 
p. This action only depends on derivatives of the gauge field. We therefore can identify the constant of 
motion d satisfying dpd = 0, 

d [dpAo) 

Evaluation of this formula and insertion of the asymptotic expansion of A reveal that this dimensionless 
constant is related to the parameters of our setup by |58j 



25/2, 



Nfy/XT^' 



(3.43) 



We can therefore think of the constant d as parametrizing the baryon density . 

The equations of motion for the background fields are conveniently obtained after Legendre trans- 
forming the action p.4ip to S = S — SS/6F40 in order to eliminate dependence on the gauge field in 
favor of dependence on the constant of motion d |58j . Varying this Legendre transformed action with 
respect to the field x gives the equation of motion for the embeddings x{p)j 



v/1 - + P^X''' 




^6/3(1 _ ^2)3 



8d2 



^6/3(1 _ ^2)3 



,2,y2 



3(1 - x') + 2pV' - 24d2 ^ X +PX 



p6J3(l_^2)3+8^2 



(3.44) 



This equation for xip) can be solved numerically for given d and initial value xo- We impose boundary 
conditions such that the branes cross the horizon perpendicularly 



xip = 1) = Xo, 



dpXip) 



p=i 



0. 



(3.45) 



Figure [3?2l shows some examples. The embeddings at finite density resemble the large p asymptotics of 
the embeddings found at zero density. For small p however, at finite particle density there always is the 
spike reaching down to the event horizon. The initial value of xo determines the position on which the 
brane reaches the horizon and in this way determines the quark mass parameter m, c. f. equation p. 371) . 
It is zero for Xo = and tends to infinity for xo — ^ 1 ■ Figure 13.31 shows this dependence of m on xo for 
different values of the baryon density d. In general, a small (large) xo is equivalent to a small (large) m. 
Foi' Xo ^ 0-5j "we nearly observe proportionality. For vanishing d = 0, we only can model quarks with 
m < 1.3, heavier quarks at vanishing density are described by embeddings of Minkowski type, which 
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we do not discuss here. (The trained eye can see that there is a maximum of m = 1.3 in figure 
before m drops to smaller values towards xo = 1- This is reflecting the existence of a phase transition 
to Minkowski embeddings at d = 0). 

The equation of motion for the background gauge field A is given by 



dpAo = 2d 



/(i-x") p^pii~x'r + 8d^ 



(3.46) 



Integrating both sides of the equation of motion from pu = 1 to some p, and respecting the boundary 
condition Ao{p — 1) =0 [58], we obtain the full background gauge field 



fVi-x' + P^x'' 



Aoip)^2d I dp 

1 W/(1-X') P^Pil~x')' + 8d^ 



(3.47) 



Examples for the functional behavior of Ao{p) are shown in figure [331 While there is a significant slope 
of A{p) near the horizon at p = 1, the gauge field tends to a constant at large p. From p.35p we recall 
that this value is the chemical potential of the field theory. We will henceforth compute p by evaluating 
the formula above for large p. Note that at any finite baryon density d cx ^ there exists a minimal 
chemical potential which is reached in the limit of massless quarks. 



3.3 Meson spectra at finite baryon density 
3.3.1 Equations of motion 

We now compute the spectral functions of flavor currents at finite baryon density d, and temperature 
T in the black hole phase. Compared to the limit of vanishing density treated in [35], we discover 
a qualitatively different behavior of the finite temperature excitations corresponding to vector meson 
resonances. 

To obtain the mesonic spectral functions, we compute the correlations of flavor currents J by means 
of the holographically dual gauge field fiuctuations about the background given by (|3.4ip . We denote 
the full gauge field by 

i^(p,x) =<50io(p) + A^(f,p). (3.48) 

According to section 13.21 the background field has a non- vanishing time component, which depends 
solely on p. The fluctuations in turn are gauged to have non- vanishing components along the Minkowski 
coordinates x only, and only depend on these coordinates and on p. Additionally, the fiuctuations are 
assumed to be small, such that it suffices to consider their linearized equations of motion. At this 
point we simply neglect the fluctuation of the scalar and pseudoscalar modes and their coupling to the 
vector fluctuations. This procedure is justified by the restriction to fluctuations with vanishing spatial 
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p d 

Figure 3.4: Examples for the background gauge field time component Aq (left) and the resulting chemical 
potential fl (right). Note that the chemical potential is not zero at asymptotically small but non-zero d 
for m > 1.3, reproducing the phase transition line in figure [5711 



momentum, which is imposed later. In this limit the vector mesons do not couple to the other mesonic 
excitations. A generalization of this work which includes these coupling and finite momentum spectra 
can be found in ref. [63] . 

The equations of motion are obtained from the action p.39p . where we introduce small fiuctuations 
A by setting 



A^ A + A, 
Fi-^ F + F. 



(3.49) 
(3.50) 



The background gauge field A is given by p.46p . The fiuctuations A now propagate on a background Q 
given by 

g = G + F, (3.51) 
and their dynamics is determined by the Lagrangian 



Vidct ig + F)\, 



(3.52) 



with the fiuctuation field strength — 29[^ A^j . In contrast background field A, the fluctuations A do 
depend on the Minkowski directions as well as on the radial coordinate of AdSs. Since the fluctuations 
and their derivatives are chosen to be small, we consider their equations of motion only up to linear 
order, and derive these equations from the part of the Lagrangian L which is quadratic in the fields and 
their derivatives. Denoting this part by L2, we get 



L2 = -i^d^ (g'^'^g^^F^pF^^ ~ ^ g'^'^g^^F^^Fp^^ 



(3.53) 



Here and below we use upper indices on g to denote elements of g ^ . The equations of motion for the 
components of A are 



,/\d^\[gt--g''^ - g^^-g-^ - g^'^^^g^^) d^^A^,] 



(3.54) 



The terms of the corresponding on-shell action at the p-boundaries are (with p as an index for the 
coordinate p, not summed) 



on-shcll 



ToTT^R^NfTr J d^x^y\detg\ 



Ao9pAo - g'^'^g^'AdpAk - Aog^'hrig-'F) 



PB 



Ph = 1 



(3.55) 



41 



From this form of the action we can derive the correlation functions by means of the procedure 
outhned in section [2.2.1l First, we Fourier transform the fields as 



A, [P, ^) = /(04 e^'^^^M (P: k). (3.56) 

As above, we are free to choose our coordinate system to give us a momentum vector of the fluctuation 
with non- vanishing spatial momentum only in x-direction, k = (w, g, 0, 0). 

To obtain the correlator Gq. with indices i, k labeling Minkowski directions, we have to consider the 
second term in the parentheses of p.55p . including all its prefactors. Denote the resulting expression by 
A{p, k). We decompose the gauge field fluctuations into a boundary and a bulk contribution, A(p, k) = 
f{p, k)A^'^y{k), where limp_j.oo f{p, k) = 1. The prescription from section [2.2. II tells us to divide out the 
boundary terms A^'^'^ from A{p, k) in order to get what was denoted by F in (|2.6ip . Once we found the 
solutions A{p, k) we can obtain this expression by evaluating 

F,,[p, k) ^ hm -^if ' ^1 , (3.57) 

p, p^^ A,{p,k)Ak{p,k) 

where the indices i, k correspond to the Minkowski indices on Gfj.. 

To evaluate this expression, we have to insert the solutions to the equations of motion for A(p, k) 
into this expression. Note that on the boundary at p — >■ cx), the background matrix Q reduces to 
the induced D7-brane metric G. Therefore, the analytic expression for boundary contributions to the 
on-shell action is identical to the one found in ref. |48| . In our case of finite baryon density, new features 
arise through the modified embedding and gauge field background, which enter the equations of motion 
(I3.54P for the field fiuctuations. 

We adopt the procedure of ref. |3H] i where the coordinates in Minkowski directions where chosen such 
that the fluctuation four vector k exhibits only one non vanishing spatial component, e. g. in z-direction 
as fc = {uj, q, 0, 0). In addition, the action was expressed in terms of the gauge invariant field component 
combinations 

Ex^i^Ax + qAa, Ey^z^LoAy^^. (3.58) 

In the case of vanishing spatial momentum g ^> 0, the Green functions for the different components 
coincide and were computed as |48) 

G« ^ G^^ = G« = G.". = ^MsZ! ( ,3 MM) , (3.59) 



E{p) J ' 



where the E{p) in the denominator divides out the boundary value of the field in the limit of large 
p. Again, the indices on the Green function denote the components of the operators in the correlation 
function, all off-diagonal correlations (as Gyz, for example) vanish. 

In the limit of q — > 0, the equations of motion for transverse fiuctuations Ey^^ match those for 
longitudinal fiuctuations E^- For a more detailed discussion see ref. |3H]- As an example, consider the 
equation of motion obtained from p.54p with cr = 2, determining Ey = UJA2, 

goo 



0=E" + dp In ( ^\de^\G''G'^ ) E' - ^ r^^^E 



Here we introduced the dimensionless frequency 

In order to numerically integrate the equations of motion (I3.60p . we determine local solutions of that 
equation near the horizon at p^^ = 1, which obey the infalling wave boundary condition. This condition 
ensures causality by demanding that the excitations can propagate in inward direction, but nothing can 
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Figure 3.5: An example for a spectral funetion at finite baryon density, compared to the zero temperature 
result. 



exit the horizon. The local solutions can be used to compute initial values in order to integrate p.60p 
forward towards the boundary. The equation of motion p.60p has coefficients which are singular at 
the horizon. According to mathematical standard methods, the local solution of this equation behaves 
as [p — PhY 1 where /? is a so-called 'index' of the differential equation 75 . We compute the possible 
indices to be 

P = ±ixo. (3.62) 

Only the negative sign will be retained in the following, since it casts the solutions into the physically 
relevant incoming waves at the horizon and therefore satisfies the incoming wave boundary condition. 
The solution E can be split into two factors, which are {p — 1)"*" and some function F{p), which is 
regular at the horizon. The first coefficients of a series expansion of F{p) can be found recursively as 
described in EH [S51 ■ At the horizon the local solution then reads 



i?(p) = (p-ir- 
= (p~i)-"- 



So, F{p) asymptotically assumes values 



F{P) 



1 + ^{P-1) + 



(3.63) 



F(p = l) = l, 



dpF{p) 



p=i 



(3.64) 



To calculate numeric values for E{p)^ we have to specify the baryon density d and the initial value 
XOi which determines the mass parameter m. These parameters determine the embeddings x appearing 
in p.60p . We can then obtain a solution E for a given frequency fo by numerical integration of the 
equation of motion (j3.60l) . using the initial values (|3.63p and p.64p . 

Spectral functions are finally obtained by combining (j3.59p and (13. 8p . 



lH(a;,0) 



NfNcT^ 



Im lim 



E{p) ) 



(3.65) 



3.3.2 Spectra 

We now discuss the resulting spectral functions at finite baryon density, and observe crucial qualitative 
differences compared to the case of vanishing baryon density. In figure 13.51 an example for the spectral 
function at fixed baryon density oc d is shown. In the limit of large tn, corresponding to asymptotically 
small temperatures, the spectral function can be derived analytically. This zero temperature result is 
given by 

JHo = TV/TVcT^ Trro^ (3.66) 

Figure 1531 shows this function as well. 
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All graphs shown here are obtained for a value of d above d* , given by p.38p , such that we investigate 
the regime in which there is no fundamental phase transition of first order. Recall that the parameters 
of our theory are given by d cx Ub/T^ and m cx raq/T . Therefore variations in the quark density at fixed 
temperature and quark mass are introduced by tuning d only. The effects of different quark masses can 
be seen by tuning m alone. The effect of changes in temperature involves changes in both ni and d. 

It is interesting to compare the spectra we obtain at finite temperature and density to the vector 
meson spectrum obtained at zero temperature and vanishing quark density. It is given by the same 
relation as the mass spectrum (|2.81l) which we encountered in the example of scalar mesons [33]. In 
our case, where the mesons do not carry spatial momentum, we can translate the mass M„ of the n**^ 
excitation into an energy a;„ = A/„. At this energy we would see a resonance in a supersymmetric setup. 
In terms of the dimensionless quantities we use here, these resonance energies are given by 



M„ / (n + l)(n + 2) 

2^^\ 2 ' " = 0,1,2,..., (3.67) 



where n labels the Kaluza-Klein modes arising from the D7-brane wrapping the S 



3 



Finite temperature effects 

We analyze finite temperature effects by choosing two values of m and d, which correspond to a given 
values of quark mass, quark density and a temperature T. A change in temperature amounts to 

T^aT (3.68) 

and thereby leads to 

~, d m ,„ ^. 

d^—, m^—. (3.69) 
a'^ a 

An example is shown in figure [3^ There we plot spectra for three different temperatures, which we call 
low (m = 5, d = 1), medium (m = 5/2, d = 1/8) and high (m = 5/3, d — 1/27) temperature. We can 
see that at high temperature there is hardly any structure visible in the spectral function. However, we 
have chosen a temperature at which already a slight excitation is visible at low energies to. Decreasing 
temperature leads to more and more pronounced peaks in the spectral function. Moreover, at decreasing 
temperature these peaks move closer to the resonance energies l|3.67p , corresponding to zero temperature 
and density (drawn as the corresponding dashed lines in the figure). 

The formation of sharp resonances at low temperature indicates the intuitively expected behavior of 
long living mesons in a cold medium, which melt, i. e. decay faster, at high temperatures. However, we 
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Figure 3.7: The dependence of the spectra on baryon density. The dashed lines again mark the super- 
symmetric spectrum. 

did not perform an analysis of the quality factor of the resonance peaks, i. e. we did not calculate the 
lifetimes of the vector mesons. From figure 13.61 we can see that the height-to-width ratio of the peaks 
seems not to improve to a great extent at low temperatures. 

Finite density effects 

To investigate the effects of finite baryon density n^, we tune d while keeping m constant. This amounts 
to varying the quark density at constant temperature and quark mass. The effect is shown in figure [3771 
We observe that the peak width is considerably influenced by baryon density. At low baryon density the 
resonances are close to line-like excitations, while they are broadened with increasing particle density. 
Additionally, increasing the particle density also causes a slight shift of the resonances to higher energies. 

These observations are interesting from a phenomenologically inclined point of view. The in-medium 
effects on mesonic bound states are important to interpret the results of heavy ion collision experiments. 
Estimations from the early 1990s based on effective models predicted decreasing vector meson masses at 
increasing densities [76 , known as Brown- Rho scaling. Experimental data from experiments at the SPS 
facility at CERN, however, is better described by models like the one found in refs. [771 [7S]. There the 
in-medium effects also are refiected in peak broadening and shifts to higher energies. 

For information on the spectral functions at vanishing particle density we refer the reader to ref. [48] . 
Where the low temperature regime for d = was investigated. 

Dependence on quark mass 

To observe the dependence on the mass of the quarks, we plotted spectra for different m at constant 
values of d in figure 13.81 We observe more and more pronounced resonances as we increase the meson 
mass. These mesons eventually nearly resemble the line spectrum p.67p known from the supersymmetric 
case of zero temperature and vanishing quark density. This observation refiects the decreasing effect of 
finite temperature and chemical potential with increasing quark mass. In a regime where the scale of the 
quark mass outweighs both additional scales T and d their effects seem to be negligible. This is the case 
when we observe a configuration which is located close to the Minkowski phase in the phase diagram, cf. 
figure 13.11 

In ref. |4] we elaborate on the spectral functions behavior at low quark masses. There we observed 
that the position of the vector meson excitations in the regime of very low quark masses decreased 
with increasing quark mass. Further increasing the quark mass lead to increasing quark masses as 
described in this section. We omit this discussion here, but resume on the topic when we discuss the pole 
structure of the spectral functions. The reason is that the peaks referred to in ref. ^ are only visible 
after subtraction of the zero temperature part *Ho from the spectral function. To interpret the spectral 
function as a probability density for the detection of a quasiparticle, we cannot subtract UIq, as we would 
otherwise produce negative probability densities, which are not well defined. 
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Figure 3.8: The dependence of the spectra on quark mass. The dashed lines again mark the supersym- 
metric spectrum. 



Our setup is a modification of the one used in ref. There, the authors considered vector meson 
spectra at vanishing baryon density. These spectra only show peaks moving to smaller frequency as the 
quark mass is increased. There is no contradiction to the results presented this work. Note that the 
authors of ref. |3H] by construction are restricted to the regime of high temperature/small quark mass. 
Nevertheless, they continue to consider black hole embeddings below the temperature of the fundamental 
phase transition where these embeddings are only metastable, the Minkowski embeddings being thermo- 
dynamically favored. At small baryon density and small m our spectra are virtually coincident with 
those of [48]. However, in our case, at finite baryon density, black hole embeddings are favored for all 
values of the mass over temperature ratio. 

3.3.3 Pole structure 

In this section we comment on the quasi normal modes (QNM) of the system under investigation. As 
discussed above, these are the poles in the complex frequency plane, i. e. the spectral function diverges 
at these locations. An impression of the continuation of the spectral function into the complex frequency 
plane is given in figure 13.91 The data to plot the graph was obtained in exactly the same way as for the 
spectral functions shown in the preceding sections, except the fact that we upgraded the numerics to 
process ir G (D. Therefore the spectral functions shown so far are given by the values along the real axis. 

Our numerics turn out to be reliable for | Im Id | < 1 and therefore cannot determine poles in the plane 
of tn G (D which lie beyond this limit. We trust the values within the regions shown in the figures of this 
work, although there possibly is room for improvement in accuracy. We checked our code for stability 
against the initial conditions and parameters, and are mainly interested in the qualitative behavior of 
the results. 

The spectra presented in preceding sections show that the first resonance peaks, i. e. those for small 
n, are very narrow, while the following peaks show a broadening accompanied by decreasing amplitude. 
The physical consequence would be a longer lifetime of the lower n excitations. This is refiected in a 
smaller imaginary part of the corresponding quasi normal mode in figure 13.91 It is a known fact that 
the quasinormal modes develop larger real and imaginary parts at higher n. So the sharp resonances 
at low It), which correspond to quasiparticles of long lifetime, originate from poles with small imaginary 
part. For higher excitations in n at larger ro, the resonances broaden and get damped due to larger 
imaginary parts of the corresponding quasi normal modes. 

Above we observed variations in the positions and widths of the peaks in the spectral function, 
depending on the changes in temperature, particle density and meson mass. This behavior can be 
translated into a movement of the quasinormal modes in the complex plane. Figure 13.101 shows the 
trajectory of the quasi normal mode corresponding to the first peak in the spectral function, parametrized 
by TO. At small densities, we can see the turning behavior of the mode which starts to move in the direction 
of decreasing real part at small to, and then turns to asymptotically large real part while converging to 
the real axis for further increasing to. 
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Figure 3.9: A contour plot of the spectral func- 
tion in the complex tr-plane. Blue shading in- 
dicates 91 > 0, brown shading indicates Dl < 0. 
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cal spectrum we plotted for several parameters 
above. 
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It would be interesting to compare our results to a direct calculation of the quasinormal modes 
of vector fluctuations in analogy to ref. [44 . There, the quasinormal modes are considered for scalar 
fluctuations exclusively, at vanishing baryon density. The authors observe that starting from the massless 
case, the real part of the quasinormal frequencies increases with the quark mass first, and then turns 
around to decrease. This behavior agrees with the peak movement for scalar spectral functions observed 
in ref. ^S] (above the fundamental phase transition at xo < 0.94) where the scalar meson resonances 
move to higher frequency first, turn around and move to smaller frequency increasing the mass further. 
These results do not contradict the present work since we consider vector modes exclusively. 

3.4 Meson spectra at finite isospin density 
3.4.1 Equations of motion 

In order to examine the case of two flavors, Aff = 2, with opposite isospin chemical potential in the 
strongly coupled plasma, we extend our previous analysis of vector meson spectral functions to a chemical 
potential with SU(2)-flavor, i. e. isospin, structure. Starting from the general action p.25p we now 
consider the non-Abclian field strength tensors 

with the Pauli matrices cr" and A given by equation p.48p . The upper index on the gauge field labels the 
component in the vector space of the SU(2) generators. The factor r^/{2Tra') is due to the introduction 
of dimensionless fields as described below (|3.35p . The totally antisymmetric f'^'"^ = e"'"^ with e^"^^ = 1 
arise from the structure constants of SU(2). 

In the non-Abelian field strength tensor, the term quadratic in the gauge field describes a self inter- 
action of the gauge field. The coupling constant for this interaction may be determined by a redefinition 
of the gauge field, such that the kinetic term of the effective four-dimensional theory has the canonical 
form. In appendix [C] (taken from ref. [3]) we show that the redefinition is given by 

A^^A, (3.71) 
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where the dimensionless constant Ca depends on the geometry of the D7 worldvolume directions along 
p and the S'^, which are transverse to the directions of the D3-brane. In particular, is independent of 
the 't Hooft coupling A. Determining the exact value of is left to further work in terms of the ideas 
presented in appendix [C] In the following we chose a convenient = The field strength tensor in 
the redefined fields is given by 



Ca 



= '^d^.K] + ^r'^A^Al (3.72) 

In order to obtain a finite isospin-charge density and its conjugate chemical potential /ij, we intro- 
duce a /o-dependent time component of the SU(2) valued background gauge field A [5J. This background 
field is defining a direction in the vector space of the SU(2) generators. We choose coordinates such that 
the direction of the background field aligns with the cr^ direction while the other SU(2) components are 
vanishing, 

^ = ^o^' = ^o(p)( J ), ^0 = ^0 = 0. (3.73) 

This specific choice of the 3-direction in flavor space as well as spacetime dependence simplifies the 
isospin background field strength, such that we get two copies of the baryonic background Fp^ on the 
diagonal of the flavor matrix, 

_ / dpAo 



The derivation of the background fleld conflguration leads to the same explicit form of the action as p. 411) . 
We can therefore make use of the background field solutions x{p) ^nd Aq{p) found in the baryonic case. 
As before, we collect the induced metric G and the background field strength F in the background tensor 

g = G + F. 

For the fluctuations, however, we encounter an additional structure. The SU(2) valued fluctuations 
in general have components along all the directions of this vector space. We make use of Tr(CTV^ ) = 2^'^ 
and apply the background field method in analogy to the baryonic case examined in section 13.31 As 
before, we obtain the quadratic action for the fluctuations by expanding the determinant and square 
root in powers of A° . The term linear in fluctuations again vanishes by the equation of motion for the 
background field. This leaves the quadratic action 

oo 

4o = 2^'^' TnTr Jdpd^x ^/\detg\ 
1 

G'-'^g''' {d^,A:^ + ^Alfja,,jac, ^.^^^^^ ^cj^^^ 

+ {g^-g^P _ g'^^gP'^) A_ Alr>^^ a^s^^. 

The factor Tr arises from the trace over the generators of SU(2). If we use the Pauli matrices as generators 
we get = 2. Another common choice for the generators is ct'^/'2,, which amounts to Tr = 1/2. We 
leave the explicit choice open, since it merely introduces an unimportant finite proportionality constant 
to the action. Note that besides the familiar Maxwell term, two other terms appear due to the non- 
Abelian structure. One of the new terms depends linearly, the other quadratically on the background 
gauge field A and both contribute non-trivially to the dynamics. The equation of motion for gauge field 
fluctuations on the D7-brane is 



(3.75) 



= 9. 



.J\d^\ {g^-g^^ - g'^-g^^) f^, - ij?/""' (^'^"a"^ - g-'^g''^) f' , (3.76) 



with the modified field strength linear in fiuctuations F^^ = '^^[p.K] +Ca/v^/"''^^o (Sop^K + SouAD r^/{2na') 
Integration by parts of (j3.75p and application of (I3.76P yields the on-shell action 



^on-shcll 



roT^Tr n^R^ Jd'^x ^Jdetg] x {C^G^'' - G^'^G^^) A^F'^^^p^ . (3.77) 
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The three field equations of motion (flavor index a — 1,2,3) for fluctuations in transverse Lorentz- 
directions /i = 2,3 can again be written in terms of the combination Ef^ = qA"^ + wA^. At vanishing 
spatial momentum g = we get 



Q = E, 



^" + 9pln(/^fetG| 



/ ^2 
/ 'o 



V27ra' 



-A? 



'q4Aq22 



G""(roa;)2 



2iroLj- 



2'Ka' 



43 p2 



(3.78) 



Q^E. 



T ^dpln[^\d^\G^^G^^j E^' 



\2-Ka' 



^0 



= i;|," + 9pln(V|detG| 



E', 



G^^G 



Ej^ 



' G44 

^AlE^ 



(3.79) 



G44 



(3.80) 



Note that we use the dimensionless background gauge field = ^o(27ra')/r-o with Tq = T-kB? . Despite 
the presence of the new non-Abelian terms, at vanishing spatial momentum the equations of motion for 
longitudinal fluctuations i?£ acquire the same form as the transverse equations p.78p . to p.80p . 

Two of the above ordinary second order differential equations are coupled through their flavor struc- 
ture. Decoupling can be achieved af@ in ref. by transformation to the flavor combinations 



X = E^^iE'^, Y = E^ 
The equations of motion for these fields are given by 



iE' 



0= + In (7^^544^22 



X' 

y 



0- y" + apin[v^d^^44^22 

- E^" + d, In f/M^e''^'' ) E 



r4 



j,4 gm 



m)'x, 



m)'r. 



(3.81) 

(3.82) 
(3.83) 



„4 nm 
A 'o » rti2p3 



with ro — Ld / (2-kT) and dimensionless but p-dependent m 
section 



(3.84) 

A^J{2'itT). Proceeding as described in 



we determine the local solution of these equations at the horizon. The indices turn out to be 

Ai{p^iy 



/3 = ±i tr T 



2ttT 



(3.85) 



Since ^q(p = 1) = we are left with the same index as in p. 621) for the baryon case, i.e. the chemical 
potential does not influence the singular behavior of the fluctuations at the horizon. The local solution 
coincides to linear order with the baryonic solution given in p.63p . 

For the special case of zero temperature the background geometry is AdSs x S^. For finite chemical 
potential in the zero temperature case we can obtain the gauge field correlators in analogy to ref. [19] . 
The zero temperature result IHo^so analog to (|3.66p is given by 



5Ko,iso = TR7Vr2^(lr±moo)' 



(3.86) 



with the dimensionless chemical potential nXoo = lim 



m. 



^At this point there is an essential difference which distinguishes this setup from the approach with a constant poten- 
tial j4q at vanishing mass followed e.g. in ref. [5]. While the metric coefficients for massless quarks are identical in both 
cases, there is a p-dependence of the background gauge field in the present setup. 
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Figure 3.11: The vector meson spectral functions of the three isospin components. For a concise image 
we emphasize the first peak of each component by stronger color saturation. 



3.4.2 Spectra 

Application of the recipe analog to the case of baryonic chemical potential yields the spectral functions 
of flavor current correlators in a medium with finite isospin density. Note that after transforming to 
fiavor combinations X and Y, given in p.8ip . the diagonal elements of the propagation submatrix in 
flavor-transverse X, Y directions vanish, Gxx = Gyy = 0. Now the off-diagonal elements give non 
vanishing contributions. However, the component E^, longitudinal in flavor space, is not influenced by 
the isospin chemical potential, such that Ge^e^ is nonzero, while other combinations with vanish [5]. 

In figure [5". Ill we compare spectral functions for the isospin case, where we emphasize the first peak 
of each of the three components. Note that the E'^ spectrum coincides with the baryonic case, as the 
equation of motion p. 841) coincides with p. 601) . 

While the qualitative behavior of the isospin spectral functions agrees with the one of the baryonic 
spectral functions, there nevertheless is a quantitative difference for the fiavor-transverse components 
X, Y . We find that the propagator for flavor combinations Gyx exhibits a spectral function for which 
the peaks are shifted to higher frequencies, compared to the Abelian case curve. For the spectral function 
computed from Gxy, the opposite is true, its peaks appear at lower frequencies. The quasiparticle 
resonance peak in the spectral function ^Kyx appears at higher frequencies than expected from the vector 
meson mass formula p.67p . The other flavor-transverse spectral function '^Kxy displays a resonance at 
lower frequency than observed in the baryonic case. 

This may be interpreted as a splitting of the resonance peak into three distinct peaks. This is due to 
the fact that we explicitly break the symmetry in flavor space by our choice of the background field . 
Decreasing the chemical potential reduces the distance of the two outer resonance peaks from the one in 
the middle and therefore the splitting is reduced. 

The described behavior resembles the mass splitting of mesons in presence of a isospin chemical 
potential expected to occur in QCD [THl HO]- A linear dependence of the separation of the peaks on the 
chemical potential is expected. Our observations confirm this behavior. Since the vector mesons are 
isospin triplets and we break isospin symmetry explicitly, we see that in this respect our model is in 
qualitative agreement with effective QCD models. Note also the complementary discussion of this point 
in ref. [SIl. 



3.5 Summary 

Two distinct setups were examined at non-zero charge density in the black hole phase. First, switching 
on a baryon chemical potential through non-zero baryon density, we find that nearly stable vector mesons 
exist close to the transition line to the Minkowski phase. Far from this regime, at small quark masses or 
high density, the spectral functions do not show distinct resonance peaks. 

Moreover, at small quark masses and particle densities we observe that the quasi normal modes move 
to positions with smaller real part in the complex w plane, in accordance with the observations in the case 
of vanishing chemical potential [IS] . Increasing the quark mass over temperature ratio beyond a distinct 
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value, the plasma adopts the behavior known from the case of zero temperature. In the spectral functions 
we computed, this zero-temperature-like behavior is found in form of line-like resonances, at low particle 
densities exactly reproducing the zero-temperature supersymmetric vector meson mass spectrum. 

Besides finite temperature effects it is especially interesting to observe the in-medium effects caused 
by finite particle density. We observe a broadening width of the resonance peaks in the spectrum as a 
function of increasing particle density. At the same time, a slight shift of the resonance position to higher 
energies occurs. This result contradicts the expectations from the effective QCD models investigated by 
Brown and Rho 176). However, experimental data from collision experiments at SPS do not support 
Brown-Rho scaling either. Instead more recent effective models, which are in good agreement with 
experimental data, also show a broadening of the p-meson resonance peaks accompanied by a small 
positive mass shift [771 IZSj- It would be interesting to investigate the mechanisms, that lead to the 
qualitative agreement of effective QCD models and the D3/D7 setup that we observed here. Other 
in-medium effects will be studied in section 14.31 There we determine the mass shift of mesons due to 
polarization of mesons by the presence of the gluonic background field in the plasma. 

Second, we switched on a nonzero isospin density, and equivalently an isospin chemical potential 
arises. The spectral functions in this case show a qualitatively similar behavior as those for baryonic 
potential. However, we additionally observe a splitting of the single resonance peak at vanishing isospin 
potential into three distinct resonances. This suggests that by explicitly breaking the flavor symmetry by 
a chemical potential, the isospin triplet states, vector mesons in our case, show a mass splitting similar 
to that observed for QCD [79]. It is an interesting task to explore the features of this isospin theory in 
greater detail in order to compare with available lattice data and effective QCD models ^82, i83l [84l [85l 
1551 1571 1551 IHni HOj- In most of these approaches, baryon and isospin chemical potential are considered 
at the same time, which suggests another promising extension of this work. Moreover, in the context of 
gravity duals, it will be interesting to compare our results for the isospin chemical potential to the work 
presented in ref. [5T] . 

Alternatively, instead of giving the gauge field time component a non- vanishing vacuum expectation 
value, one may also switch on i?-field components and combine the framework developed in [HH [Ml [S3] 
with the calculation of spectral functions for the dual gauge theory. 

Our spectra also show that for given quark mass and temperature, lower n meson excitations can be 
nearly stable in the plasma, while higher n excitations remain unstable. At vanishing baryon density, 
the formation of resonance peaks for higher excitations has also been observed in [94] . 

4 Diffusion in the holographic plasma 

Based on the observation that the many particle system observed at RHIC is well described by hydro- 
dynamics with very low viscosity, the quark-gluon plasma is widely regarded as an almost perfect liquid 
[551 ISni [SI]. In principle, the theory of hydrodynamics should be capable of modeling the collective 
dynamics of the plasma. This requires knowledge of initial conditions and the equations of state of 
the system. The hydrodynamic description then yields the dynamics of the system in terms of collective 
quantities such as currents, densities and entropy. For a comprehensive understanding of the quark-gluon 
plasma and related systems it would be pleasing to be able to derive the thermodynamic and kinetic 
properties of the system from first principles. In this section we present work which was conducted with 
the motivation to advance towards this goal. Adopting different points of view, we contemplate one 
particular attribute characteristic for fluids: diffusion. 

The diffusion coefficient is a transport coefficient which parametrizes the ability of a fluid to reach an 
equilibrium state by transport of some initially unevenly distributed quantity through currents. In the 
quark-gluon plasma, these currents are the color and flavor currents, which account for the transport of 
quarks and gluons through the plasma. A particle of high diffusivity (high mobility) looses only a small 
part of its energy while traversing a given distance the medium and will transport its associated charge 
much faster than a particle which looses much momentum by interaction with the medium. 

Close to thermodynamic equilibrium, transport coefficients such as that for diffusion can be derived 
from first principles by so called Kubo-formulae, which describe the coefficients in terms of correlation 
functions of the current which accounts for equilibration of the system. This approach has been used 
successfully in the past to derive transport coefficients and conductivities from holographic models [SH 
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In this section, once more the adjoint and fundamental matter described by the gauge multiplet of 
A/" = 4 and the Af — 2 hypermultiplet of thermal SYM theory in the limit of a large number of gauge 
degrees of freedom serves as a model for the quark-giuon plasma. 

The energy loss of heavy quarks and mesons in media has been a subject of intense experimental 
interest |109l IllOl lllll 11121 11131 11141 1115] . The suppression of charm and bottom quarks observed at 
RHIC motivated several groups to utilize the gauge-gravity duality [HI [131 HH H] to compute the drag of 
fundamental heavy quarks in A/" = 4 super Yang-Mills theory at strong coupling |103[ 11001 1102) . In this 
approach, the heavy quark is given by a classical string attached to the D7-brane probe. First studies 
of flavors in thermal AdS/CFT beyond the quenched approximation, i.e. with non-zero Nf/N, were 
performed in [1161 1117j . 

We pick up the previous efforts and generalize them by including the effects of finite density, respec- 
tively chemical potential, on diffusion properties. The first part of this section very shortly addresses 
consequences of finite baryon density on a certain holographic method to derive the diffusion coefficient 
of baryons in a holographic plasma. Subsequently, in section [4.21 we study the gauge/gravity dual of 
a finite temperature field theory at finite isospin chemical potential. The isospin chemical potential is 
obtained by giving a finite vacuum expectation value to the time component of the non-Abelian gauge 
field on the brane, as in the previous sections. In order to obtain analytical results, we restrict our 
attention to the limit of massless quarks. 

The consideration of an isospin chemical potential is an interesting field to study since it is still easier 
accessible by lattice calculations than setups at baryonic chemical potential. Hopefully holographic 
models one day allow for comparison to e. g. large N lattice calculations. Moreover, isospin diffusion has 
been measured in heavy ion reactions [1181 161] . 

Eventually in section H31 we adopt a different point of view on diffusion in the holographic quark-giuon 
plasma. From the technical point of view, we do not pursue the approach of Kubo to obtain the diffusion 
coefficient directly from current correlation functions. Instead, we make use of a stochastic Langevin 
model, which determines the momentum broadening of particles due to random kicks from interactions 
with the medium. 

Moreover, conceptually we extend the area of research on diffusion processes in the holographic 
QGP to mesons traversing the plasma. This effort bears two interesting aspects. One is the aim for a 
description of the kinetics of heavy mesons in the QGP, since observation show that heavy mesons like the 
J/V' survive the deconfinement transition. The other is the estimation of the effects of strong coupling 
on the plasma. The particular effective model we use does not rely on any weak or strong coupling limit 
for the interaction of the mesons with the medium. This allows for comparison of perturbative results 
for momentum broadening at weak coupling with holographic results for the strong coupling regime. In 
this way we present a method that may allow to estimate the effect of strong coupling on dynamic effects 
in the QGP. 

4.1 Baryon diffusion 

In this section, we calculate the baryon diffusion coefficient D and its dependence on the baryon density 
in the thermal holographic plasma. The coordinates and parameters we use are the same as in section [31 
they are discussed in detail in section [3^ 

In the context of holography, the idea is to describe the conserved current J of the gauge field theory 
in terms of the dual gauge field in the supergravity theory, as we did in section [31 This current in our case 
is the current J which transports baryon charge and therefore is non-zero when baryon diffusion occurs. 
The dual supergravity field was identified as the gauge field fluctuation A on the probe brane. Any 
gauge field configuration of A that satisfies the equations of motion also generates a conserved current 
in terms of the field strength F = dA, as dF — 0. As in electrodynamics one can identify the columns 
of the field strength tensor with vector currents. Pick's law of diffusion J — D'VJ'^ can be shown to 
be satisfied in the long distance limit for the on-shell field strength tensor of a supergravity gauge field. 
The constant D is then identified with the diffusion constant 24 . This constant describes how strong 
the currents J are which drive a hydrodynamic system into equilibrium, as a reaction on gradients in 
the charge distribution J*^. 

To solve for D, one therefore has to find the solutions of the gauge field fluctuations which are 
holographically dual to the relevant current. As soon as a solution is found and gauge/gravity duality 
is invoked to compute the current J, one may solve Pick's law for D. At vanishing particle density. 
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the gauge field solution can determined from the DBI action in terms of metric coefficients alone [24] . 
Because we use the same coordinate system as the authors of ref. [48j, we arrive at the same explicit 
form of the induced metric p.lSp on the probe D7-brane in our setup. We therefore reuse the result for 
the diffusion coefficient derived there. 

A very concise formula for the diffusion coefhcient D of R-charges was derived in reference to the 
membrane paradigm in ref. [24j . It was later directly translated to the diffusion of flavor currents in 
the D3/D7 setup in ref. [15]. The name "membrane paradigm" does not refer to D-branes but instead 
alludes to the fact that the analogies between black hole physics and thermodynamics very often can 
be expressed in terms of events taking place at the event horizon (or slightly outside the horizon, then 
referred to as the "stretched horizon"), which has no materialistic manifestation but still appears as a 
significant surface, a membrane, to an observer or in the relevant formulae. 

We extend previous efforts by introducing finite baryon density. This quantity explicitly enters the 
solution to the gauge field p.47p which in turn explicitly enters the equation of motion p.60p for the 
fluctuations. Therefore, one should expect a modified result for the diffusion coefficient including the 
explicit occurrence of the baryon density d. We rely on the fact that these terms vanish for d — >■ 0, 
restrict to the small density regime and stick to the expression found in ref. [15]. Nevertheless we still 
cover finite density effects in this way, since the probe brane embeddings x ^re different for vanishing and 
finite density, as we see from their equation of motion (|3.44p . This difference should therefore translate 
into a dependence of the diffusion constant on the parameter d, which is proportional to the baryon 
density by (j3.43p . We know that in the case of finite baryon density, black hole embeddings describe the 
entire parameter range of temperature and quark mass |60l 158 1 . 

Finally, the formula for the diffusion coefficient D found in ref. [48 is given by 



where the metric coefficients G^^ can be obtained from p.l8p with the square root of the absolute value 



of the metric determinant ^—G given by (I3.19P . Insertion of these coefhcients and Tq — TttB? yields 




The embeddings x a-i'^ determined as in section [3] by solving equation (|3.44p in terms of the parameters 
d for baryon density and initial value x(l) = Xo, which determines the quark mass normalized to 
temperature, c. f. figure 1331 

The results for D are shown in figure l4m where we compare to the result at vanishing baryon density 
found in in ref. |48| . There is a phase transition, at approximately m = 1.3 which we briefly address in 
section [O] 

The diffusion coefficient never vanishes in the medium with non-zero density. Both in the limit of 
T/mq — > and T/niq — )• oo, D converges to l/(27rT) for all densities, i.e. to the same value as for 
vanishing baryon density, as given for instance in j24) for R-charge diffusion. In the regime of moderate 
to low temperatures the diffusion constant develops a nonzero minimum. 

In order to give a physical explanation for this behavior, we focus on the case without baryons 
first. We see that the diffusion coefficient vanishes at the temperature of the fundamental deconfinement 
transition. This is simply due to the fact that at and below this temperature, all charge carriers are 
bound into mesons not carrying any baryon number. 

For non-zero baryon density however, there is a fixed number of charge carriers (free quarks) present 
at any temperature. This implies that the diffusion coefficient never vanishes. Switching on a very small 
baryon density, even below the temperature of the phase transition, where most of the quarks are bound 
into mesons, by demanding c? 7^ there will still be a finite amount of free quarks. By increasing the 
baryon density, we increase the amount of free quarks, which at some point outnumber the quarks bound 
in mesons. In the large density limit the diffusion coefficient approaches Dq = 1/ (27rT) for all values of 
T/rUq, because only a negligible fraction of the quarks is still bound in this limit. 

As a final comment, we point out that after the publication of these results in refs. [U [3] a more 
careful analysis of the calculation of the diffusion coefficient was performed in ref. |119| . Here the 
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Figure 4.1: Approximate baryon diffusion coefficient D as a function of the quark mass to temperature 
ratio m. The curve for d = reproduces the result from ref. [48]. The exact results derived in ref. [119) 
show deviations most notably in the limit of large to. 



diffusion coefficient was identified with the proportionality coefficient D in the dispersion relation for the 
hydrodynamic quasi normal mode uj (the so called diffusion pole or hydrodynamic pole, cf. section [3. ip 
of the gauge field fluctuations at finite spatial momentum k, given by 

uj = -iDk'^ + 0{k^) . (4.3) 

The inclusion of finite spatial momentum introduces several new aspects we circumvented in the limit 
of fc — > 0. In accordance with the results obtained from this analysis, the results for D in the low quark 
mass/high temperature regime agree well with our nai've approach [119j . 

4.2 Isospin diffusion 

We study the diffusion coefficient of particles charged under isospin chemical potential. The results 
where originally published as ref. [S] . Physically, the isospin chemical potential corresponds to the energy 
necessary to invert the isospin of a given particle. Within nuclear physics, such a chemical potential is of 
relevance e. g. for the description of neutron stars. In two-flavor QCD, effects of a flnite isospin chemical 
potential have been discussed for instance in refs. |120l [85l Il21j . 

In the following paragraphs we outline the procedure and comment on the restrictions we imposed. 
Recent work revealed that some of these can actually be considered as shortcomings. In particular the 
consideration of a constant background gauge fleld on the brane can be justifled only for very small 
chemical potentials. Due to the assumed smallness of the chemical potential we neglected second order 
terms in the background gauge fleld on the brane. Both limitations however can be cured, based on the 
insights published in refs. [60l[58]. A reviewed version of our results can be found in ref. 122_. 

As in section [31 we introduce an isospin chemical potential /i by defining the vacuum expectation 
value of the SU(2) gauge field on the coinciding probe D7-branes to be 

For simplicity, we work with a constant background field configuration. The above is a solution to 
the D7-brane equations of motion and is present even for the D7-brane embedding corresponding to 
massless quarks. However, We consider /i to be small, such that the Bose-Einstein instability observed 
in ref. |54j, which is of order 0(/z^), does not affect our discussion here. 

Even though, this constant gauge field given by (|4.4p is a solution to the equation of motion, it does 
not represent the thermodynamically preferred configuration [60j . Instead we should rather make use 
of the solution presented in p.47|) with constant embedding function % = 0, for which the authors of 
ref. [BD] present an analytical solution. Nevertheless, for simplicity we stick with (14.41) in this section, 
which may be justified in the case of very small densities, where the derivative of Ao{u) is very small, 
c. f. figure [3111 
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Again, we work in the D3/D7 setup. For simplicity, we consider only the D7 probe embedding for 
vanishing mass m = 0. This embedding is constant and terminates at the horizon. This simplification 
allows us to derive our results purely analytically. We establish the SU(2) non-Abelian action for a probe 
of two coincident D7-branes and obtain the equations of motion for fluctuations about the background 
(|4.4|) . These fluctuations are dual to the SU(2) flavor current J^"'. We find an ansatz for decoupling 
the equations of motion for the different Lorentz and flavor components, and solve them by adapting 
the method developed in refs. [HI [53]. This involves Fourier transforming to momentum space, and 
using a power expansion ansatz for the equations of motion. We discuss the approximation necessary 
for an analytical solution, which amounts to considering frequencies with uj < fi < T. With this 
approach we obtain the complete current-current correlator. The key point is that the constant chemical 
potential effectively replaces a time derivative in the action and in the equations of motion. In the 
Fourier transformed picture, this leads to a dependence of physical observables on the square root ^/uJ 
of the frequency. This non-linear behavior goes beyond linear response theory. We discuss the physical 
properties of the Green functions contributing to the current-current correlator. In particular, for small 
frequencies we find a frequency-dependent diffusion coefficient D{uj) oc ^^/uj/fl. Whereas frequency- 
dependent diffusion has — to our knowledge — not yet been discussed in the context of the quark-gluon 
plasma, it is well-known in the theory of quantum liquids and therefore may possibly also apply to 
the quark gluon plasma. For instance, for small frequencies the square-root behavior we find agrees 
qualitatively with the results of refs. [123lll24j for liquid para-hydrogen. Generally, frequency-dependent 
diffusion leads to a non-exponential decay of time-dependent fluctuations [125j . 

The approach used in this section is related but different from the procedure we implemented to 
obtain the spectral functions at finite isospin density in section [3] The non-analytical behavior we derive 
is due to the limits we consider. Coming from equation p.85p . the difference between a constant non- 
vanishing background gauge field and the varying one becomes clear. Here, the field is chosen to be small 
and constant in p, such that terms quadratic in the background gauge field Ag <C 1 can be neglected. 
This implies that the square (to =F n^)^ in (13.82p and p.83p is replaced by =F 2rDm, such that we obtain 

the indices (3 = ^^\f^~T^^^- instead of (|3.85p . If we additionally assume to <C Ag, then the 1 under 
the square root can be neglected. In this case the spectral function develops a non-analytic structure 
coming from the ^/uj factor in the index. 

This section is organized as follows. We start with a comment on frequency-dependent diffusion within 
hydrodynamics and the method we use to compute the diffusion constant holographically. Thereafter, 
we establish the D7 probe action in presence of the isospin chemical potential, derive the corresponding 
equations of motion and solve them. Finally, we obtain the associated Green functions in the hydro- 
dynamical approximation. From their pole structure we can read off the frequency-dependent diffusion 
coefficient. We comment on our results briefiy where it is appropriate in this section and leave a sum- 
marizing discussion for section [4.41 including an interpretation of our results. 

4.2.1 Diflfusion coeflRcients from Green functions 

Thermal Green functions have proven to be a useful tool not only to derive spectra, as above, but also 
for analyzing the structure of hydrodynamic theories and for calculating hydrodynamic quantities such 
as transport coefficients. In this section we once more use the gauge/gravity dual prescription of ref. [23] 
for calculating Green functions in Minkowski spacetime. These correlators can be thought of as being 
determined by their pole structure, in the way discussed in connection with the spectral functions, c. f. 
section 13.11 From these poles at complex frequencies we derive the diffusion coefficient D of isospin 
charge with charge density J° and conserved current four vector J = ( J°, J). 

Considering systems governed by hydrodynamics, such as the quark-gluon plasma, we are eager for 
solutions to the hydrodynamic equations of motion. Regarding diffusion, we pay special interest to the 
Green function for the diffusion equation 

9oJ°(i,x) = L'VV°(i,x), (4.5) 

with J*^ the density, given by the time component of a diffusive current four vector J, and D is the 
diffusion constant we are interested in. In Fourier space this equation reads 

J"(w, k) = L»kV°(a;, k). (4.6) 
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This determines the dispersion relation lo = —iD\<? of the mode with energy w. In the language of Green 
functions we will observe solutions of the form 



k) « , 2 , ^„ 3^ • (4-8) 



G(a.,k)(x: (4.7) 

ILO — UK 

Here, the dispersion relation determines the poles of the Green function. Finding the correct dependence 
of the poles of the correlator on w and k therefore allows to determine the diffusion constant D. 

The diffusive quantity we are interested in the isospin charge with density J*^ as part of the four vector 
J — (J", J). These currents are holographically dual to the gauge field on the brane. We therefore have 
to solve the equations of motion for the gauge fields to obtain the relevant field theory Green function G 
of isospin diffusion of J° by following the recipe for retarded correlation functions (3^ . We are interested 
in the hydrodynamic properties and therefore restrict to the hydrodynamic long wavelength/low energy 
limit such that we restrict to the lowest order in k. 

For the non-Abelian case with an isospin chemical potential, in section l4.2.4l we will obtain retarded 
Green functions of the form 

1 

L'(a;)l? + 0(k3) 

Retarded Green functions of this type have been discussed for instance in ref. |125j . Equation 
describes frequency-dependent diffusion with coefficient D{lo), such that (j4.6l) becomes 

zwJ°(w, k) = D{uj)\i^.f{uj, k) . (4.9) 

In our case, J*^ is the isospin density at a given point in the liquid. 

This is a non-linear behavior which goes beyond linear response theory. In particular, when Fourier- 
transforming back to position space, we have to use the convolution for the product D ■ Jq and obtain 

t 

doJ°{t,x)+V^ J dsJ°{s,x)D{t- s) ^0 (4.10) 

— oo 

for the retarded Green function. This implies together with the continuity equation OqJ'^ + V • J = 0, 
with J the three- vector current associated to J°, that 

J = -V{D* j"), (4.11) 

where * denotes the convolution. This replaces the linear response theory constitutive equation J = 
~-DVJ°. Note that for D{t - s) = DS{t - s) with D constant, ffTJU\i reduces to 

4.2.2 Holographic setup 

In this subsection we use the AdSblack hole coordinates given in appendix |21 equation (|A.4p , 

f 1 2 ^ p2 (4.12) 

< M < 1, e K 

with the metric dil| of the unit 5-sphere. This geometry is asymptotically AdSs x S'^with the boundary 
of the AdSpart at m = 0, the black hole horizon is located at w = 1. 

Into this ten-dimensional spacetime we embed A'^^ = 2 coinciding D7-branes, hosting the flavor gauge 
fields A. We choose the same embedding as in the previous calculations, which extends the D7-branes 
in all directions of AdSs space and wraps an on the S^. Here, we restrict ourselves to the most 
straightforward case, that is the trivial constant embedding of the branes through the origin along the 
AdSradial coordinate u. This corresponds to massless quarks in the dual field theory. On the brane, the 
metric simply reduces to 

ds^ = (~/(w) dt^ + dxj + dxj + dxj) + ^^^^^ du^ + R'dnl (4.13) 
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Due to the choice of our gauge field in the next subsection, the remaining three-sphere in this metric will 
not play a prominent role. We use labels fi,iy, . . . to denote any direction, ... to refer to Minkowski 
directions, u is used as a label for the radial coordinate, and a will be used to refer to the x^'^ directions. 



Introducing a non-Abelian chemical potential 

A gravity dual description of a chemical potential amounts to a non-dynamical time component of 
the gauge field in the action for the D7-brane probe embedded into the background given above. There 
are essentially two different ways to realize a non-vanishing contribution from a chemical potential to the 
field strength tensor F = 29[^yl^] -I- f°''"'A''^Al. The first is to consider a u-dependent baryon chemical 
potential, as we did in the preceding section. We work with a constant chemical potential of the form 

A^^A^T", (4.14) 

where we sum over indices which occur twice in a term and denote the gauge group generators by T". 
The brane configuration described above leads to an SU(iV/) gauge group with Nf = 2 on the brane, 
which corresponds to a global SU(iV/) in the dual field theory. For Nf = 2, the generators of the gauge 
group on the brane are given by = with Pauli matrices cr°. We will see that (j4.14p indeed 
produces non-trivial new contributions to the action. 

Using the standard background field method of quantum field theory, we consider the chemical 
potential as a fixed background and study gauge field fluctuations around it. We single out a particular 
direction in flavor space by taking Aq — fi as the only non- vanishing component of the background field. 
From now on we use the symbol A"^ to refer to gauge field fiuctuations around the fixed background, 

Al ^ fiS,f,6^^ + Al. (4.15) 

We gauge the component along the radial coordinate to Au — and assume that yl^ = for fi — 5,6, 7. 
Due to the symmetries of the background, we effectively examine gauge field fiuctuations depending 
on the five-dimensional subspace on the brane spanned by the coordinates x*''^'^''^ and by the radial 
AdScoordinate u. The magnitude of all components of A and the background chemical potential fi are 
considered to be small. This allows us to simplify certain expressions by dropping terms of higher order 
in A and in the chemical potential fj,. 



4.2.3 Equations of motion and their solutions 

The action describing the dynamics of the fiavor gauge fields in the D3/D7 setup is the Dirac-Born-Infeld 
action. Since we work with vanishing gauge field components in all of the directions perpendicular to 
the D3-branes, there are no contributions from the Chern-Simons action. As mentioned, we consider the 
constant D7 probe embedding corresponding to vanishing quark mass, m — 0. The metric on the brane 
is then given by (|4.13p . Since we are interested in two-point correlators only, it is sufficient to consider 
the DBI action to second order in a' , 

u=l 

5d7 = -TjTr {2Tr^a')^R^ J dud^x yf^G^"" G"^ F^^ F^p , (4.16) 

ti=0 

where we use the D7-brane tension Ty as in (j2.14l) . performed the integration over the 5, 6, 7-directions, 
which are the directions along the , and the factor Tr arising from the trace over the representation 
matrices T°, 

Ti{T°-T^)=Tj^5°-\ (4.17) 

hi our case we have Tr = 1/2. 

Evaluating the DBI action given in (|4.16l) with the substitution rule (|4.15l) . we arrive at 

u=l 

5-07= -TjT^X2TT'^a'fR^ [ dud'^x V^G^"'G''^ 

uio (4-18) 

X Ud^^A:^ - sdoJopr'^d^.A;^ Al ^i ^ 
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where we use the short-hand notation /^"^ = fiS^'^ and neglect terms of higher than hnear order in /i, and 
higher than quadratic order in A since both are small in our approach. 

Up to the sum over flavor indices a, the first term in the bracket in (j4.f 8p is reminiscent of the 
Abelian super-Maxwell action in five dimensions, considered already for the R-charge current correlators 
in ref. |22j . The new second term in our action arises from the non- Abelian nature of the gauge group, 
giving terms proportional to the gauge group's structure constants f"^'"^ in the field strength tensor 

We proceed by calculating the retarded Green functions for the action (|4.18l) , following the prescrip- 
tion of ref. [23 as outlined in section [2.2.1l According to this prescription, as a first step we consider the 
equations of motion obtained from the action (j4.18L which are given by 

r / M (4.19) 

+ r'"' V^G^^G"/^ {dpAl - 2^0^^) + S°''d^ f V^G^G^'^A^Ai'') ■ 

It is useful to work in momentum space from now on. We therefore expand the bulk gauge fields in 
Fourier modes in the directions, 

M^, ^) = e~^-■^o+^l^-A,{u, k). (4.20) 

As we work in the gauge where A^ = 0, we only have to take care of the components Ai with i = 0,1, 2, 3. 

For the sake of simplicity, we choose the momentum of the fluctuations to be along the direction, 
so their momentum four-vector is fc = {uj,0,0,q). With this choice we have specified to gauge fields 
which only depend on the radial coordinate u, the time coordinate x'^ and the spatial x^ direction. 

Equations for A"^- and Aj'Components 

Choosing the free Lorentz index in the equations of motion (|4.19p to be = a = 1, 2 gives two identical 
differential equations for Ai and A2, 



/ " ^ Up ^" + ^\/2^ 2^r^ 



= AS" + + tF-^o + ^^^r'TTT^^A^c., (4.21) 



where we indicated the derivative with respect to u with a prime and have introduced the dimensionless 
quantities 

"'^2^' ''=2^' "^=2^- (^-^^^ 

We now make use of the structure constants of SU(2), which are /"'"^ = e"'"^, where £°'"^ is the totally 
antisymmetric epsilon symbol with e^'^'^ — 1. Writing out (|4.2ip for the three different choices of a = 1, 2, 3 
results in 

= Ai" + ^Ai' -f ^^^Ai - 2.^A^ , (4.23) 
= Ar + fAr + ^^i^Ai + 2.^Ai, (4.24) 

= A3 " + ^Ar + "^^^Al ■ (4.25) 

The first two of these equations for the gauge field directions transverse to the background field are 
coupled, the third longitudinal one is the same equation that was solved in the Abelian Super-Maxwell 
case [22]. Note that these equations are influenced by one of the addressed oversimplifications in a way 
that will turn out to be crucial at the end. We neglect terms quadratic in m. If we would allow for these 
quadratic terms we could complete the square and avoid the non-analytical behavior we observe later 
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Equations for Aq- and ^g-components 

The remaining choices for the free Lorentz index = 0, 3, u in (j4.19p resuh in three equations which are 
also not independent. The choices v = and u = u give 

^-^^ uf^^ 2^T^^' ^^-^^^ 

= lr< + q/Ar+*r''=^<. (4.27) 

Solving (I4.27P for A^' , differentiating it once with respect to u and using (|4.26p results in equation (j4.19p 
for v = 3, 

= Af + ^A%' + ^A% + ^A^r, + i-^L/'^fcc + 2i-^r'"'^ A% . (4.28) 

3 / ^ up ^ up ° up ■' 2ttT ° up-' 2ttT ^ ^ ' 

We will make use of the equations (j4.26p and (j4.27p which look more concise. These equations of motion 
for Aq and A'^ are coupled in Lorentz and flavor indices. To decouple them with respect to the Lorentz 
structure, we solve (|4.27p for A'^' and insert the result into the differentiated version of (|4.26p . This gives 



= Ao +2z— / —A,. (4.29) 



The equations for a = 1,2 are still coupled with respect to their gauge structure. The case a = 3 was 
solved in ref. (22]. We will solve (|4.29p for A^ and can obtain A'^ from (|4.27p . Note that it is sufficient for 
our purpose to obtain solutions for the derivatives of the fields. These contribute to equations (|2.58p to 
(|2.6ip that give the retarded thermal Green functions, while the functions A = f{u, k)A^'^y(k) themselves 
simply contribute a factor of f{u, —k) which merely gives a factor of unity at the boundary. 

Solutions 

Generally, we follow the methods developed in ref. |22| , and our differential equations are very similar to 
the ones considered there. Additionally, we need to respect the flavor structure of the gauge fields, only 
the equations for flavor index a = 3 resemble those analyzed in ref. 22]. Those for a = 1, 2 involve extra 
terms, which couple the equations. Coupling occurs not only via their Lorentz indices, but also with 
respect to the flavor indices. We already decoupled the Lorentz structure in the previous section. As 
typical for an explicitly broken SU(2), the equations of motion which involve different gauge components 
will decouple if we transform to the variables 

X, = A\+iAj, 

; (4.30) 

Here, the A], Af are the generally complex gauge field components in momentum space. Note that up 
to SU(2) transformations, the combinations p.8ip are the only ones which decouple the equations of 
motion for a = 1, 2. These combinations are reminiscent of the non-Abelian SU(2) gauge field in position 
space. 

The equations of motion for the flavor index a = 3 were solved in ref. [22]. To solve the equations 
of motion for the fields Af with a = 1,2, we rewrite them in terms of Xi and Xi. Applying the 
transformation (|3.8ip to the equations of motion (|4.23l) and (I4.24p and the a = 1,2 versions of (I4.29P 
and ((i??7)) leads to 



/ " up 

uf uf^ 
= (trTm)X^ + q/X^, (4.34) 



= + + i^Lf^x^ , (4.32) 

^ X- + ^ X^' + " (4.33) 
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where again a — 1,2, and the upper signs correspond to X and the lower ones to X. 

As in the section on spectral functions, we observe that some coefficients of these functions are 
divergent at the horizon u = 1. We hark back to the ansatz 

X, = il- uf F(u), X, = (1 - uf Fiu), (4.35) 

with regular functions F(u) and F(u). To cancel the singular behavior of the coefficients, we have to 
find the adequate /3 and /?, the so-called indices, given by equations known as the indicial equations for 
/? and /?. We eventually get for all Xi and Xi 



^ , jtD / 2m ~ itn / 2m „„s 

Note that these exponents differ from those of the Abelian Super- Maxwell theory [22] by a dependence 
on ^/m in the limit of small frequencies (tn < m). In the limit of vanishing chemical potential m — ?> 0, the 
indices given in ref. |22) are reproduced from p.62p . Again, however, if we include the quadratic order 
in m it is possible to complete the square to get /3 = ^i/{2w =F tn) |122) . 

In order to solve (j4.32L (j4.33p and (|4.34p analytically, we introduce a series expansion ansatz for 
the function F in the momentum variables to and q. In section [3.3.11 we solved the resulting equations 
up to first order in the radial coordinate in order to get initial conditions for the subsequent numerical 
integration. Here we are not interested in the dependence along u, but in the dependence on the lowest 
order in in and q in order to extract the dispersion relation that determines the poles in the according 
correlators. In fact, the physical motivation behind this expansion is that we aim for thermodynamical 
quantities which are known from statistical mechanics in the hydrodynamic limit of small four-momentum 
k. So the standard choice would be 

F{u) =Fo + mFi + q^Gi + . . . . (4.37) 

On the other hand, we realize that our indices will appear linearly (and quadratically) in the differential 
equations' coefficients after inserting (I4.35P into (I4.32p . (|4.33p and (|4.34p . The square root in /3 and (3 
mixes different orders of in. In order to sort coefficients in our series ansatz, we assume tn < m and keep 
only the leading tn contributions to /3 and (3, such that 



This introduces an additional order 0(tn^/^), which we include in our ansatz (|4.37p . giving 

F(u) = Fo + tni/2^1/2 + tnFi + q^Gi + . . . , (4.39) 

and analogously for the tilded quantities. If we had not included 0(tn^/^) the resulting system would be 
overdetermined. On the other hand this procedure of including non-analytical square root terms would 
be superseded by including the second order terms in m from the beginning. The results we obtain by 
using the approximations (j4.38l) and (|4.39p have been checked against the numerical solution for exact 
/? with exact F{u). These approximations are useful for fluctuations with q, tn < 1. 
Note that by dropping the 1 in p.62p we also drop the Abelian limit. 

Consider the indices (j4.36p for positive frequency first. In order to meet the incoming wave boundary 
condition, we restrict the solution (3 to the negative sign only. For the approximate /? in (j4.38p we 
therefore choose the lower (negative) sign. This exponent describes a mode that travels into the horizon 
of the black hole. In case of /?, we demand the mode to decay towards the horizon, choosing the lower 
(positive) sign in (|4.38p consistently. Note that for negative frequencies a; < the indices (3 and /3 
exchange their roles. 

Using (|4.38|) in (|4.35p and inserting the ansatz into the equations of motion, wc find equations for 
each order in q^ and tn separately. After solving the equations of motion for the coefficient functions Fq, 
Pi/2i Fi and Gi, we eventually can assemble the solutions to the equations of motion for X as defined 
in Km . 

X(u) = {\~uf F{u) 

' ^ ' (4.40) 

= (1 - uf (Fo + %/mFi/2 + tnFi + q^Gi + ...). 
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and a corresponding formula for X(u) from the ansatz (|4.35p . 

Illustrating the method, we now write down the equations of motion order by order for the function 

Xa- To do so, we use (I4.40p with (|4.38p in (|4.32p with the upper sign for Xa- Then we examine the 
result order by order in ro and q^, 

O(const) : = F^' + ^F^, , (4.41) 

0{V^: Q^F';,, + ljF[,,-^F^-^^^F,, (4.42) 

OM ■ 0-F,+jF,- Y^F,/, - ^j^jF,,^ - m F, , (4.43) 



0(q2): Q^G'i+^-^G\-—Fo. (4.44) 

At this point we observe that the differential equations we have to solve for each order are shifted 
with respect to the solutions found in ref. [22]. The contributions of order tn" in ref. [32] now show up 
in order tn"/^. Their solutions will exhibit factors of order /x"/^. Again, we emphasize that inclusion of 
O(m^) terms would just result in a shift of tr M- tr ± m. 

Solving the system (|4.4ip to (I4.44p of coupled differential equations is straightforward in the way 
that they can be reduced to several uncoupled first order ordinary differential equations in the following 
way. Note that there obviously is a constant solution Fq — C for the first equation. Inserting it into 
(I4.42P and (j4.44D leaves us with ordinary differential equations for F[^^ and G'l respectively. Using the 
solutions of Fq and in (j4.43D gives one more such equation for F{. 

To fix the boundary values of the solutions just mentioned, we demand the value of F{u — 1) to 
be given by the constant Fq and therefore choose the other component functions' solutions such that 
limu-j.ii^i/2 = 0, and the same for Fi and Gi. The remaining integration constant C is determined 
by taking the boundary limit u — )• of the explicit solution (j4.40p . making use of the second boundary 
condition 

lim X{u) = X^'^y, (4.45) 

see appendix IB. II Eventually, we end up with all the ingredients needed to construct the gauge field's 
fluctuations X{u) as in (I4.40p . 

We solve the equations (|4.32p with lower sign for Xa and (|4.33p for Xq and its tilded partner in 
exactly the same way as just outlined, only some coefficients of these differential equations differ. The 
solution for Xg is then obtained from (j4.34p . 

All solutions are given explicitly in Appendix IB. II together with all other information needed to 
construct the functions Xa, Xa, Xq, Xq, X^ and X'^. 

4.2.4 Current correlators 

In this section we obtain the momentum space correlation functions for the isospin currents by means 
of the holographically dual gauge field component combinations X and X defined in equation (j4.30p . 
Recall that the imaginary part of the retarded correlators essentially gives the thermal spectral functions 
(c. f. section The following discussion of the correlators' properties can therefore be related to the 
discussion of the corresponding spectral functions. 

First note that the on-shell action gets contributions from the non-Abelian structure, 



Su7 = - TtTr (27r2a')'i?' 



X 2 



d^g 

(M 



11=1 




(4.46) 



where j,k — 0,1,2,3, and the index u denotes the radial AdS-direction. Up to the sum over flavor 
indices, the first term in the bracket is similar to the Abelian Super-Maxwell action of ref. [52]. The 
second term is a new contribution depending on the isospin chemical potential. It is a contact term which 
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we will neglect. The correlation functions however get a structure that is different from the Abelian case. 
This is due to the appearance of the chemical potential in the equations of motion and their solutions. 
Writing function of X and X results in 

d^g 



(x/Xfe + X'^Xk) + AfA 



Aq^i I dwV-GG°°G33 2X[oX3] 



(4.47) 

In order to find the current correlators, we apply the method outlined in section [2.2.11 to (|4.47l) . 
with the solutions for the fields given in appendix IB. II As an example, we derive the correlators Gqq = 

(^Joiq)Joi~'t)'^ and Gqq = (^Jo{q) J(){—q)'^ of the flavor current time components Jq and Jq, coupling 

to the bulk fields Xq and Xq, respectively. Correlation functions of all other components are derived 
analogously. 



Green functions: Calculation 

First, we extract the prefactor of {OAq)'^ from the action (j4.16p and call it B{u), 

B{u) = -TjTn (27r2a')^i?^V^G™ G™. 
We need this factor below to calculate the Green function, 



(4.48) 



= lim B{u)f{u, ~k)duf{u, k). 



(4.49) 



The second step, finding the solutions to the mode equations of motion, has already been performed 
in section 14.2.31 In the example at hand we need the solutions Xq and Xq . From (j4.40p and from 
appendix IB] we obtain 



Xq' 



(l-^.)V- 



\/2mxv + tnm In 2 + 

2^2 



(4.50) 



2 \u + l 



/ 2 

12 



31n2 2 + 31n2(l 
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l + u 



12Li2(l -u) + 12Li2(-M) - 12Li2 



l-u 



q^ln 



u + 1 
2u 



irqX 



bdy 



i V2mrD + rom In 2 — q^ 

2 



(4.51) 



m , ( 2u 
in 

2 \u + l 



tr — TT^ + 3 In^ 2 + 3 In^ (1 + u) + 6 In 2 In 



l + u 



12Li2(l -u) + 12Li2(-u) - 12Li2 



l-u 



q^ln 



u+l 
2u 



Note that we need the derivatives to apply (12.59^ and (j2.60p . 

Now we perform the third step and insert (|4.48p . (j4.50l) and (|4.5ip into (|4.49l) . Our solutions Xq and 
Xq replace the solution /(u, k) and /(u, —k) in (j2.60l) . The resulting expression is evaluated at Ub ~ 0, 
which comes from the lower limit of the u- integral in the on-shell action (j4.47p . At small u = e ^ 1, 
(|i3Dl) and ((i3T|) give 



lim Xq' 

u->0 



lim Xq' 

u-)-0 



q2XQ^'^^ + trqX3^'^^ 
V2mrD + rom In 2 + q^ e 



lim (q^Xo^'^y + XX) 



2 ybdy 

q ^0 



jV2mrD + rom In 2 — q^ « 



lim {fX^^ 



q^a'") Ine, 



(4.52) 
(4.53) 
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In the next to leading order of (I4.52p and (|4.53p there appear singularities, just like in the Abelian 
Super-Maxwell calculation ^22, equation (5.15)]. However, in the hydrodynamic limit, we consider only 
the finite leading order. 



2 



Green functions: Results 

Putting everything together, for the two Green functions for the field components Xq, Xq given in (|4.30p 

by ^ 

Xo=Al + iAl, Xo = Al~iA, 

we obtain 

NcT 2ttT q 

Stt iy/2mtv - q2 + rom In 2 



Goo = ^ T^. — - (4-54) 



NcT 2ttT q2 

Svr -^/2mtv - - rom In 2 



goo=^^ ^ (4-55) 



These are the Green functions for the time components in Minkowski space, perpendicular to the chemical 
potential in flavor space. All Green functions are obtained considering hydrodynamic approximations in 
0(tD-^/^, ttJ, q^), neglecting mixed and higher orders 0(03"^/^, ro^/^q^, q"*). 

The prefactor in (|4.54p . (I4.55P is obtained using Tj as in (|2.14l) . Tr from (|4.17p . and carefully inserting 
all metric factors, together with the standard AdS/CFT relation R"^ = AnggNa' . As in other settings 
with flavor [108] . we concordantly get an overall factor of N, and not N'^, for all correlators. Contrary 
to those approaches, we do not get a factor of Nf when summing over the different flavors. This is due 
to the fact that in our setup, the individual flavors yield distinct contributions. Most striking is the 
non-trivial dependence on the (dimensionless) chemical potential m in both correlators. Note also the 
distinct structures in the denominators. The first one, (j4.54l) . has an explicit relative factor of i between 
the terms in the denominator. In the second correlator, (j4.55p . there is no explicit factor of i. The 
correlator (|4.54[) has a complex pole structure for w > 0, but is entirely real for uj < 0. On the other 
hand, (j4.55p is real for a; > but develops a diffusion structure for w < 0. So the correlators Ggg and 
Gqq essentially exchange their roles as w changes sign. We find a similar behavior for all correlators G^j 

and with j, 1 — 0,1, 2, 3. Once more, this behavior is a consequence of the insertion of 0(txi^/^) and 
neglecting of terms of order O(tn^) in the hydrodynamic expansion (|4.39p . 

We assume tn to be small enough in order to neglect the denominator term of order O(rDm) ^ 
0(-\/lrm., q^) . Moreover, using the definitions of lT),q and m from (j4.22l) we may write (j4.54l) and (|4.55p 



NT 

8TTy/2jI uj + q^D{uj) 
NT 

%'K\/2[i iiu — q'^D{uj) 



Goo^-^^T^-^TTKTT' (4-56) 
Goo= ^ ...Tr.r.. ^ (4-57) 



where the frequency-dependent diffusion coefficient D{uj) is given by 

, X / u} 1 

We observe that this coefficient also depends on the inverse square root of the chemical potential /i. Its 
physical interpretation is discussed below in section [4.2.51 
In the same way we derive the other correlation functions 



NT uj^/^co-fi) NT uj^/^uj + f,) 



NT y^qjuj-fi) NT ^ q[u + 



_ NT co^/'q NT u^l^q 
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with the short-hand notation 



Q{uJ,q)=iu^q'^D{u), Q{u,q) =uj + q^D{u). (4.62) 

Note that most of these functions are proportional to powers of q and therefore vanish in the hmit of 
vanishing spatial momentum q ^ 0. Only the 33-combinations from (j4.59p survive this limit. In contrast 
to the Abelian Super- Maxwell correlators from ref. [22] given in appendix IB.2[ it stands out that our 
results (I4.56p . (|4.57p and (|4.59p and (|4.6ip have a new zero at to — ±/i. Nevertheless, bear in mind 
that we took the limit w < /i in order to obtain our solutions. Therefore the apparent zeros at ±/i lie 
outside of the range considered. Compared to the Abelian case there is an additional factor of y/uj. The 
dependence on temperature remains linear. 

In the remaining X-correlators we do not find any pole structure to order y/uj, subtracting an 0(^q^^ 
contribution as in ref. [21], 

V2N,T 

GiT =^22 = — ^ VJ^, (4-63) 

an 



iV2NcT 
"8^ 



Gn =^22 = -^^^VJ^ ■ (4-64) 



We can see that the GaS (with a ~ 1,2) are purely imaginary for negative uj and real for positive uj. 
The opposite is true for Gsq, as is obvious from the relative factor of i. 

The correlators of components, pointing along the isospin potential in flavor space (a — 3), are found 
to be 

N,T q^ 

iVcT ujq ^ ' 

NrTiuj 



G A3 — G A3 A3 — 



2 (4.66) 



47r iuj — Daq'^ ' 

with the diffusion constant Dg = l/(27rT). Note that these correlators have the same structure but differ 
by a factor 4/7V from those found in the Abelian super-Maxwell case [52] (see also (jB.36|) and (|BT38|). 
In particular the correlators in equation (|4.65p do not depend on the chemical potential. 

To analyze the novel structures appearing in the other correlators, we explore their real and imaginary 
parts as well as the interrelations among them, 

ReG„,(.>0)= R.G5„(„<0) = -^-^^=^l^, («7) 



G.,(. < 0) ^ - 1„ G,„(. > 0) = ^J^-^^^^, . (4.69) 

ImGoo(a; > 0) = ImG^^iuj < 0) = 0. (4.70) 

Now we see why, as discussed below (|4.57p . Gqq and Gqq exchange their roles when crossing the origin 
at a; = 0. This is due to the fact that the real parts of all G^j and Gjj are mirror images of each other by 
reflection about the vertical axis at a; = 0. In contrast, the imaginary parts are inverted into each other 
at the origin. The real part shows a deformed resonance behavior. The imaginary part has a deformed 
interference shape with vanishing value for negative frequencies. All curves are continuous and finite at 
UJ = 0. However, due to the square root dependence, they are not differentiable at the origin. Parts 
of the correlator which are real for positive uj are shifted into the imaginary part by the change of sign 
when crossing a; = 0, and vice versa. 
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To obtain physically meaningful correlators, we follow a procedure which generalizes the Abelian 
approach of ref. |126) . In the Abelian case, gauge- invariant components of the field strength tensor, 
such as Ea — ujAa, are considered as physical variables. This procedure cannot be transferred directly 
to the non- Abelian case. Instead, we consider the non-local part of the gauge invariant TrF^ which 
contributes to the on-shell action (|4.46|) . In this action, the contribution involving the non- Abelian 
structure constant — as well as /i — is a local contact term. The non-local contribution however generates 
the Green function combination 

G^iAi + Ga^a^ + Ga^a^ ■ (4-71) 

We take this sum as our physical Green function. This choice is supported further by the fact that it 
may be written in terms of the linear combinations p.8ip which decouple the equations of motion. For 
example, for the time component, written in the variables Xq, Xq given by p.8ip . the combination (I4.7ip 
reads (compare to (|4.47p ) 

<^oo + Gqo + Ga^aI ■ (4.72) 

The contribution from Ga^a^ is of order O^/i*^), while the combination for the first two flavor directions, 
Gqq + Gqq, is of order 0{fi). 

We proceed by discussing the physical behavior of the Green function combinations introduced above. 
Their frequency dependence is of the same form as in the Abelian correlator obtained in ref. [52] , as can be 
seen from (jB.36l) . Since we are interested in effects of order 0(/i), we drop the third flavor direction a = 3 
from the sum (|4.72p in the following. It is reassuring to observe that the flavor directions a = 1, 2, which 
are orthogonal to the chemical potential, combine to give a correlator spectrum qualitatively similar 
to the one found in ref. [32] for the Abelian Super-Maxwell action. However, we discover intriguing 
new effects such as the highly increased steepness of the curves near the origin due to the square root 
dependence and a kink at the origin — which have to be seen with skepticism because they vanish upon 
reinstating terms of order O(m^). 

We observe a narrowing of the inverse resonance peak compared to the form found for the Abelian 
Super-Maxwell action (and also compared to the form of our G^a^a . At the origin, the real and imaginary 
part are finite and continuous, but they are not continuously differentiable. However, the imaginary part 
of Ga^aI has finite slope at the origin. The real part though has vanishing derivative at a; = 0. 

The correlators G33, G33, Gqj and Gg3 have the same interrelations between their respective real and 
imaginary parts as Gpg and Gqq. Nevertheless, their dependence on the frequency and momentum is 
different, as can be seen from (|4.59|) to (j4.6ip . A list of the 33-direction Green functions split into real 
and imaginary parts can be found in appendix IB. 31 



4.2.5 Isospin diflFusion coefficient 

The attenuated poles in hydrodynamic correlation functions have specific meanings (for exemplary dis- 
cussions of this in the context of AdS/CFT see e.g. refs. [5511127] ). In our case we observe an attenuated 
pole in the sum Gqq + Gqq at w = 0. The pole lies at Rew = 0. This structure appears in hydrodynamics 
as the signature of a diffusion pole located at purely imaginary u>. Its location on the imaginary w-axis 
is given by the zeros of the denominators of our correlators as (neglecting 0(^u>, g^)) 



Squaring both sides of (|4.73p we see that this effect is of order O(g^). On the other hand, looking for 
poles in the correlator involving the third flavor direction G^s^a, we obtain dominant contributions of 
order O(g^) and 0(At°) (neglecting 0(w^,g^)) 

(4.74) 



This diffusion pole is reminiscent of the result of the Abelian result of ref. 122] given in appendix IB. 21 
As discussed in section 14.2.41 we consider gauge invariant combinations Gpg -|- Gqq -|- G^a^a. In order 
to inspect the non- Abelian effects of order 0(/x) showing up in the first two correlators in this sum, we 
again drop the third flavor direction which is of order ©(/i*'). 
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Motivated by the diffusion pole behavior of our correlators in flavor-directions a = 1,2 corresponding 
to the combinations X, X (see l|4.73p ). we wish to regain the structure of the diffusion equation given in 
(|4.9p . which in our coordinates (fc = (w, 0, 0, g)) reads 

iujJo = D{uj)q^Jn. (4.75) 

Our goal is to rewrite (I4.73P such that a term of 0{uj) and one term of order 0[q^) appears. Furthermore 
there should be a relative factor of —i between these two terms. The obvious manipulation to meet these 
requirements is to multiply (j4.73p by ^/uJ in order to get 



Comparing the gravity result (j4.76p with the hydrodynamic equation (|4.75p , we obtain the frequency- 
dependent diffusion coefficient 

Our argument is thus summarized as follows: Given the isospin chemical potential as in (j4.4p . (j4.15p . Jo 
from (|4.9p is the isospin charge density in (j4.75l) . According to (j4.75l) , the coefficient (I4.77P describes the 
diffusive response of the quark-gluon plasma to a gradient in the isospin charge distribution. For this 
reason we interpret D{ui) as the isospin diffusion coefficient. 

Near the pole, the strongly coupled plasma behaves analogously to a diffractive medium with anoma- 
lous dispersion in optics. In the presence of the isospin chemical potential, the propagation of non-Abelian 
gauge fields in the black hole background depends on the square root of the frequency. In the dual gauge 
theory, this corresponds to a non-exponential decay of isospin fluctuations with time. 

The square root dependence of our diffusion coefficient is valid for small frequencies. As long as 
ijj/T < 1/4, the square root is larger than its argument and at ui/T = 1/4, the difference to a linear 
dependence on frequency is maximal. Therefore in the regime of small frequencies uj/T < 1/4, which is 
accessible to our approximation, diffusion of modes close to 1/4 is enhanced compared to modes with 
frequencies close to zero. 



4.3 Meson diffusion at strong and weak coupling 

In this section, we consider heavy mesons moving slowly through high temperature non-Abelian plasmas. 
In the context of transport properties of the holographic quark-gluon plasma we are mainly interested in 
the diffusion behavior of mesons. The central quantity we discuss here will however not be the diffusion 
coefficient, but its inverse, the momentum broadening coefficient k, which determines the square of the 
momentum transfer per unit time, as we will see below. 

The motivation for considering meson diffusion is twofold. First, future experiments at RHIC promise 
to measure the elliptic flow of J/Tp mesons, and it is important to support this experimental program with 
theoretical work. To this end, various groups have studied the thermal properties of heavy mesons within 
the context of the AdS/CFT correspondence |128Lll29lll30LllOJ] . However, in spite of this progress, the 
transport properties of these mesonic excitations are not well understood. Although the kinetics derived 
in this work are not directly applicable to the heavy ion experiments, we believe that the results do hold 
some important information for phenomenology. 

The second motivation for this work is theoretical. After the quark drag was computed using the 
correspondence, it was realized that the drag of quarkonium is zero in a large Nc limit [1281 11291 1130j . 
Since within a thermal environment the drag and diffusion of these mesonic states is certainly not zero, 
it remained as a theoretical challenge to compute the kinetics of these states using the AdS/CFT setup. 

As a central result, we compare the diffusion of mesons at weak and strong coupling. Using a simple 
dipole effective Lagrangian which does not rely on the value of the coupling, we calculate the in-medium 
mass shift and the drag coefficient of the meson in A/" = 4 Super Yang Mills theory. At weak coupling 
we use perturbative methods, at strong coupling holographic models are employed. In the large N limit 
the mass shift is finite while the drag is suppressed by We reach the conclusion that relative to 

weak coupling expectations the effect of strong coupling is to reduce the momentum diffusion rate and 
thereby increase the relaxation time, which measures the time until the mesons in the plasma equilibrate 
their momentum spectrum to that of the thermal medium. 
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We also briefly pick up the discussion of in-medium efi^ects on meson spectra, subject of section [31 
There, the width of mesons in hot dense media was holographicahy determined by extending the analysis 
of meson melting to finite baryon density. In general the meson lifetime determined in this way is 
suppressed by the density of heavy quarks. However, we do not address the effects of finite density 
in this section. We are concerned with the thermal effects which capture the rescattering between the 
meson and the surrounding M = A medium. 

We focus on heavy mesons where the binding energy is much greater than the temperature. In this 
tight binding regime, mesons survive well above the critical temperature Tc for deconfinement and the 
meson width is sufficiently narrow to speak sensibly about drag and momentum diffusion. This behavior 
was observed for holographic models in section |3l 

For real charmonium, the binding energy can be estimated from the mass splitting AAf//^j^g « 
589MeV between the 2s and Is states, and for bottomonium from the 3s and Is states with AM^^_ig « 
895MeV respectively |131] . Therefore it is not really clear that real quarkonia above Tc « 170-190MeV 
[1321 1133j can be modeled as a simple dipole which lives long enough to be considered a quasi-particle. 
Indeed weak coupling hot QCD calculations of the spectral function show that over the temperature 
range gl^^M-g^^^M , the concept of a meson quasi-particle slowly transforms from being well defined 
to being increasingly vague [1341 11351 11361 11371 1138) . There is lattice evidence based on the maximal 
entropy method (which is not without uncertainty) that J/^ and T survive to 1.6 Tc and approximately 
3Tc respectively [l39l l49l [50l [T40l [TiTl [142] . However, model potential calculations which fit all the 
Euclidean lattice correlators indicate that the J/i/' and T survive only up to at most 1.2 Tc and 2.0 Tc 
respectively [1431 1144] . Clearly, the word "survive" in this context is qualitative and means that there 
is a discernible peak in the spectral function. Given these facts, our assessment is that the dipole 
approximation might be reasonable for Ti^ but poor for charmonium states and other bottomonium 
states. 

An overview of this section is as follows. First, in section [4.3.11 we review the computation of drag 
and diffusion of heavy quark-antiquark bound states within the setup of perturbative QCD. This will 
outline a two step procedure to determine the drag coefficient at strong coupling. 

The first step is to determine the in-medium mass shift 5M (it is finite at large N in the quantum 
field theory) , which determines the polarizabilities of the meson. As expected from the dipole effective 
theory, the mass shift scales as T'*/A'^, with T the temperature and Kb the inverse size of the meson. 
In the perturbative calculation. As is the inverse Bohr radius, while in the AdS/CFT computation the 
meson mass plays this role. In the A/" = 4 field theory the dipole effective Lagrangian couples the heavy 
meson to the stress tensor and the square of the field strength, which we denote by the operator Opi. 
In AdS/CFT we obtain the mass shifts from the linear response of the meson mass by switching on 
the dual operators. This amounts to consider a black hole background or a non-trivial dilaton fiow, 
respectively. For the dilaton flow we consider the D3-D(— 1) gravity background of Liu and Tseytlin 
[145| . This background and the AdS-Schwarzschild background allow for an analytic calculation of the 
meson polarizabilities. 

The second step is to compute the force-force correlator on the meson using the previously computed 
polarizabilities. This determines the drag coefficient rjo and the momentum broadening k as reviewed in 
section [4.3.31 This step requires the calculation of two-point functions involving gradients of the stress 
tensor and the field strength squared. Within gauge/gravity duality, these are obtained by considering 
graviton and dilaton propagation through the AdS-Schwarzschild black hole background. 

Finally, we compare our results to perturbation theory and reach some conclusions for the heavy ion 
collision experiments in section [4.41 

A few passages of this section are adopted from ref. [2\ as they stand. The phenomenological input 
and perturbative calculations as well as the numerical calculation of the Green functions by holographic 
methods where performed by D. Teaney, K. Dusling and C. Young during our collaboration on ref. 21. 
The main contribution of the author of this work was the AdS/CFT calculation in section [4.3.31 which 
is described in detail. 

4.3.1 EfTective model for heavy meson diffusion 

We make use of a model for heavy mesons and their interaction with the quark-gluon plasma, which 
was introduced in ref. [146| . This effective model describes the interaction with the medium by a dipole 
approximation. It relies on the large mass of the meson relative to the external momenta of the gauge 
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fields, i. e. the momentum scale given by the temperature of the medium, but does not rely on the 
smallness of the coupling constant. It was used previously to make a good estimate for the binding of 
J/ip to nuclei [146^. 

Because the model does not rely on the weakness of the coupling constant, we can make use of it in 
both the strong and weak coupling regime. At weak coupling we will refer to results from perturbation 
theory, while the results at strong coupling can be calculated from holographic duals. Since the exact 
dual to QCD is not known, we once more have to be satisfied with results ioi M — A SYM theory. 
Therefore, we have to rephrase the model in terms of supersymmetric fields. 

Diffusion in large N QCD 

The heavy meson field 4> describes a scalar meson which has a fixed four- velocity = (7,7V). Then 
the effective Lagrangian for this meson field interacting with the gauge fields is jl46) 

Leff = ~(j>hu ■ 90 + ;^<^^Ob0 + (4.78) 

where we refer to the last two terms as the interaction Lagrangian Lint , and 

Oi; = -^F'^™F,""ii^u,, (4.79) 

Ob = u^u, + -^F'^P'^F^p'^ . (4.80) 

Here, F is the non-Abelian field strength of QCD, with Greek letters /z, i/, . . . as Lorentz indices and 
gauge index a, The Ce and Cg are matching coefficients (polarizabilities) to be determined from the QCD 
dynamics of the heavy quark-antiquark pair. In inserting a factor of \/N'^ into the effective Lagrangian 
we have anticipated that the couplings of the heavy meson to the field strengths are suppressed by N"^ 
in the large limit. 

In the rest frame of a heavy quark bound state with u = (1,0) the operators Oe and Ob simplify to 

Oi5 = iE"-E\ (4.81) 

C>B = ^B''-B% (4.82) 

where E" and B° are the color electric and magnetic fields. If the constituents of the dipole are non- 
relativistic it is expected that the magnetic polarizability cb is of order O(v^) relative to the electric 
polarizability. For heavy quarks, where Cb is neglected, and large N these matching coefficients were 
computed by Peskin [14711148) . 

CB = ^, Cb=0. (4.83) 

Here A^ :~ l/flo = {'mq/'^)CFOts is the inverse Bohr radius of the mesonic bound state. It is finite at 
large A^ since with Cp — N/2 and finite A we have Kb = mqA/(167r). 

The effective Lagrangian can be used to calculate the in-medium mass shift. We will do so in the 
subsequent by simply consulting first order perturbation theory which says that 

6M = (iJint) - - (Lint) . (4.84) 



Translating the model to = 4 Super Yang-Mills theory 

Our aim is to calculate the heavy meson diffusion coefficient from gauge/gravity duality. Subsequent 
to this subsection we explain the Langevin dynamics we use to describe this process, it requires the 
calculation of the two-point correlators as well as of the associated polarizabilities Ce and c^. Because 
we do not now the gravity dual to QCD we translate the effective meson model to A/" = 4 Super Yang-Mills 
theory, our standard toy model. 

The formalism in A/" = 4 SU(A^) Super Yang-Mills theory is not different from the one we introduced 
in the preceding section. In general all operators in A/" = 4 SYM which are scalars under under Lorentz 
transformations and SU(4) R-charge rotations will couple to the meson at some order. The contribution 
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of higher dimensional operators is suppressed by powers of the temperature to the inverse size of the 
meson. The lowest dimension operator which could couple to the heavy meson field is Ox^ = TrX*X', 
where X' denotes the scalar fields of the theory. However, the anomalous dimension of this operator is 
not protected, and the prediction of the supergravity description of A/" = 4 SYM is that these operators 
decouple in a strong coupling limit [8]. The lowest dimension gauge invariant local operators which are 
singlets under SU(4) and which have protected anomalous dimension are the stress tensor 7^^ which 
couples to the graviton, and minus the Lagrangian Op2 = — L7V'=4, which couples to the dilaton. (Since 
we can add a total derivative to the Lagrangian, the operator — L is ambiguous. The precise form of the 
operator coupling to the dilaton is given in ref. |149) . We neglect this ambiguity here.) There also is the 
operator Of*f = TrF^''*F^,y + . . ., which couples to the axion. An interaction involving Of*f breaks 
CP-symmetry, which is a symmetry of the Lagrangian of the Af — 2 hypermultiplet of the A/" = 4 SYM 
gauge theory. Thus interactions involving 0f*f can be neglected. 

Summarizing the preceding discussion, we find that the effective Lagrangian describing the interac- 
tions of a heavy meson coupling to the operators in the field theory is 

Leff - -(/.^ {t, ^)iu ■ d(t>{t, x) + ^ 0^ {t, x) Or (t>{t, x) + ^ 0^ {t, x) Of- cj,{t, x) , (4.85) 

which is a linear perturbation of A/" = 4 Super Yang-Mills theory. The two composite operators Oj- and 
O F- in the interaction Lagrangian Lint are 



=0 



Or - r^" u^u, i V", (4.86) 
Of2 = F^'^F,^. (4.87) 

They account for the interaction of the mesons with the background. In gauge/gravity duality the 
modification of the Lagrangian described by Oj- is achieved by considering the AdS-Schwarzschild black 
hole background where (Of-) — 0. On the other hand, a finite (Of-) ^ is dual to a non-trivial dilaton 
fiow described by Liu and Tseytlin in ref. [145] . Details follow below. 

The polarization coefficients Cr and Cp will be determined below from meson mass shifts in gauge/gravity 
duality. This requires breaking some of the supersymmetry. We work in the linearized limit of small 
contributions from Oj- and Of-- This allows to investigate the effects of finite temperature and back- 
ground gauge fields separately. Additionally, this justifies the use of first order perturbation theory to 
compute the meson mass shifts in the medium as above by setting SM — — (Ljnt). For the contribution 
of the energy-momentum tensor, this is achieved by switching on the temperature. Then, the mass shift 
of the meson is given by expectation value of the stress tensor. Again we consider the rest frame of the 
mesons, 

SM^^^irn, (4.88) 
In contrast, for the meson response to (Of-) the mass shift of a heavy meson is given by 

SM = --^ {Of-} . (4.89) 



Langevin dynamics 

We now turn to the kinetics of the slow moving heavy meson with mass M in the medium. The kinetic 
energy i?kin = pv /2 of the meson can be assumed to be of order of the temperature T of the medium, 
such that pv ~ T. With p — Mv we can estimate the velocity and momentum to be 



p« Vmt, 




(4.90) 



For time scales which are long compared to medium correlations, we expect that the kinetics of the 
meson can be modeled as Brownian motion and can be described by Langevin equations. These are valid 
for times which are long compared to the inverse temperature but short compared to the lifetime of the 
quasi-particle state. We model viscous force and random kicks in spatial directions Xi by 

^ = m - VnP. , im ^At')) = % S{t - t') . (4.91) 
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Here, is a component of the random force ^ with second moment k and rjo is the drag coefhcient. The 
solution for pi (t) is given by 

t 

p,it)^ /dt' e.(i'), (4.92) 



supposed that rjot 1 |150j . This ahows to relate the drag and fluctuation by 



3MT = = f dhdt2 e"°(*i+*^) (e.(t2) Ut2)) - (4.93) 

— oo 

This leads to the Einstein relation 

(4.94) 

One of the aims of this section is the calculation of the diffusion coefficients or k, equivalently. 
From (I4.9ip we can obtain these coefficients once we know the microscopical phenomenological force 

:F,{t) = ^ (4.95) 

acting on the quasiparticle state. We can then compare the response of the Langevin process (j4.91l) to the 
microscopic theory (|4.95p . Over a time interval At which is long compared to medium correlations but 
short compared to the time scale of equilibration we can neglect the drag, which is small for the heavy 
meson with rjn oc 1/M. Since the considered time interval is long compared to medium correlations we 
can however equate the stochastic process, the random kicks ^, to the microscopic theory. We average 
(HUD 

Jdt Jdt' (6(0 ^j{t')) = AtnS.j = Jdt Jdt' {F,{t)F,{t')) . (4.96) 

At At 

In a rotationally invariant medium we have for i — j 



dt {Tj{t)Tj{0)) . (4.97) 



We now identify the force with the negative of the gradient of the potential V that we read off from the 
Lagrangian or our theory, i.e. the interaction Lagrangian V = —Lint- For the case of QCD with only 
Oe switched on we get 

J^{t) = J d'x 0t(t,x) ^VO£(t,x) <^(t,x), (4.98) 

which is the usual form of a dipole force averaged over the wave function of the meson. 

In our case k is a constant in space and time, i. e. we consider situations with constant diffusion 
parameters in a homogeneous medium, for instance slight deviations from equilibrium. From the point 
of view of a more general description in Fourier space with k{uj) we therefore are only interested in 
the hydrodynamic limit of a; — >■ 0. The fluctuation dissipation theorem relates the spectrum of (|4.97p 
(with the specified time order of operators) to the imaginary part of the retarded force-force correlation 
function cx (J%, (i) J^, (0)) on the right hand side. In the hydrodynamic limit we get 

1 2T 

— lim — ImG"(w), (4.99) 
3 w-i-o ll) 

where the full form of the retarded correlator is 

= -i [ dt e+'"* e{t) ([J-j (t) , J-j (0)]) . (4.100) 



Integrating out the heavy meson field as discussed in detail in ref. [lOOj . which treated the heavy quark 
case, we obtain a formula for the momentum diffusion coefficient 

1 4 f d^g ^2 (^_^iniG«(t.,q)) , (4.101) 



3 iV4 7 (27r)3 \ oj 
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with the retarded OeOe correlator given by 

G^(w,q) = -z j A-'x e+''^'-'^-^ 6{t) {[OE{t,Q) ,Oe{QM) ■ (4-102) 

We can understand this result with simple kinetic theory. Examining the Langevin dynamics we 
see that 3k is the mean squared momentum transfer to the meson per unit time. The factor of three 
arises from the number of spatial dimensions. In perturbation theory this momentum transfer is easily 
computed by weighting the square of the transferred momentum of each scattering with the transition 
rate for any gluon in the bath to scatter with the heavy quark, 

Here, p is the spatial momentum of the incoming gluon, p' is the momentum of the outgoing gluon and 
q is the momentum transfer q = p — p', and \A4\^ is the gluon meson scattering amplitude computed 
with the effective Lagrangian in (j4.78l) and weighted by the appropriate momentum distributions n of 
the incoming and outgoing gluons. 



IMl'^^u^ (l + cos2(V))- (4.104) 

Alternatively (as detailed in appendix A of ref. |2j), we can simply evaluate the imaginary part of the 
retarded amplitude written in (j4.10ip to obtain the same result. 
In A/" = 4 theory the generalized force is given by 

J^{t)^- yd3a;(/.t(i,x)v(^0T(i,x) + ^0;..(t,x))</)(t,x) (4.105) 

which results in a momentum broadening 

1 r c\^n IT f \ 

^ = -3 iTo/ (2^^V(#I^^^-(^'^) + #I"^^-(-'^V' ^'-'"'^ 
where the retarded correlators at vanishing velocity are 



G^r^-^ yd*xe+*"*-''i-"e(t)([r"°(t,x),r°"(0,O)]), (4.107) 
G^F - / d'^^; e+*"*-*i-'' 9{t) {[Op^ (t, x) , Op^ (0, 0)]) . (4.108) 



In writing (I4.106P we have implicitly assumed that there is no cross term between Op^ and O-j-. In the 



(5 5- 

gauge/gravity duality this is reflected in the fact that at tree level in supergravity ^^oo (_x)F<^{y) 
4.3.2 Weak coupling — perturbative results 

We begin with the results for perturbative QCD (pQCD). The mass shift is obtained from first order 
perturbation theory as SM = (iJint) = — (Lint), yielding 

JMpQCD = {Oe)t 

^ /^yi4 (4.109) 
Us J 45- 

In the second line we have calculated the thermal expectation value {Oe)t = fo-^^-^^ ™^ ^ ^""^^ gluon 
gas and used (|4.83p . 

The importance of this result is that it is finite at large N and that it is in general suppressed 
by {T/Ab)^j i.e. by powers of the hadron scale to the temperature. If higher dimension operators were 
added to the effective Lagrangian their contributions would be suppressed by additional powers of T / Kb- 



0. 
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For QCD the integrals written in (|4.103p are straightforward and yield the following result for the 
rate of momentum broadening 



1 2 6477^^ 



T3 /ttTV 501767r 



6 (4.110) 



Af2 \AbJ 1215 ■ 

The high power of temperature arises since the dipole cross section rises as w'*. The matching 
coefficient is directly related to the mass shift of the dipole and the inverse Bohr radius by (I4.83P 
and (|4.109l) . It encodes the coupling of the long distance gluon fields to the dipole. By taking the ratio 
between the momentum broadening and the mass shift squared, we find a physical quantity which is 
independent of this coupling 

'''' (4.111) 



{6My 



pQCD 



The large numerical factor 1280/3 originates from the cross section which grows as lu'^. A similarly large 
factor appears in A/" = 4 SYM below. 

For comparison with the AdS/CFT result we list the results for A/" = 4 super Yang-Mills theory in 
the limit of small 't Hooft coupling A, again computed in ref. [5], 

SAh^o^c^^-, (4.112) 
2 623271^ 

^.^o^cl^-, (4.113) 

= ^36.9. (4.114) 



4.3.3 Strong coupling — holographic calculation 

We first determine the polarizabilities cp and ct from the mass shifts of the meson in two different 
backgrounds using (|4.88p and (|4.89l) . To accomplish this, we will switch on the perturbations of the 
A/" = 4 Lagrangian which correspond to finite Oj- and finite Op2 . Again we consider the linear limit such 
that we can investigate the effects of finite temperature and finite background field strengths separately. 

Subsequently we will compute the correlators in (j4.107l) and (|4.108l) for strongly coupled A/" = 4 theory 
at finite temperature. The results we obtain in this section will be put together in the next subsection 
using (I4.106[) to deduce the rate of momentum broadening and compare it to the weak coupling result. 



Backgrounds dual to finite temperature and field strength 

The gravity background dual to A/" = 4 SYM theory at finite temperature is given by the AdS- 
Schwarzschild black hole with Lorentzian signature. This background is needed below both for cal- 
culating the necessary two-point correlators {Or Or) and {Op20p2), as well as for obtaining the meson 
polarizability ct, which accounts for meson mass shifts due to finite temperature. 

We make use of the coordinates derived in appendix |X] as (IA.3I) to write the AdS-Schwarzschild 
background in Lorentzian signature as 



d.^ = J (-^ dt^ + /» dx^) + 5 (d.^ + vW,) 



(4.115) 



■'^ ' 4w4 ' ■'^ ' 4^4 



In this way we can identify the transverse part to Minkowski space as nothing else than and we write 
it as 

6 

dv^ + v'^dnl = = + y^dnj + dvl + dvl . (4.116) 

IR''("l....,4) IR2(«5,6) 
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with the metric dfi| of the unit 3-sphere, and v'^ = y^+vf+Wg. The boundary is reached at asymptotically 
large y while the horizon is located at ro/\/2. Notice that the black hole radius Tq is related to the 
expectation value {T^^) by jl51j 

(ro") In^N^T^ , To = TnR^ (4.117) 
8 



The field configuration dual to {Op^) 7^ and {T^^^) = is a non-trivial dilaton background with 
has been given by Liu and Tseytlin and consists of a configuration of D3-branes with homogeneously 
distributed D(— 1) instantons |145] . The type IIB action in the Einstein frame for the dilaton $, the 
axion C, and the self-dual gauge field strength F5 — reads 



1 



'IIB 



'9 



The ten-dimensional Newton constant is given by 

1 27r 



(4.118) 



(4.119) 



As can be seen from the supersymmctry transformations of the A/" = 4 fermions, such a background 
breaks the supersymmctry to M — 2. Solving the equations of motion derived from (I4.118p . Liu and 
Tseytlin obtain the metric ^145] 



ds^ 

"^■^ string 



^ "^''Einstein 



(dr^ + r^dVt 



and axion-dilaton solutiord 



e* = l 



The expectation value {O p-i ) is given by 



{Op2) = lim 



r->oo (5<i)(r, x) 27r2i?8 



(4.120) 
(4.121) 

(4.122) 



Computing the polarization coefficients from meson mass shifts 

We again identify mesons with fluctuations of a D7-brane embedded into the background dual to the 
field theory under consideration. Stable embeddings are obtained if the D7-brane spans all Minkowski 
directions as well as the radial AdS coordinate and a 3-sphere in the remaining angular directions. 
Consider the metric (|4.115l) with (I4.116P as an example. Here, the D7-brane shall be embedded such 
that it spans all directions except W5 and vq. The meson mass M is then obtained by solving the equation 
of motion for the fluctuations (p [43] , as outlined in section 12.31 Read as an eigenvalue equation, the 
equation of motion for the fluctuation gives the meson mass as the eigenvalues M to the corresponding 
eigenfunctions (p. The discrete values of M describe the Kaluza-Klein mass spectrum of mesons for any 
given quark mass. 

To see this explicitly and generalize to the backgrounds of interest below, we rephrase this procedure 
for the vacuum case (T™) = (Of^) = in a notation suitable for the subsequent generalization. 
Subsequently we will introduce a non-zero {Op^) and (T'^^), respectively. 

In the case of a D7-brane embedded in a ten-dimensional background, the brane embedding is de- 
scribed by the locations and vq in the two directions transverse to the brane. This setup was introduced 
in section B.2.21 In general these locations depend on all eight coordinates of the eight-dimensional 
D7-brane worldvolume and are determined by extremizing the DBI-action (j2.65p . We rephrase it here 
including the dilaton which we where free to set to unity in (j2.65p . 

d«e e-*y^ditG , Gab = -Q^-Q^ 9%. (4.123) 

^Regarding conventions, note that in our notation q = ^<?lt, wliere Qlt is used in the paper of Liu and Tseytlin |145| . 
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with Tj the D7-brane tension and g^^ is the string frame metric of the ten-dimensional background with 
coordinates X^. It is related to the Einstein metric as in (|4.120p . The distinction between the Einstein 
and string frame is ultimately important below to account for the effects of the non-trivial dilaton flow. 
The puUback G contains the functions W5(0 ^J^d ve{^), which are determined by solving their equations 
of motion, derived from S'dbi- 

The background AdSs x S^dual to (r°°) {Op^) = is obtained e.g. from (|4.115l) with ro = 0. It 
is well known that for this background a probe brane embedding is given by the functions 

V5 = 0, (4.124) 
vq — niy = const, (4.125) 

and the constant rriy determines the quark mass ruq — my / {2'Ka'). In terms of the example in section [!j.2.2l 
this is instantly derived from the equation of motion (|2.72p for the embedding L in the zero temperature 
limit To — >■ 0. Note however that there is a factor \/2 differing in the definition of the quark mass due to 
the choice of coordinates. 

We are interested in the meson spectrum, which can be obtained from the brane fluctuations as in 
the previous sections, c. f. for instance page [27] in section 12.31 We thus allow for small fluctuations if> 
around this solution. Here we consider the fluctuation of the radial part tig, dual to the scalar meson 
excitations, 

Vf, ^ vq + 2iTa (f{x,w). (4.126) 

By the symmetries of the setup, the fluctuations only depend on the Minkowski directions x and on 
the coordinate y, denoting the radial coordinate on the part of the D7-brane which is transverse to 
the Minkowski directions. The resulting equation of motion is analog to formula (|2.78p (modulo the 
mentioned factors of \/2). The solutions were found by plugging in the ansatz 

(^ = ¥'(y)e^''">'(S3), (4.127) 

where y'(S'^) are the scalar spherical harmonics on the S'^ wrapped by the probe D7-brane and k denotes 
a four vector. The resulting equation of motion for the function ip{y) for I — reads 

- dpp^d^p) = (^2 + 1)2 M^vip)- (4.128) 

Here we introduced the following dimensionless quantities 

V - 

p=—, M= M, (4.129) 

rriy niy 

where the dimensionless p is not to be mistaken with the same symbol of different meaning defined 

in section [31 We moreover identified the meson mass squared with the square of the momentum 
four-vector k of the fluctuations, 

= -P. (4.130) 

The eigenfunctions ipn solving the Sturm-Liouville equation (|4.128l) are given in terms of the standard 
hypergeometric function 2F1, 

^-^P) = (^2 + i)n+i (-(« + i);-";2;-p'), (4.131) 

where n = 0, 1, 2, . . . and Cn is a normalization constant such that 

oo 

'^P ^p2 2)2 ^P^ ^P^ ^^rim- (4.132) 


The lowest mode </3o is given by 

Mp) = ^1- (4.133) 
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The corresponding eigenvalues A/„ to the functions ipn are given by 



A?„ = 2v/(" + !)(" + 2). (4.134) 

We note that the mass of the lowest state with ?i = is 

m„ r- Ana' r- 47rm„ 

Afo = -^272 = 771,-^72=-^, (4.135) 

which will appear frequently below. For a more detailed derivation of these results the reader is referred 
to ref. |43j. 

Mass shift in the dilaton background Let us now calculate the polarizability Cf which determines 
the change SM of the meson mass at a given value of the gauge condensate {Op2 ) with respect to the 
meson mass at {Op2) = 0, 

To find Cp we will determine the mass shift SM and identify Cp with the proportionality constant in 
front of {Op2). 

We are interested in the eigenvalues of fluctuations in the case oi q (x {Op2) ^ 0. The ten-dimensional 
background geometry dual to this scenario is given in (j4.120p and the equation of motion for D7-brane 
fluctuations analog to (j4.128p was derived in ref. [lS2. to be 

- - ^(^) - 4^" (,. + l)(/+(p^ + l)^) ^'-'"'^ 

with the dimensionless 

g- = . (4.138) 

To obtain analytical results, we consider the case of small q and linearize in this parameter. Therefore 
the equation of motion to solve is 

- dpp'dMp) = M^J^^J^ Vip) + MpMp): (4.139) 



where the operator A(p) is given by 



A(p) = -49p^9,. (4.140) 



It is this term that describes the difference between the equation of motion at non vanishing background 
perturbation to (|4.128p . which is reproduced for q = 0. 

To flnd the solution (po{p) corresponding to the lightest meson with n = we set up a perturbative 
expansion. Any deviation Stpo from the solution ipo of the case q = may be written as a linear 
combination of the functions ipn, which are a basis of the function space of all solutions, 

oo 

4>{p) = Mp) + "-nMp), an < 1, (4.141) 

n=0 

M2 ^ ^ gj^2^ g]^2 ^ ^ (4.142) 

Plug this ansatz into the equation of motion (I4.137p . make use of (j4.128p and keep terms up to linear 
order in the small parameters a„, q and 6Mq to get 

2 00 3 ;^ 00 

J^TW X] °"^^"'^"('°) ^ '^^^" (p2 + i)2 ^o{p) + A^0 (^2%)2 E «" ^"('°) + ^ip)Mp)- (4.143) 

We now multiply this equation by ^o{p), integrate over p £ [0, 00] and make use of (|4.132D and (j4.133p 
to see that 

00 

6M^ = - I dpMp)Mp)Mp) = -^9. (4.144) 
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From 5M^ = 2Mq5Mq we obtain 

where in the last step we used (I4.135P for the mass and (|4.122p and (|4.138p to relate q and Op2. By 
comparison with (j4.136p we identify the polarizability 



8 / 27r ^ ^ 



57r V M, 



(4.146) 



Mass shift in the finite temperature background The calculation of the polarizability Ct is 
completely analogous. We are now looking for the proportionality constant of meson mass shifts with 
respect to deviations from zero temperature, which we noticed to be given by 

5M = -|j(r°"). (4.147) 

The background dual to the finite temperature field theory is the AdS black hole background given in 
(|4.115p with (|4.116p . 

Again we calculate the meson mass spectrum to identify the polarizability by comparison with (I4.147P . 
The embedding functions and vq in this background are given by 

V5 = 0, (4.148) 
VQ = VQ{y), (4.149) 

where the quark mass is determined by mq = limj,_j.co vq/ (27ra'). The function VQ{y) has to be computed 
numerically |37) . Some examples of such embeddings are shown in figure 12.41 
We introduce small fiuctuations tf{p)e'^^'-^ in the W5 direction, 

V5 ^ vz{y, x) + ^{y)e'^^ . (4.150) 

The linearized equation of motion for the fluctuations ip{y) in the limit of vanishing spatial momentum 
and = — fc^ can be derived from the DBI action (|4.123p to be 



0= dy 



1 + [OyVQ) 



((y2 +t,^)2 -^4) 



(4.151) 



where we abbreviated 

g^Q^l ^° j\. (4.152) 

V 16(y^+^i)V 

In the regime of small temperatures, we may linearize in which is the leading order in Tq. Further- 
more, as may be seen from figure I^T^ in the regime of a small temperature T compared to the quark mass 
TOq, or respectively small ratios of Tq/ limj,_^oo ?Je(j/), the embeddings become more and more constant. 
So for constant embeddings vq = and up to order T** cx the equation of motion simplifies to 

- dpP^dMp) = ^'^^^J^^viP) + ^(P)^(P). (4-153) 
where we made use of the dimensionless quantities (|4.138l) and identify 
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For the lightest meson, the ansatz (|4.14ip this time leads to 



oo 

t2 



5M^ = - dpipo{p)A{p)ipo{p) 

(4.155) 



'o 



40 n4 

Reinstating units and solving for SMq leads to 



12 / 27r \^ 1 



From this we can read off the polarizability Ct as 



3 



Computing finite temperature correlators 

According to (I4.106P we need to compute the correlators (|4.107p and (|4.108l) at finite temperature. We 
do so by once more performing the calculation along the lines of refs. (22^,23^, sketched in section B.2.11 
The dual supergravity field to the energy momentum tensor T is the graviton h, and the corresponding 
field to the operator Op2 is the dilaton $. Therefore, correlators of 7^° are associated with graviton 
propagation and obtained from the supergravity field solution to h. Respectively, correlators of Op2 are 
associated with dilaton propagation and obtained from the supergravity solution to the dilaton The 
calculational procedure for these two correlators is standard and has been discussed in refs. ^511126) . and 
first applied in ref. [22], in order to find two-point Minkowski correlators as discussed in section [2. 2 .11 

On the gravity side, both field correlators are computed in the black hole background to account 
for finite temperature correlation functions of the dual gauge theory operators. In this subsection it is 
convenient to work in the coordinates derived as (IA.4I) . with radial AdScoordinate u. For an explicit 
calculation we are more specific here as in (|2.6ip by writing 

G"(w, q) = -^ / d\ e-'^-^9{t) {[Op^ , Op^]) 

■> ^ (4.158) 

AdS/CFT ^^^^ ^^^^ _ duf{u, k). 

Here once more the function /(u, k) relates the boundary and bulk values of a gravity field to each other. 
For example the dilaton field $ is related to its value at the boundary (1)°'^^ by 



<^{u,k)^f{u,k)cf}'''y{k), (4.159) 

and / is normalized to one at the boundary, i. e. /(O, fc) = 1. For fiuctuations hP'^ of the metric component 
(7*^°, $(m, k) is replaced by h^^{u, k) and the field theory operators Op^ in (|4.158p are replaced by T*^*^. 
The factor A{u) can be read off from the classical supergravity action 

Sci^ dud-^xAiu) (9„$)2 + . . . . (4.160) 

The classical five-dimensional gravity action for the graviton and dilaton is obtained from (j4.118l) as 



^ci = ^ y dud^x^, (^(7^ - 2A) - i(a$)2 



(4.161) 



where 

L-^-Jfl. (4162) 
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So comparing to (|4.160p we get 

A^-^g--. (4.163) 
The equation of motion derived from (j4.16ip in momentum space reads 

with dimensionless frequency tn = uj/{2ttT) and spatial momentum component q ~ q/{2T:T). The prime 
denotes the derivative with respect to the radial coordinate u. Note that in momentum space the function 
<I>(m) is akin to the radial part we denoted by f{u) above. The equation of motion (j4.164p has to be solved 
numerically with incoming wave boundary condition at the black hole horizon. Computing the indices 
and expansion coefficients near the boundary u ~ 0, as done in refs. |64[ 165] . we obtain the asymptotic 
behavior as a linear combination of two solutions ^1^2 with asymptotic behavior for small u as 

$i = (l + ...), (4.165) 

$2 = (w^ + ...)■ (4.166) 
The general solution to the equation of motion therefore is given by the linear combination 

$(w) = $i+6$2, (4.167) 

where we normalized the functions $ such that the coefficient for $1 is 1, which we are free to do since 
(|4.164p is a homogeneous equation. In this way we achieve that the correct normalization for the radial 
part lim„_>.o ^(u) = 1 is implemented. 

At the horizon the asymptotic solution satisfying the incoming wave boundary condition is 

$(w) (1-m)-"5'(1 + ...). (4.168) 

As discussed in refs. [Ml [53] we find the coefficient B by integrating the two boundary solutions 
from (j4.167p forward towards the horizon and by matching the linear combination of the numerical 
solutions $(m) — $5"^™ + B^2^™ to the solution (j4.168p at the horizon. We recognized that the solution 
of the radial part $(u) found in this way is equivalent to the radial part of the total solution <1> and can 
therefore be plugged in for / in (j4.158p . 

The imaginary part of the retarded correlator then is given by 

ItclGpp = —5 . (4.169) 

Solving (|4.164D and matching the asymptotic solutions as described above, thus enables us to obtain 

l™n /t^^ f-^ I™ =^'^'67.258. (4.170) 

The corresponding result for the energy-momentum tensor correlator is obtained in an analogous way 
but the analysis is significantly more complicated. Fortunately it has been extensively and carefully 
analyzed in ref. [126) . The final result is 

lim /-^^ f^Im Gf^(w,g)) ^iV^T^ 355.169. (4.171) 
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4.3.4 Comparing vireak and strong coupling 

We now have the results to compare momentum broadening of the heavy meson in a hot medium at weak 
and strong coupling. Over the duration of the lifetime of the heavy meson state it will loose momentum 
on average and simultaneously receive random kicks as codified by the Langevin equations of motion 
(|4.9ip . The drag and momentum broadening rates are related by the Einstein relation 
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with Mq the meson mass. 

For strongly coupled A/" = 4 SYM theory we obtain our principle result by collecting the results for 
polarizabilities (|4.146p . (|4.157p and force correlators (I4.170p . (|4.17ip . and use (|4.106p . 



«A^=o = — ( — I (^1 67.258+1-1 355.169 
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) 224.726. 

The mass shift in strongly coupled M — ^ SYM is given by the sum of the mass shift due to dilaton 
and graviton contributions, respectively. However, the exact value of the dilatonic mass shift (I4.136P is 
determined by (|4.122p in terms of g, which we do not want to speculate about here. Nevertheless, from 
the definition of q in ref. |145| we know that it is positive for positive instanton numbers. In this case 
the dilatonic mass shift contributes with the same sign as the graviton mass shift (I4.147P and we write 



5Mo 



-T 



f2.TV%. (4.174) 
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Here we made use of relations (|4.157p and (I4.117P . Comparing these formulas with the analogous formulas 
in weak coupling large N QCD given in equations (j4.110|) and (|4.109p . 

r3 (^T\'' 50176 ,^ 
-pQco = ^(^) (4.175) 



and 



JiVWo = -T ( ^ y , (4.176) 



we see that the meson mass Mo plays the role of the inverse Bohr radius As = {mq/2)asCF in the 
strong coupling dipole effective Lagrangian. This is as expected for relativistic bound states. Since these 
prefactors are different we do not compare the numerical values. Below we will compare the values of 
the ratio k/{5M'^) at strong and weak coupling. We moreover observe that the drag coefficient for heavy 
mesons is suppressed by N'^ in the large N limit. 

The phenomenological model we used modeled the interaction of the mesons with the medium as 
dipole interaction terms, which capture short distance phenomena. From the field theory side one may 
wonder if these dipole interactions indeed are the dominant interaction mechanism of the medium and 
the mesons. Considering AdS/CFT, we admit that a dipole picture of a meson is a short distance 
description an we can not guarantee to describe such UV effects of the field theory well by holographic 
models. On the other hand, even if we cannot describe the interaction from first principle, the diffusion 
of the mesons and the scattering of gluons which then propagate with modified momentum into the 
medium is a long distance effect, which we can hope do describe by our approach. We therefore consider 
a quantity that is independent of short scale dipole interpretation of the medium, which is parametrized 
by the connection coefficients c. To construct such a quantity, we remember from (|4.106p that (for the 
simple case of Cp = 0) 

'-f „„J#^f£fz^„„Gf,(„.c,), ,4,177) 



Kn'^J JTo J (27r)3 3 V 

Together with (|4.88p we can then construct a quantity that does not depend on the connection coefficient 

K 1 ,. /■ /'-2T 



,r olim/^^— ImG^T-(w,q) . (4.178) 

((5M)2 (7-00)2 7 (-27r)3 3 V TT\ J V ; 

It is independent from the connection coefficient . AH quantities in this expression have been calculated 
previously for weak coupling as well as for strong coupling. The strong coupling result is obtained from 
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inserting the relations (|4.117p and (|4.171l) and yields 



((5M)2 



T,A->oo 



(4.179) 



To include the dilaton contribution, we make use of the above results (|4.173l) and (j4.174p to get 



{5Mf 



(4.180) 



It is then reasonable to use AdS/CFT to estimate to what degree strong coupling physics modifies this 
ratio in QCD. In the free finite temperature A/" = 4 theory, the result was computed in appendix A of 
ref. [2] and cited in (|4.112p to be 

(4.I8I) 



(<5M)2 



iV2 



36.9. 



Thus, comparing the strong coupling result (j4.180p with the weak-coupling result (|4.I8ip . we conclude 
that strong coupling effects actually reduce the momentum transfer rate relative to the mass shift. 



4.4 Summary 

We studied transport properties of baryon charge, isospin charge and heavy mesons with three different 
methods. For the charge diffusion we derived the diffusion coefficient in dependence on temperature and 
particle density (or chemical potential, equivalently). In the case of meson diffusion in the quark- gluon 
plasma we observed how strong coupling effects the equilibration. 



Baryon diflFusion We made use of the most naive formulation of the membrane paradigm and derived 
a dependence of the diffusion coefficient for baryons that was qualitatively confirmed by recent and more 
comprehensive studies fll9j . The results, shown in figure BTTl indicate that for very large values of the 
baryon density, the diffusion constant asymptotes to its maximal value D = 1/{2'kT). This reflects the 
fact that in this case, the free quarks outnumber the quarks bound in mesons. For flnite densities the 
diffusion coefficient is reduced and we observe a minimum at values of the quark mass to temperature 
ratio TO, which lie above the fundamental phase transition at zero density. We will draw a connection 
between the phase transition and the behavior of the diffusion in section [S] 



Isospin diffusion To obtain the isospin diffusion constant we analyzed the dispersion relation of 
the lowest lying quasi normal mode of the isospin current. We have considered a relatively simple 
gauge/gravity dual model for a flnite temperature fleld theory, consisting of a constant isospin chemical 
potential /i obtained from a time component vacuum expectation value for the SU(2) gauge field on 
two coincident brane probes. We have considered the constant D7-brane embedding corresponding to 
vanishing quark mass. 

Within the strong restrictions of our analytic derivation, the main result is that this model, despite 
its simplicity, leads to a hydrodynamical behavior of the dual field theory which goes beyond linear 
response theory. We find in particular a frequency-dependent diffusion coefficient with a non- analytical 
behavior. Frequency-dependent diffusion is a well-known phenomenon in condensed matter physics. Here 
it originates simply from the fact that due to the non-Abelian structure of the gauge field on the brane 
probe, the chemical potential replaces a time derivative in the action and in the equations of motion 
from which the Green functions are obtained. For a more comprehensive study one should include terms 
of quadratic order in the chemical potential m, which then would cancel the non-analytic behavior. 



Quarkonium diffusion In the last part of this section we studied the diffusion of heavy mesons in 
the holographic plasma, by setting up a Langevin model. The forces on the meson where deduced from 
an effective dipole model for mesons that allowed to carry out our computations at strong and weak 
coupling. We therefore where able to estimate the effects of strong coupling on meson diffusion. 

On the field theory side, we deduced the meson couplings to the stress tensor T and the operator Opi 
at strong coupling, which are the only relevant operators coupling to the heavy dipole. The couplings 
where deduced from mass shifts, which can be computed on the gravity side of the correspondence as a 



80 



change in the normal vibrational modes of the D7-brane in the presence of an external gravitational, or 
respectively dilatonic, field. 

Because the gravitational and dilatonic fields shift the spectrum of the D7-brane excitations, gradients 
such as those from fiuctuations of these fields give rise to a net force on the mesonic modes. We used 
the fluctuation dissipation theorem to relate the spectrum of these fluctuations in the plasma to the 
momentum broadening of the meson. 

The result for «;/(5M)^, the momentum broadening relative to the square of the in medium mass 
shift, for strongly coupled A/" = 4 SYM theory is roughly four to five times smaller than the result in the 
weakly coupled limit. We therefore conclude that in this model the effect of strong coupling is to reduce 
the momentum broadening and increase diffusion of mesons relative to weak coupling. 

From a phenomenological perspective the current calculation was limited to very heavy mesons, 
which survive above Tc, where dipole interactions between the meson and the medium are dominant. 
It is certainly unclear if this is the relevant interaction mechanism above Tc even for bottomonium. 
Furthermore, the dipole coupling between a heavy meson and the medium is dominated by short distance 
physics which is not well modeled by AdS/CFT. 

From the supergravity perspective, a better understanding of how gravitational and dilatonic fields 
fiuctuate in bulk would give a straightforward procedure to calculate the drag of a finite mass meson. 
Specifically, fiuctuations in the bulk would force motion of meson wave functions which extend into the 
holographic fifth dimension. 

5 Exploring the phase diagram 

While the theory of quantum chromodynamics is concise in its mathematical formulation, it exhibits a 
rich phenomenological structure, including several phase transitions. For instance, we frequently referred 
to the confinement / deconfinement transition quark matter is supposed to undergo when it is heated up 
and/or exposed to high chemical potential. 

Since the development of QCD it was discovered that quark matter can exhibit numerous qualitatively 
different behaviors. The most prominent example in QCD is the change in the coupling constant (7ym 
of the strong interaction with respect to the momentum scale at which the theory is probed [1531 1154j . 
It eventually accounts for the transition from the zero temperature regime (ground state with mini- 
mum momentum), where quarks are confined to the hadronic color singlets that make up the nuclei of 
atoms, to a deconfined state of matter at asymptotically high temperatures (and high thermal momenta) 
where quarks and gluons roam freely throughout spacetime. The value of the momentum interchange 
between quarks and gluons determines their interaction potential, which may change qualitatively, e. g. 
by increasing temperature, from a confining shape of infinite depth to a potential with finite binding 
energy. 

Another example is the discovery of the so called color-fiavor-locked phase of QCD with three color 
and three flavor degrees of freedom at high chemical potential, equivalent to high particle densities. 
It was observed that the thermodynamically favored ground state of QCD changes with the chemical 
potential. While we observe color and flavor symmetry separately for QCD at low chemical potential, 
there is some critical value of the chemical potential at which symmetry breaking occurs and a relation 
between color and flavor degrees of bound states of quarks is established |155l I156| . 

The motivations to explore the properties of QCD at high temperature and large chemical potentials 
include the struggle for a deeper understanding of the features of QCD itself, the formation of hadronic 
matter during the evolution of the universe and the description of matter inside dense astrophysical 
objects such as neutron stars. 

We may wonder which parameters determine if we observe confinement and/or color flavor locking, 
which other symmetries and observables undergo qualitative changes and whether these changes occur 
abruptly or smoothly in the parameter space. This raises the question about the structure of the phase 
diagram of the theory. It describes the regions in the parameter space of the state variables in which 
the thermodynamic potentials and their derivatives behave analytically. The change from one phase to 
another is often accompanied by symmetry breaking and can be described by an order parameter which 
assumes flnite values in one phase and vanishes in the other. For QCD, the quark condensate (^/''0) is 
frequently used as an order parameter. Non vanishing values break chiral symmetry. We can distinguish 
the hadronic phase, where ('(/'V') 7^ 0, from the quark-gluon plasma with a higher amount of symmetry. 



81 



T 

T 



QGP 



Tc 



\ 







phase ,,,,,, ,gj 



hadronic 




Figure 5.1: The conjectured QCD phase diagram in the plane of temperature T and baryon chem- 
ical potential /i. The black lines indicate first order phase transitions. The first order confine- 
ment/deconfinement transition ends in a critical point at unknown (/^c, 7c)- Here the transition becomes 
second order. Lattice simulations suggest 170MeV ^ < 190MeV, nuclear matter has [i ~ IGeV |158| . 
CFL is the color-flavor-locked phase. 

here (V'V') = 0. However, strictly speaking the chiral condensate is non-zero in both phases for finite 
quark masses and should only be used as an order parameter in models with vanishing masses for the 
light quarks. 

The exact overall structure of the phase diagram is not known yet. Especially at low temperature 
and chemical potential where the coupling is strong, theoretical treatments rely on lattice gauge theory, 
which on the other hand has its problems with modeling the QCD dynamics at finite temperature and 
finite baryon density. Figure ISTTI shows a sketch of the most basic features of the theoretically conjectured 
phase diagram of QCD based on combinations of analytical and numerical predictions |157l 1158] . 

In huge volumes that exist for long periods of time, such as the cubic kilometer volumes of quark 
matter in neutron stars, one can expect the thermodynamic limit to be an appropriate approximation 
to describe matter, i. e. thermodynamics is applicable. As hydrodynamic models have successfully been 
applied to describe the collective motion of the fireball produced in heavy ion collisions, it is reasonable 
to assume that even the quark-gluon plasma observed in experiments reaches thermal equilibrium. The 
equilibration time is estimated to approximately Ifm/c, the plasma state then exists for about another 
4fm/c |159| . Though in equilibrium, the multiparticle system in a collision experiment certainly does 
evolve along some trajectory in the phase diagram, which leads from some point in the QGP phase to a 
system of hadrons, which can eventually be detected. 

For a thermodynamic description, the state variables can be chosen to be for instance the temperature 
T, the pressure P, and the chemical potentials [ij for the different species of particles, labeled by the 
index j. The chemical potential /ij is the conjugate variable to the particle number Nj. It is therefore 
only well defined if the number of particles of the species j is well defined. 

There are only a few charges that are conserved by all standard model processes. They allow for the 
definition of particle numbers and thereby determine the parameters of the phase diagram. These are 
quark number (which can be translated into baryon number), lepton number, electric charge, and color 
charge. Each of them has a chemical potential associated to it. In the processes we discuss in heavy 
ion collisions, there are no leptons (and in neutron stars they are radiated off by neutrino emission), 
i. e. the lepton chemical potential can be set to zero. The system furthermore is color-neutral, i. e. the 
chemical potentials associated to color charge can also be set to zero. The Gell-Mann-Nishijima-relation 
allows to rephrase the remaining parameters, quark number and electric charge, in terms of isospin and 
hypercharge. Hypercharge in turn is determined by the number of particles of the various quark species. 

One may argue that in an equilibrium state weak interactions could account for flavor changing 
processes, and therefore there would be no well defined particle number associated to each quark fiavor. 
As a result, the notion of a chemical potential would not be well defined. In heavy ion collisions however, 
the system in equilibrium state has not enough time to undergo weak interactions and therefore the 
flavor numbers are conserved, and the different quark flavors have to be assigned individual chemical 
potentials. 
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Thus, the remaining degrees of freedom in the phase diagram are temperature and the chemical 
potentials of the interacting quark species. In two-flavor setups, which we elaborate on in this work, 
we are free to express these potentials in terms of the baryon and isospin chemical potential. This is 
why we used and continue to use these charges as the parameters throughout this work. In principle we 
therefore consider a three-dimensional phase diagram in (T, For sake of simplicity, however, we 
restrict to cases of either non- vanishing baryonic or isospin chemical potential. 

In the holographic context, several publications where dedicated to the investigation of the structure 
of the phase diagram of theories with gravity duals. For instance, the meson melting transition at finite 
temperature for the fundamental matter introduced by D7-branes in a background generated by D3- 
branes, was considered at zero as well at finite particle density. Studies of the behavior of D7-brane 
probes in the AdSSchwarzschild black hole background revealed a phase transition, which occurs when 
the D7-brane reaches the black hole horizon. This transition was shown to be of first order in ref. [35] , see 
ref. [52] for a similar transition in the D4/D6 system, further details can be found in refs. [55] [55]. Related 
phase transitions appear in refs. [TTl 11601 1161j . Subsequently this phase transition was investigated at 
finite particle density ^3 [Ml \S2\ ■ The transition is characterized by a stable quasiparticle spectrum in 
one of the phases and melting mesons at finite density in the deconfined phase. Figure IXTl illustrates the 
two phases. A way to derive the values of the transition line was mentioned in the same section, see the 
caption of figure and ref. [SI] for details. For the phase diagram in the model given by D8-D8 probes 
in a near horizon limit of D4-branes with non vanishing chemical and isospin potential see ref. [57] . 

In this section we contribute to these studies by different means. In the subsequent section we observe 
a phase transition in the baryon diffusion coefficient and determine the critical baryon density at which 
this transition vanishes. The observations we make can be related to the results from refs. [ISl and 
confirm the observations described therein. 

In section we describe the occurrence of a new phase transition at finite isospin density, which we 
published in ref. [3]- 

5.1 Phase transition of the baryon diffusion coefficient 

In this section we pick up the discussion from section [4. 11 where we introduced a simple yet incomplete 
derivation of the baryon diffusion coefficient. As we mentioned there, the qualitative behavior of the 
result captures the physics correctly and the quantitative evaluation bears only little error in the regime 
of masses up to the temperature scale [119| . On the other hand, the evaluation of the simple formula 
(|4.2p is a considerable simplification, compared to the exact treatment, described in ref. [119 . We will 
therefore stick with our approximation. 

We focus on the temperature regime near the phase transition of fundamental matter. In the zero 
density limit, a phase transition of fundamental matter was observed to take place simultaneously to 
the geometric transition from Minkowski embeddings at small temperature to black hole embeddings at 
high temperatures. In the low temperature phase the mesonic spectral functions exhibit discrete delta 
peaks, and thereby describe stable mesons, while the spectrum becomes continuous with finite width 
peaks in the high temperature phase. The value of the quark mass to temperature ratio at which the 
phase transition can be observed was found to be near m = 1.3, c. f. figure [mi The transition can be 
seen as a first order phase transition in the quark condensate |37| . In the zero density limit the authors 
of ref. j48j observed another manifestation of the fundamental phase transition a,t m — 1.3 as a phase 
transition in the baryon diffusion coefficient D. We are interested in the finite density effect on this 
transition. 

At finite density, we know that black hole embeddings capture the physics at all temperatures, i. e. 
the entire parameter regime of m. The fundamental phase transition in this case is a transition between 
two different black hole embeddings [SS] . As discussed in ref. [S^ , the baryon density affects the location 
and the presence of the fundamental phase transition. The transition is of first order only very close to 
the separation line between the regions of zero and non-zero baryon density shown in figure 13.11 Note 
that as discussed in refs. [50] [SHI [iS] there exists a region in the {d,T) phase diagram at small d and T 
where the embeddings are unstable. This instability disappears for large d. 

We study the baryon diffusion coefficient at different baryon densities. Figure shows the Diffusion 
coefficient Z? as a function of the ratio of quark mass to temperature m. By fixing the quark mass we 
may think of m as the inverse of the temperature. 
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Figure 5.2: The normalized baryon diffusion coefficient as a function of normalized inverse temperature. 
At densities below d — 0.00315 we observe a multivalued dependence on m, signaling a phase transition. 
For the behavior of the coefficient in a larger range of m see figure 14.11 

We find that the phase transition is slightly shifted towards smaller temperatures when we increase 
the density. At a critical density of d* = 0.00315 the phase transition temperature is given by m = 1.31. 
Beyond the critical density the transition vanishes, in agreement with the critical density d* for the phase 
transition in the quark condensate, discussed in ref. |58) . 

5.2 A new phase transition at finite isospin potential 

In this section we have a different look at the mesonic spectral functions, introduced in section [3] There 
we where interested in the behavior of the quasiparticle resonances and the modification of the particle 
spectrum when we leave the limit of zero temperature and vanishing particle density. Here, we focus on 
a new phenomenon occurring at high densities. 

We recall that the three solutions X, Y and to the fluctuation equations of motion (j3.82l) to p.84|) 
constitute the isospin triplet of mesons which may be constructed out of the isospin 1/2 quarks of the 
field theory. This is analog to the p-meson in QCD. We discovered that the mode coincides with the 
solution in case of a pure baryonic chemical potential, while the other two solutions have peaks in the 
spectral function at lower and higher values of to, c. f. figure [57TT] The magnitude of this splitting of the 
spectral lines is determined by the chemical potential and the undetermined coupling Ca- 

In the limit of zero frequency to 0, equations (j3.82p and (j3.83p coincide and will result in identical 
solutions X and Y . In this limit the solution E^ ^ though, differs from X and y, by means of the last 
term. So for small frequencies tr, we expect differences between the solutions E^ and X^Y. All three 
equations of motion depend on the particle density d parametrically, since the density has infiuence on 
the background fields. 

Spectral functions and quasi normal modes at high densities 

We work in the canonical ensemble and will now investigate the effects of variations in d. Spectral 
functions for various finite baryonic and isospin densities d are shown in figure 15.31 As in section 13.31 
the peaks in these spectral functions indicate that quarks form bound states. At low baryon densities 
we recognized the positions of the peaks to agree with the supersymmetric result (j3.67p . Increasing 
the quark density leads to a broadening of the peaks, which indicates decreasing stability of mesons at 
increasing baryon density [STl 1162] . At the same time the positions of the peaks change, which indicates 
a dependence of the meson mass on the baryon density. Now, further increasing the quark density leads 
to the formation of a new structure at to < 1. We will discuss this structure below together with the 
results at finite isospin density. 

We now turn to the effects of finite isospin density on the spectrum. The peaks in the spectral 
functions again correspond to mesons. An interesting feature at finite isospin chemical potential is the 
formation of a new peak in the spectral function in the regime of small tn at high density/high chemical 
potential, see figure [Ol Notice that compared to the baryonic case, the density at which the new peak 
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Figure 5.3: Spectral functions for various baryon densities (left) and isospin densities (right), again 
normalized to NfNcT'^ /A. At increasing densities d the peaks are smeared out, as we saw in the discussion 
of the spectral functions in section |31 At very high densities a new structure forms at small tr. 



forms is about two orders of magnitude smaller. As in the baryonic case, the excitations related to the 
supcrsymmetric spectrum broaden, the corresponding mesons become unstable. 

We pointed out that the structure of the spectral function is determined by the pole structure of the 
retarded correlator, see section 13.11 The poles of this function are located in the complex w-plane at 
positions r2„ G C. The spectral functions show the imaginary part of the correlator at real valued w. 
Any pole in the vicinity of the real axis will therefore introduce narrow peaks in the spectral function, 
while poles far from the real axis have less influence and merely introduce small and broad structures. 

In section [3.11 we outlined how the imaginary part of the quasinormal modes describes damping, as 
long as Im f2„ < 0. The short note on the pole structure demonstrated the dependence of the position 
of the quasinormal modes on the chemical potential/particle density. From figure we deduce that at 
higher densities than studied so far, a quasinormal mode approaches the origin of the complex lo plane 
as the particle density is increased. We observe a pole at ro = for a certain particle density dcrit, the 
value depends on m. An impression of the variation in the spectral function is given in figure [531 

In figure lSTSl we qualitatively sketch the result from the investigation of the behavior of the quasinormal 
modes closest to the origin of the complex m-plane. These modes do not produce the peaks corresponding 
to the spectrum (|3.67p . At low densities all quasinormal modes are located in the lower half plane. When 
increasing the isospin density, the lowest frequency modes of the solutions X and Y to (|3.82p and p.83p 
move towards the origin of the frequency plane. At the same time two quasinormal modes of move 
towards each other and merge on the negative imaginary axis, then travel along the axis towards the 
origin as one single pole. At the critical value of d = dcrit the modes from X and Y meet at the origin, the 
quasinormal modes from still reside in the lower half plane. This observation matches the discussion 
at the beginning of this section, where we expected X and Y to behave similarly at small ro, while E^ 
should differ from this behavior. Upon further increasing the isospin density, the modes from X and Y 
enter the upper half plane, maintaining their distinct directions. The sign change in Im n from Im 51 < 
to Im O > indicates that a damped resonance changes into a self-enhancing one, and thus introduces an 
instability to the system. Figure 1^751 illustrates the transition of a quasinormal mode of Y from the lower 
half plane to the upper half plane. The E^-mode does not enter the upper half plane at any value of d 
we considered. Compare this to the values of d in figure 15.31 at which the pole induces visible structures 
at small fo. A comparable movement of poles in a different but related setup was found in ref. |163) . 
There the quasinormal modes of correlation functions of electromagnetic currents were investigated as a 
function of temperature. 

In the following we interpret the observation of decaying mesons and the emergence of a new peak 
in the spectral function in terms of field theory quantities. In particular we speculate on a new phase in 
the phase diagram for fundamental matter in the D3/D7 setup. 

In the far UV, the field theory dual to our setup is supcrsymmetric, thus containing scalars as well 
as fermions, both of which contribute to the bound states we identified with mesons, even when super- 
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Figure 5.4: Plot of the spectral function for the 
mode Y around tr = 0. At a value of = 15.35 
a pole appears at the origin. This behavior is 
due to the movement of poles in the complex 
ro-plane, illustrated in figure [ 
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Figure 5.5: A sketch of the positions and move- 
ments of the quasinormal frequencies under 
changes of . Color indicates the function: 
red — F, green = X, blue ~ . The symbols 
indicate the range oi d^: o < dcrit, • = dcrit, 
■ > dcrit- Poles in the gray region introduce 
instabilities. 
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Figure 5.6: Contour plots of the spectral function for the mode Y around ir = in the complex ro-plane. 
The density increases from the left plot at sub-critical density to the right one at super-critical density. 
Here, the pole in the upper half plane introduces an instability. 
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Figure 5.7: the (/i^, T)-plane. In the blue shaded region D7-branes have the topology of Minkowski 
embeddings, the white and brown regions are modeled by black hole embeddings. These become unstable 
in the brown region. The boundary of the unstable region asymptotically seems to agree with the thin 
gray line of constant density = 20.5. 

symmetry is eventually broken. The meson decay at non-vanishing particle densities may be explained 
by the change of the shape of the potential for the scalars in the field theory upon the introduction of a 
non- vanishing density. As outlined in appendix |D1 a chemical potential may lead to an instability of the 
theory, since it induces a runaway potential for the scalar fields at small field values |164) . Nevertheless, 
interactions of (j)*-type lead to a Mexican hat style potential for larger field values. In this way the theory 
is stabilized at finite density d while the scalar fields condensate. This squark condensate presumably 
contributes to the vev of the scalar fiavor current, 



In the AdS/CFT context, the presence of an upside-down potential for the squark vev has been shown 
in ref. [54] using an instanton configuration in the dual supergravity background. 

The occurrence of a pole in the upper half plane of complex frequencies at finite dcrit indicates an 
instability of the theory. A comparable observation was made in ref. |81) , where in fact the vector meson 
becomes unstable by means of negative values for its mass beyond some critical chemical potential. 
The difference between this work and ref. [5T] is that our model includes scalar modes in addition to 
the fundamental fermions. Nevertheless, in both models an instability occurs at a critical value of the 
chemical potential. The theory may still be stabilized dynamically by vector condensation |165j . In this 
case the system would enter a new phase of condensed vectors at densities larger than dcrit, in accordance 
with the expectation from QCD calculations |166[ 1551 1167] . 

We perform the analysis of the pole structure at tn = for various m, and interpret the phenomenon 
of the transition of poles into the upper half plane at finite critical particle density as a sign of the 
transition to an unstable phase. We relate the critical particle density dcrit to the according chemical 
potential jl^ by = limp_>oo (p) ^^'^ ^^e the pairs of m and critical dimensionful fj,^ to trace the line 
of the phase transition in the phase diagram of fundamental matter in the D3/D7 setup. The result is 
drawn in figure [5771 The picture shows the (^^, T)-plane of the phase diagram and contains three regions, 
drawn as blue shaded, white, and brown shaded, as well as solid lines, separating the different regions. 

The blue shaded region marks the range of parameters, in which fundamental matter is described by 
D7-branes with Minkowski embeddings. The line, delimiting the blue region, marks the line of phase 
transitions to the black hole phase, where fundamental matter is described by D7-branes which have 
black hole embeddings. Using the symmetry of the DBI action, this phase transition line can be mapped 
to the line of phase transitions between Minkowski and black hole embeddings, present at finite baryon 
chemical potential [551 1501 HT] . 

The brown shaded region in the phase diagram in figure 15.71 marks the observation made in this 
section. The line delimiting the brown region marks the values of dcrit at which the pole in the spectral 
function appears at to = 0. Beyond this line we enter the brown shaded unstable region. 

We observe that the separation line of the unstable phase asymptotes to a straight line at high 
temperatures. Within the values computed by us, this line agrees with the asymptotic behavior of 




(5.1) 
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the contour of particle density with d « 20.5, drawn as a thin gray line in the phase diagram. We 
thus speculate on a finite critical particle density beyond which the black hole phase is unstable. This 
interpretation is supported by analogous studies of the phase diagram of A/" = 4 super- Yang-Mills theory 
with R-symmetry chemical potentials, where a similar line in the phase diagram was discovered |168II169] . 
The remaining question is whether the brown shaded phase in figure 15.71 indeed is unstable in the sense 
that it inaccessible for any physical setup, or if there is a way to stabilize the system in the parameter 
range of question. Recent publications revealed that the introduction of a further vev for a different 
gauge field component on the stack of probe branes leads to a stabilization of the system (1701 1171j . The 
resulting setup exhibits a second order phase transition to the new phase, which bears analogies to the 
theories of superfluidity and superconductivity [1701 11721 11731 1171j . 

Note that the location of the transition line to the unstable phase in figure 15.71 as well as the results 
shown in figure [5751 and figure [5?5] are obtained from the analysis of poles in the spectral functions. These 
functions in turn are obtained as solutions to equations (I3.82p to p.84p . which do depend on the so far 
unknown factor in determining the self coupling of the gauge field on the brane. The computation 
of this factor is left to future work. It will determine the exact position of the boundary of the brown 
shaded region in figure 13.11 This will answer the question whether there is a triple point in the phase 
diagram and if the color shaded regions meet at a common border. Moreover, other poles than the ones 
investigated here may have influence on the stability of this system. 

5.3 Summary 

We made two observations concerning the thermodynamic behavior of fundamental matter in the D3/D7 
setup. 

First, in our simple approximation of the baryon diffusion coefficient, we observe the fundamental 
phase transition and its dependence on the baryon density. We find that increasing the density from 
zero, where the transition temperature is given by m = 1.3, the transition temperature is lowered slightly 
until the transition vanishes at a critical value of d* = 0.00315, where the transition occurs at to = 1.31. 
This confirms the results of ref. |58j . where the transition was observed in the quark condensate. 

Second, we observe a new phase transition which renders the D3/D7 setup unstable at values above a 
critical isospin density. This becomes manifest by quasinormal modes of the fluctuations which develop 
positive imaginary parts in this region of the phase diagram. The exact position of the phase transition 
line cannot be determined yet. However, we speculate that the instability is due to a modification of 
the potential for the scalar fields in the field theory. This instability can be cured by vector meson 
condensation. 

It is tempting to compare figures 15.71 and 15.11 However, we point out that we cannot interpret the 
brown shaded region in fig. 15.71 as a direct analogon of the color-flavor-locked phase (CFL) in fig. 15.11 
since the parameter range scanned by us only allows to observe the phase transition to the new phase 
only at finite isospin chemical potential and not at finite baryon potential. Also the critical point at 
finite (T^,/i^) is not reproduced so far from the thermodynamics of the D3/D7 model. Nevertheless, 
the sheer appearance of the phase diagram in figure 15.71 may serve as a motivation for further efforts in 
exploring the phase diagrams of holographic models. 

6 Conclusion 

We considered different generalizations of the AdS/CFT correspondence in order to shed light on in- 
medium effects on the fundamental matter in holographic models for the quark-gluon plasma. The 
influence of the medium was parametrized by the values of the temperature and particle density. 

The first aspect we considered in |2]was the influence of the medium on bound states of quarks, in 
particular vector mesons. The description of mesons from first principles is interesting in its own right, 
because the strong coupling parameter forbids to apply well established perturbative methods in QCD. 
In the holographic setups, mesonic excitations arise more or less naturally as vibrational modes of open 
strings on D-branes. In the low energy limit, they account for fluctuations of supergravity fields. We 
presented the capabilities of a certain realization of a D3/D7 brane configuration by deriving the spectral 
functions for vector mesons from it. In the limit of vanishing temperature and density the derived spectra 
agree with the previously known results. The main achievement of our efforts, however, was the extension 
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of the spectral description of vector mesons into the finite density and temperature regime for all values 
of quark masses and temperature. We observe the melting of mesons at high temperature and at the 
same time studied the effects of finite particle density. Technically, we related the characteristics of the 
spectra to the behavior of the quasi normal modes of the excitations that holographically account for 
mesonic bound states of quarks. 

The main contribution to a better understanding of in-medium effects from this project is the deriva- 
tion of in-medium effects on the spectra. We observe a destabilization of mesonic bound states with 
increasing particle density in the quark-gluon plasma, which is simultaneously accompanied by a slight 
shift of the meson masses to higher energies. Due to the fact that the holographic models are too com- 
plex to be solved by analytical methods alone, the precise mechanisms that account for this behavior 
are difficult to reveal. A probable physical explanation for the destabilization certainly can be seen in 
the fact that in the strongly coupled medium the surrounding free quarks alter the binding interquark 
potential of the mesonic bound state. The closer a quark of the medium comes to a constituent quark 
of the meson under consideration, and the higher the amount of such perturbing spectator quarks is, 
the more influence can be expected from the medium on the mesons. Therefore increasing the baryon 
density (which can be seen as a measure for the amount of free quarks in the medium which are not 
bound into mesons) accounts for accumulating perturbation of the binding quark-antiquark potential, 
eventually leading to a dissociation of the meson. The shift in the meson mass may also be a consequence 
of the modification of the interquark potential, which in turn shifts the binding energies and therewith 
the energy content of a meson. 

Without speculating further on the mechanisms that lead to the observation we made, we note that 
our results are in qualitative agreement with phenomenological models and observation from experiment. 
The fact that our result is a non-trivial consequence derived from the D3/D7 setup can be seen as an affir- 
mative answer to the question whether string theory motivated models can capture phenomenologically 
relevant physics. 

Another such example was also derived in the context of meson spectra. Namely, we have shown 
that the introduction of finite isospin chemical potential indeed leads to a mass splitting of the different 
components of the isospin triplet, constituted by the three possible isospin one combinations of quark- 
antiquark pairs. While this is a success on the one hand side, the quantitative evaluation of the mass 
diff'erence remains as a task for future investigation, as the magnitude of the mass splitting heavily relies 
on meson coupling constants which are not determined yet and where chosen arbitrarily in our setup. 
The qualitative observation of the mass splitting, however, can be explained entirely analytically. We 
notice that the degenerate spectrum at vanishing isospin density stems from the fact that we have an 
SU(2) isospin symmetry in our system. By introducing a finite vaciium expectation value for one single 
generator, we break this symmetry and thereby suspend the degeneracy. From the equations of motion 
we can read off that the vev does not aff'ect the longitudinal component in flavor space but shifts the 
energy eigenvalues of the transverse modes by identical absolute amount with opposite sign. 

The second observable we studied was the diffusion coefficient of both baryon and isospin charge 
as well as the diffusion coefficient of mesons. The motivation to consider baryon and isospin diffusion 
apparently is to understand the transport processes of quarks and antiquarks in the QGP and quark 
matter as e. g. expected to exist in neutron stars. The interest in mesons stems from the fact that there 
is experimental evidence for mesons to survive the deconfinement transition to the QGP. We capture this 
eff'ect in our setup, as we have shown by observing discernible peaks in the mesonic spectral functions at 
finite temperature, discussed above. 

The results for the baryon diffusion coefficient where derived in an extremely simplified manner 
by plugging in the results for the embedding functions of the D7-brane in our background into the 
formula for the diffiision coefficient derived from the membrane paradigm. The main purpose of this 
task was to show that the D3/D7 setup at finite density is able to yield baryonic diffusion parameters 
for the plasma for ratios of the quark mass to temperature in both regimes, below and above the phase 
transition for fundamental matter. At vanishing density, the diffusivity of quarks normalized to the 
inverse temperature was known to be almost independent from the mass in the regime of light mesons 
(compared to the deconfinement temperature) and to be reduced monotonically with increasing quark 
mass. In our simple extension, we have shown that the effect of finite density on the normalized diffusion 
of baryon charge is leading to a dependence on the diffusion coefficient that exhibits a minimum for 
quarks with masses niq slightly heavier than the scale determined by the critical melting temperature 
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Tc as TO^'''* = ^ Y^ ^c- We identify the origin of this behavior as the fact that neither asymptotically 
heavy nor massless quarks will be influenced from the thermal momentum scale. A particle, however, 
with intermediate mass is certainly sensitive to momentum transfer by e. g. collisions with particles in 
the medium. 

Moreover, we observe that with increasing density, the mass dependence of the diffusivity becomes 
smaller. We address this to the fact that an extremely high density, accompanied by a likewise high 
chemical potential, outweighs the energy scale set by the finite temperature and in this way suppresses 
the intermediate mass scale dependence. 

Although we knew that our simple ansatz could not capture all effects of finite density in this way, 
our results where proven to be qualitatively correct in later publications. These, by the way, support 
the above comment on the rivaling scales of quark mass and chemical potential. The results of ref. [119] 
show that at large quark masses (which then outweigh the energy scale of the chemical potential) the 
diffusion constant indeed depends on the quark mass but is almost independent of the particle density. 

The situation at non-vanishing isospin chemical potential was analyzed by means of the dispersion 
relation for particles carrying isospin charge. Here, we restricted to massless quarks and small chemical 
potential. Within the tight restriction of our setup we where able to derive a frequency dependence 
of the diffusion coefficient, which can be interpreted as a dependence of the diffusion coefficient on the 
energy of the diffusion massless particles. 

The most extensively investigated transport coefficient however is the diffusion coefficient for heavy 
scalar mesons in the quark-gluon plasma. We set up a kinetic model that allowed for a derivation of 
the (inverse of the) diffusion coefficient at both strong and weak coupling. This enabled us to compare 
the perturbatively obtained weak coupling result to the holographic strong coupling result. Moreover, 
we where able to derive the polarizability of mesons from holographic duals. The latter results for the 
polarizabilities have to be read with care, as they rely heavily on the short-distance dipole approximation 
of the underlying effective model. As the short distance dipole interaction with the medium most likely 
rely heavily on large momentum transfer and thereby on weak coupling contributions, the validity of the 
AdS/CFT contributions may be vulnerable to serious criticism. Nevertheless, the long range effects on 
the momentum distribution of the scattered medium particles should be captured by our gauge/gravity 
model. We therefore divide out the effects due to polarizability, and parametrized by (5M)^ and compare 
the quotient of momentum broadening (inverse diffusion coefficient) k and (SM)'^. As a result, within 
the limits of the validity of our assumptions, we observe a reduction of the momentum transfer from the 
mesons to medium particles at strong coupling compared to weak coupling. This has consequences on 
the equilibration of the meson momentum distributions into thermal equilibrium, which we expect to be 
slowed down at strong coupling. Hopefully, measurements of the heavy meson momentum distributions 
in heavy ion collisions at RHIC and LHC will allow for a comparison of experimental data with our 
theoretical expectations. 

Finally, we devoted one to the analysis of the phase structure of fundamental matter in the holographic 
description of quark matter. 

One result was the observation of a phase transition in the baryon diffusion constant, which shows 
parallels to a previously observed phase transition in the the quark condensate, which vanishes at a 
critical baryon density d* . We observe a dependence of the value of the quark mass to temperature 
ration at which this phase transition occurs on the density. And we identify a critical density at which 
the phase transition vanishes. This density matches the value d* mentioned above, which makes us 
believe that we observe the same physical transition in just an other parameter. 

More important is the observation of a new phase transition at finite isospin density. Above we 
described the observation that the mass eigenvalues of two of the three components of the isospin triplet 
vector mesons experience shifts due to finite chemical potential. Increasing the isospin chemical potential 
n' , we observe an instability of our setup at a critical value fi' (T) of the chemical potential. This value 
depends on the temperature of the medium. By numerical evaluation of n'{T) we are able to trace out 
the boundary of the stable phase in the (fx' , T)-plane of the phase diagram. Recent publications indicate 
that the theory can be stabilized even beyond this line if additional gauge field components on the flavor 
branes acquire finite vacuum expectation values. This indicates that the boundary we traced out in the 
pase diagram marks the border between two different phases. The exact position, however, is subject to 
the same open questions we addressed when we discussed the splitting of the vector meson spectrum. 

As a general conclusion, we ascertain that the D3/D7 setup for holographic duals to strongly coupled 
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gauge theories, provides the capabiUty to describe a rich amount of phenomenology of the dual field 
theory. In this article we highlighted a small part of it. Alluding to the general motivation behind 
applications of AdS/CFT to bridge the gap between string theory and phenomenologically relevant field 
theories, we finally end with a satisfactory statement. The analyses and observations described in this 
work show that the gauge/gravity duality is not a one-way street. In one direction we where able to 
confirm many observations in holographic models by known results and expectations from established 
field theories such as QCD, and even experiment. In this way our confidence in the applicability of the 
correspondence to real world phenomena was strengthened. In the other direction, using the example 
of meson diffusion, we discovered ways to derive results in regimes of field theories, which so far where 
inaccessible, and hopefully bear at least qualitative truth when compared to field theory results or 
experiments in future. 
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A Coordinates for the AdS black hole background 

Numerous local coordinate systems arc used to parameterize the AdSs x S'^Schwarzschild black hole 
background. The problem at hand determines which of them is most useful. Here we list some of the 
common coordinates and the transformations between them. For all of the following coordinate systems 
we use the same symbols 

= AngsNa'^, r^ = TttR^. 



i?2 1 



ds2 = _ (-/(r) dt^ + dx2) + ——dr^ + R^dQl (A.l) 



Coordinate system 1 

^2 

i, x^, x^, x'^ € ]R, r > 0, dfig — metric of the unit 5-sphere 



with 



and 



horizon at r — ro 

boundary at r — > oo 



Coordinate system 2 

Introduction of a new radial coordinate 

= r"^ + \/r'^ ~ 

transforms (jA.ip into 

d.s^ = ^ (-^ dt^ + /(,) dx^) + ^ {dg' + gW,) (A.2) 
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with 

Q Q 

horizon at Q = ro 

boundary at g ^ oo 

Parametrization of the radial part 

We can identify the part in the last pair of parentheses of (|A.2I) as nothing else than and we write it 
as 

6 

i~{ ' — ^ ^ ' ' 

iR*(ei,...,4) lR=^(e5,6) 

with radial coordinate g = QiY'^- We write this space as a product space of a four-dimensional 
IR'* in polar coordinates with radial coordinate w > and a two-dimensional with radial coordinate 
i > 0, such that g^ = uP' + L"^ . The subspace parametrized by (w, L) is the first quadrant of a Cartesian 
coordinate system and can also be parametrized by its radial part g and an angle < < 7r/2, 

L ~ gcosd, 
w — gsinO 

such that 

clL^ + dw^ =dg^ + g^de^. 
Finally, we introduce x = cos 6 and thus can write (jA.2[) as 

d.^ = ^ (-^ dt^ Kg) dx^) + i?2 + (1 - x') dill + (1 - xT" dx' + x^d^^) . iQ) 

Coordinate system 3 

Introduction of a new radial coordinate 



transforms (lA.ll) into 



with 



= ( - 4^ + /(^,) dx^) + {dv' + v^dnl) (A.3) 



4 4 



horizon at v — ^ 

boundary at w — > cxi 



Coordinate system 4 

Introduction of a new dimensionless radial coordinate 

„2 



transforms (|A.1I) into 



u = 



„2 n2 



ds2 = -2- (-/(u) dt2 + dx2) + ^ d^^ + i?2df]2, (A.4) 



where 

f{u) = l-u^, 

horizon at m = 1, 

boundary at u = 0. 
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Coordinate system 5 

Introduction of a new radial coordinate 



transforms (lA.ll) into 



with 



horizon at z = Zo 

boundary at z — 



ds^ ^ ^ (-f^z) + dx2 + -i- dz2 



/(^) 
,4 



(A.5) 



B Isospin diffusion related equations 
B.l Solutions to equations of motion 

Here we exphcitly write down the component functions used to construct the solutions to the equations 
of motion for the gauge field fluctuations up to order tn and q^. The functions themselves are then 
composed as in (I4.40p . 

The solutions for the components with flavor index a = 3 where obtained in ref. |22) . 

B.1.1 Solutions for Xa, and A^^ 



The function Xa{u) solves (|4.32l) with the upper sign and is constructed as in (|4.40l) from the following 
component functions, 



xo m 



OH 











Fl/2 


"^V 2" 








- r— 








+ 31n(l - 



2 ' 

TT^ - 9 In^ 2 + 3 ln(l - u) (In 16-4 ln(l + u)) 
+- u) (ln(4(l + u)) - 41nu' 

12(Li2(l-M)+Li2(-w)+Li2 



Gi = 



C 
2 



12 



1 + M 



+ lnuln(l + u) + Li2(l - u) + Li2(-M) 



(B.l) 

(B.2) 
(B.3) 



(B.4) 

(B.5) 
(B.6) 



where the constant C can be expressed it in terms of the field's boundary value X^'^^ = limK_j.o X{u, k), 

c = x^'^y 

^2 ln2 2^ 



X ( 1 + a/— ln2 + mro I — 



^q2 + 0U3/2 n4 



(B.7) 



The solutions of the equations of motion (j4.32l) with lower sign for the functions Xa (u) are given by 
P= -^J^+Oiu;), (B.8) 
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= c, 






~ rm 


Fl/2 


= iCJ- 


2 










-^12 





TT^ - 9 In^ 2 + 3 ln(l - u) (In 16 - 4 ln(l + m)) 
+ 3 ln(l + u) (ln(4(l + u)) - 41nu) 
- 12 ( Li2(l -u)+ Li2(-M) + Li2 



1 + u 



— +lnuln(l + u) +Li2(l - m) + Li2(-u) 



(B.9) 
(B.IO) 

(B.ll) 

(B.12) 
(B.13) 



with C given by 
C = X^'^^ 



X ( 1 - t\l —^ln2 - mw I — 



In 2 



^q2 + 0(t^,3/^q^)] , (B.14) 



that lini„^o X{u, k) = X'^'^^ . 

The solution for Aj^ solves (14.25^ up to order ro and with boundary value 



vbdy 



It is 



8-4itr In2 + 7r2q2 

.113 , 1 + M 
1 + ^2^^^ 

+ y f — + lnuln(l + u) + Li2(l - u) + Li2(-u) 



(B.15) 



B.1.2 Solutions for X'^, X' and A^' 



Here we state the solutions to (j4.29l) . This formula describes three equations, differing in the choice of 
= 1,2, 3. The cases a = 1, 2 give coupled equations which are decoupled by transformation from AJ' ^ 
to Xo and Xq. The choice a = 3 gives a single equation. 

The function X'^ is solution to (|4.33l) with upper sign. We specify the component functions as 



Fo 


= G, 






_ fm _ 


Fl/2 




- In 


Fi 


m 

- ~^12 












+ 12 (hhil 






u 








2u 



7r'^+31n^2 + 31n^(l + u) + 61n21n " 



1 + u 



1-u 



(B.16) 

(B.17) 
(B.18) 

(B.19) 

(B.20) 

(B.21) 
(B.22) 
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where the constant C can be expressed in terms of the field's boundary value X^'^'^ ~ lini„_j.o X{u, k), 



C 



\/2mo + mro In 2 + 
To get the function Xq, we solve (|4.33p with the lower sign and obtain 



(B.23) 



(B.24) 



^^i/2 = .CW-ln^ 



TT^ + 3 In^ 2 + 3 In^ (1 + u) + 6 In 2 In 



Gi = Cln 



+ 12 ( Li2(l - u) + Li2(-u) - Li2 

1 + u 
2u 



l + u 
l-u 



(B.25) 
(B.26) 

(B.27) 

(B.28) 

(B.29) 
(B.30) 



where the constant C can be expressed it in terms of the field's boundary value X^'^^ = lim^^o X{u, k), 



C 



q'X^^^y + yoqX^/^ 



iVSrnro + mro In 2 — 
The solution for (j4.29p with a = 3 is the function Aq' , given by 



iKi — q'^ 



1 + — In ■ 

2 l + u 



q2 In 



l + u 
2u 



B.1.3 Solutions for X!^, X'^ and A^' 
We give the derivatives of X3 and X3 as 



"J + m ~, 



The solution for A^ is 



^0 • 



(B.31) 



(B.32) 



(B.33) 
(B.34) 

(B.35) 



B.2 Abelian Correlators 

For reference we quote here the correlation functions of the Abelian super-Maxwell theory found in 
ref. |22j . The authors start from a 5-dimensional supergravity action and not from a Dirac-Born-Infeld 
action as we do. Therefore there is generally a difference by a factor N /A. Note also that here all Nf 
flavors contribute equally. In our notation 



Gab ^ah 
11 — ^22 — ^ 



16^? 



°° 167r(^CJ - Dq^) 



(B.36) 
(B.37) 
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where D = l/(27rr) 



B.3 Correlation functions 

In this section we hst the real and imaginary parts of the flavor currents in the first two flavor- 
directions a — 1,2 and in the third Lorentz-direction coupling to the supergravity- fields and 
(as defined in ^^). 



Im G^^(uj > = - Im G^Juj < = / ,^ ' \ ^' J/ , B.41 

- ' 33y J Svr 2^ w 27rT 21 ' ^ ^ 



ReG33(w < 0) = ReG^^iuj > 0) = . ^""^ — j , (B.42) 

ImG33(cj < 0) = 0,IniG'33(w > 0) 0. (B.43) 

C Coupling constant for vector meson interaction 

In this section we show how the coupling constant for the interaction of vector mesons can be computed, 
extending the ideas presented in |43J. This coupling constant in the effective four-dimensional meson 
theory can be determined by redefinition of the gauge fields such that the kinetic term has canonical 
form. This coupling constant depends on the geometry of the extra dimensions. 

(2) 

First we consider the eight-dimensional theory determined by the DEI action 6*13131 expanded to 
second order in the fluctuations A, 



l\2 



(2) _ r7(2W) 



J d^^v^ g^^'^g^^F^^Fp^ , (c.i) 



where g contains the background fields and we simplify the analysis by considering only Abelian gauge 
fields. Defining the dimensionless coordinate p — gj R and integrating out the contribution of the S'^, we 
obtain 



4 

To obtain a four-dimensional effective theory we have to integrate over the coordinate p. This contribution 
depends on the geometry induced by the p dependence of the metric factors. However, we expect that 
it is independent of the 't Hooft coupling A. We parametrize this contribution by c\. The kinetic term 
of the effective theory is then given by 



4 

where the prefactor may be written as 

r7(27ra')2vol(S3)i?*c', 



(C.4) 



where the numerical values independent of the 't Hooft coupling are grouped into the coefficient c^- 
From this we can read off that a rescaling of the form 

i^^i (C.5) 

Va 



casts the Lagrangian into canonical form with a prefactor of l/g. 



2 

YM • 
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D Chemical potentials in field theories: Runaway potential and 
Bose-Einstein condensation 



In our setup we consider a field theory which is supersymmetric in the far UV. Its fundamental matter 
consists of complex scalars (squarks) and fermionic fields (quarks). In this section we describe the effect 
of the chemical potential on the field theory Lagrangian and on the vacuum as e. g. in ref. |164] . We 
consider a theory with one complex scalar (j) and one fermionic field ip with the same mass niq coupled to 
an U(l) gauge field A^. The time component of the U(l) gauge field has a non-zero vev which induces 
the chemical potential ^, 

A^=fiS^o. (D.l) 

The Lagrangian is given by 

L = - {D^cj,)* D^4> - ml4>*4> - i>{Ip+ m^)^ - -Fp.vF'''' , (D.2) 

where — ~ iA^ is the covariant derivative and = d^A^ — d„Af^ the field strength tensor. 
Expanding the Lagrangian around the non-zero vev of the gauge field, it becomes 

L = -df,4>*d^'(t> - {ml - i?)(t>*(t> + MJ(f - V^(^ + 'mq)'4i + ^lJ^ (D.3) 

where = i - (9^0)), and (J^)^ = -i^l^j'^il: are conserved currents. These conserved 

currents are the population densities Ns for the scalar field and Np for the fermionic field, such that the 
linear terms in the Lagrangian are uNs and ^Np. 

The mass term — (m^ — ij?)4>^ of the Lagrangian (jP.Sp introduces an instability if /i > mq since 
the corresponding potential V = (m^ — + • • • is not bounded from below. In some systems this 

runaway potential is stabilized by higher interactions and becomes a Mexican hat potential such that the 
scalar condenses and the scalar density becomes non-zero. This condensation is known as Bose-Einstein 
condensation (BEC). 
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